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ABSTRACT 


A  compilation  has  been  made  of  numerous  topics  in  the  theory  of 
compressible-fluid  dynamics.  The  presentation  is  designed  as  a  guide 
for  students  as  well  as  a  source  of  reference  material  for  the  research 
worker.  The  choice  of  material  has  been  influenced  by  the  needs  and  pro 
jects  at  the  Los  Alamos  Scientific  Laboratory,  where  much  of  the  theoret 
ical  work  in  fluid  dynamics  is  carried  out  by  numerical  techniques  on 
high-speed  electronic  computers.  There  is,  however,  no  discussion  of 
the  numerical  methods  themselves. 
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INTRODUCTION 

The  major  emphasis  in  fluid-dynamics  studies  has  gradually 
shifted  over  the  last  few  years.  During  the  early  days  of  aircraft 
almost  all  processes  of  interest  took  place  at  speeds  far  less  than 
that  of  sound.  Fluids  could  therefore  be  treated  as  incompressible,  or, 
as  having  a  density  which  varied  with  time  or  space  only,  generally  in 
a  predetermined  manner.  As  aircraft  speeds  increased,  other  high-speed 
phenomena  became  of  interest,  and  assumption  of  an  incompressible  fluid 
was  no  longer  valid  for  many  practical  situations.  As  a  result,  new 
analytical  techniques  were  developed  and  applied  to  a  variety  of  design 
problems,  especially  in  aerodynamics. 

Meanwhile,  interest  has  grown  in  the  dynamics  of  fluids  adjacent 
to  high- intensity  explosions,  and  in  the  theoretical  techniques  neces¬ 
sary  for  understanding  the  processes  involved.  The  complexity  of  the 
analysis  has  often  required  a  change  in  the  theoretical  approach  from 
analytical  to  numerical.  At  Los  Alamos,  various  numerical  methods  have 
been  evolved  and  applied  by  means  of  high-speed  electronic  computers. 

As  a  background  for  this  work,  it  has  been  necessary  to  gather  from 
numerous  parts  of  the  literature  analytical  techniques  by  which  solu¬ 
tions  could  be  obtained  for  comparison  with  the  numerical  results. 


-7- 


INTRODUCTION 


The  analytical  methods  sometimes  required  extension,  modification,  or 
reinterpretation,  or  new  techniques  had  to  be  developed.  It  is  the  pur¬ 
pose  of  this  manual  to  present  this  analytical  material  in  the  hope  that 
it  will  be  a  useful  background  to  others  whose  compressible-fluid  pro¬ 
blems  are  similar  to  ours. 

This  compilation  is  for  reference  as  well  as  teaching.  It  thus 
contains  many  detailed  and  alternative  formulations  which  may  at  first 
confuse  the  beginner;  but  at  the  same  time  it  omits  certain  types  of 
elegant  argument  for  which  referaxce  to  original  papers  is  more  appro¬ 
priate.  The  student  interested  in  further  details  can  profit  from  sev- 

* 

eral  excellent  books;  these  have  often  served  as  basic  references  herein. 

It  is  suggested  that  the  student  read  quickly  the  first  four  chap¬ 
ters  of  this  manual  and  any  subsequent  sections  of  special  interest, 
and  then  return  to  the  earlier  chapters  for  detailed  study  or  reference 
as  the  need  arises.  The  newcomer  to  fluid  dynamics  will  find  it  useful 
to  have  a  qualitative  feeling  for  the  phenomena  involved,  before  pro¬ 
ceeding  to  the  details  of  formulation.  The  rest  of  this  Introduction 
offers  such  a  qualitative  picture . 

* 

These  include: 

Richard  Courant  and  K.  D.  Friedrichs,  "Supersonic  Flow  and  Shock 
Waves,"  Interscience  Publishing  Company,  New  York,  1948. 

R.  von  Mises,  ’Mathematical  Theory  of  Compressible  Fluid  Flow," 
Academic  Press,  Inc.,  New  York,  1958. 

L.  D.  Landau  and  E.  M.  Lifshitz,  "Fluid  Mechanics,"  Pergamon  Press, 
Ltd.,  Iondon;  Add! son-We sley  Publishing  Company,  Inc.,  Reading, 

Mass;  1959* 
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The  first  thing  that  is  needed  in  a  theoretical  treatment  of  fluid 
dynamics  is  a  means  for  characterizing  the  fluid.  Thus*  for  example,  if 
we  have  material  in  a  hox  and  wish  to  examine  the  dynamics  of  it,  then 
it  Is  necessary  to  specify  certain  field  variables  and  try  to  determine 
the  manner  in  which  they  change  with  time.  Such  field  variables  include 
mass  density,  velocity  of  each  part  of  the  fluid,  pressure,  internal 
energy,  temperature,  and  various  other  similar  quantities.  We  shall 
attempt  to  define  these  quantities  in  the  following  chapters  but  will 
here  proceed  on  the  basis  of  our  intuition  of  their  meaning. 

To  describe  how  these  quantities  change  with  time  from  a  given  ini¬ 
tial  configuration,  subject  to  prescribed  conditions  at  the  boundary  of 
the  fluid,  we  must  derive  a  set  of  governing  equations.  One  of  these 
will  be  a  statement  of  the  characteristics  of  the  fluid  itself,  the 
"equation  of  state"  which  describes  in  mathematical  terms  the  relation¬ 
ship,  at  any  given  instance  and  place,  among  such  quantities  as  tempera¬ 
ture,  density,  and  pressure.  If  the  material  is  a  simple  gas,  then  the 
form  of  the  equation  may  be  quite  simple,  but  if  the  material  is  anything 
else,  then  the  relationship  may  be  extremely  complicated  or  even  impos¬ 
sible  to  write  in  analytical  form. 

Concerning  the  dynamics  of  the  gas  itself,  we  shall  show  that  all 
processes  take  place  according  to  the  dictates  of  the  three  fundamental 
conservation  laws  of  physics,  namely,  those  of  mass,  momentum,  and  energy. 

It  is  useful  to  examine  the  motion  of  a  fluid  from  at  least  two 
-PnnriflTnfmt.Rl  iy  different  points  of  view.  In  one  case  the  coordinate 
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system  will  "be  fixed  relative  to  the  observer,  while  in  the  other  the 
coordinate  system  is  considered  to  follow  the  motion  of  the  fluid. 

The  former  viewpoint  is  called  Eulerian;  the  latter  is  called  Lagrangian 
The  nature  of  the  problem  will  dictate  which  form  of  the  equations  will 
be  more  desirable. 

For  many  applications  it  is  reasonable  to  neglect  the  effects  of 
heat  conduction  and  viscosity.  As  a  result  there  are  two  rather  remark¬ 
able  consequences.  The  first  of  these  is  that  the  entropy  of  each  ele¬ 
ment  of  the  fluid  remains  forever  constant.  The  second  is  that  the 
equations  may  lead  to  discontinuous  solutions  even  if  the  initial  and 
boundary  conditions  are  apparently  very  well  behaved.  The  constancy  of 
entropy  is  discussed  in  some  detail  in  Chapter  II.  The  treatment  of  dis 
continuity  will  receive  much  attention  in  many  parts  of  the  manual. 

The  mathematical  discontinuities  are  called  "shocks."  Their  be¬ 
havior  is  not  governed  by  the  differential  equations,  which  become 
meaningless  at  a  discontinuity;  but  properly  treated,  the  mathemati¬ 
cally  discontinuous  shock  is  useful  because  it  closely  represents  a 
true  phenomenon  in  which  such  quantities  as  density,  pressure,  and  ve¬ 
locity  may  change  very  rapidly  over  a  very  small  distance. 

The  equations  that  are  derived  in  Chapter  I  show  that  weak  signals 
are  propagated  through  a  fluid  with  a  certain  characteristic  velocity 
called  the  "sound  speed."  The  sound  speed  will  be  shown  to  vary  with 
temperature  and  density  in  a  manner  which  depends  upon  the  equation  of 

state  of  the  fluid.  Generally,  with  an  increase  in  temperature  and  den- 
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sity,  the  sound  speed  goes  up.  Thus,  if  a  compressive  signal  is  formed 
by  a  piston  pushing  into  a  gas,  the  gas  near  the  piston  is  compressed 
and  heated  so  that  the  sound  speed  is  increased  above  that  which  is 
ahead  of  the  disturbance.  As  the  piston  velocity  increases,  sound  sig¬ 
nals  begin  to  pile  up  on  one  another,  and  the  front  of  the  disturbance 
wave  becomes  sharper  and  sharper.  This  produces  a  shock.  If,  on  the 
other  hand,  the  piston  is  withdrawn  from  the  gas,  then  the  sound  speed 
adjacent  to  the  piston  goes  down  because  of  the  expansive  temperature 
decrease,  and  signals  from  the  piston  are  not  able  to  overtake  the  dis¬ 
turbance  front  moving  into  the  fluid.  The  piston  motion  forms  a  rare¬ 
faction  and  no  discontinuity  occurs  in  the  fluid.  We  shall  later 
demonstrate  more  precisely  the  fact  that  shocks  are  formed  only  at  the 
fronts  of  compressive  waves. 

The  two  basic  elements  of  fluid  flow,  shock  and  rarefaction,  will 
prove  to  be  the  basic  building  stones  for  many  useful  analytical  solu¬ 
tions.  We  shall  see,  however,  that  it  is  seldom  possible  to  solve  the 
equations  simply,  because  they  are  nonlinear  partial  differential  equa¬ 
tions  which  lack  the  powerful  solution  techniques  available  for  linear 
equations.  Instead  we  shall  often  have  to  generate  solutions  to  prob¬ 
lems  by  a  careful  piecing  together  of  known  parts  of  solutions;  many 
of  the  examples  in  this  manual  have  been  formed  in  that  way. 

Almost  ail  theoretical  treatments  of  fluid  dynamics  are  based  on 
the  model  of  a  continuous  fluid  whose  elements  remain  forever  contiguous. 

In  reality,  a  fluid  is  formed  of  numerous  molecules  which  may  rapidly 
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move  about  so  that  neighboring  particles  at  one  instant  may  be  widely 
separated  a  short  time  later.  It  may  thus  seem  remarkable  that  the  sim¬ 
plified.  model  gives  useful  results  at  all.  In  Chapter  VIII  we  have  indi¬ 
cated  hov  the  success  of  the  continuous  model  can  be  explained,  and  also 
have  shown  how  the  failings  of  the  continuous  model  can  be  remedied  some¬ 
what  by  inclusion  of  the  effects  of  viscosity  and  heat  conduction.  In 
the  continuous  model,  viscosity  is  interpreted  as  a  drag  stress,  whereas 
in  reality  it  arises  from  diffusion  across  velocity  gradients.  The  equa¬ 
tions  with  viscosity  and  heat  conduction  do  not  tend  to  give  discontin¬ 
uous  solutions;  indeed  the  detailed  structure  of  a  shock  can  be  derived, 
and  the  results  agree  closely  with  experimental  measurements. 

These  qualitative  features  are  illustrated  in  this  manual  by  giving 
solutions  to  numerous  problems.  In  this  way,  the  reader  learns  that 
some  situations  in  fluid  dynamics  are  easily  analysed;  he  will  also  quickly 
realize,  as  he  tries  additional  problems,  that  there  are  many  apparently 
simple  situations  for  which  analytical  solutions  are  presently  impossible. 

The  progress  of  this  manual  has  profited  considerably  from  a  careful 
study  of  the  entire  manuscript  by  Cyril  Hirt,  from  a  detailed  checking  of 
several  sections  by  Jacob  Fromm,  and  from  numerous  discussions  with  my  as¬ 
sociates.  The  hard  jobs  of  editing  the  raw  manuscript,  removing  the  numer¬ 
ous  inconsistencies  in  nomenclature,  adding  a  coherence  of  style  among  the 
chapters,  and  tending  to  numerous  details  leading  to  the  final  form,  all 
with  exceptional  care  beyond  the  call  of  duty  were  accomplished  by  Marion 
Richardson.  Readability  of  the  final  form  owes  much  to  the  patient  and 
careful  typing  of  Grace  Cole. 
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CHAPTER  I 

THE  BASIC  EQUATIONS 


Introduction 


By  popular  concept,  a  fluid  is  a  mobile  substance  into  which  a 
solid  object  can  be  immersed  with  ease  but  which  does  not  sustain  the 
deformation  when  the  object  is  removed.  More  precisely,  a  fluid  may  be 
defined  as  a  material  for  which  the  stress  produced  by  a  shearing  strain 
is  independent  of  the  strain.  This  definition  does  not  exclude  those 
materials  for  which  the  stress  depends  upon  the  rate  of  shearing  strain; 
such  fluids  are  called  viscous. 

No  precise  description  of  the  dynamics  of  any  macroscopic  section 
of  fluid  is  ever  likely  to  be  possible;  such  a  description  would  involve 
the  exact  history  of  every  molecule.  There  are,  however,  numerous  situa¬ 
tions  in  which  theoretical  results  of  high  accuracy  can  be  obtained  with¬ 
out  such  detailed  knowledge. 

The  first  basic  requirement  for  any  theoretical  study  of  fluid  dyna¬ 
mics  is  a  means  for  characterizing  the  macroscopic  features  of  the  fluid. 
The  second  requirement  is  a  formulation  of  the  dynamical  problem  itself. 
In  this  chapter,  we  describe  these  for  rather  idealized  fluids,  namely, 
for  those  in  which  the  effects  of  viscosity  and  heat  conduction  are 
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negligible,  and  for  which  there  is  no  internal  source  of  energy  (as  would 
arise,  for  example,  from  chemical  reactions).  We  shall  call  such  a 
fluid  a  "simple  fluid;"  likewise,  we  shall  often  refer  to  a  "simple  gas." 
In  spite  of  these  idealizations,  the  resulting  equations  will  be  very 
useful  for  describing  quite  accurately  numerous  problems  of  dynamics. 

We  shall  be  limited  to  those  situations  in  which  the  fluid  can  be 
considered  as  macroscopically  continuous.  A  rigorous  discussion  of  the 
circumstances  under  which  this  is  possible  is  beyond  the  scope  of  this 
manual;  however,  we  can  formulate  the  conditions  loosely  as  follows. 
Consider  that  the  space  occupied  by  the  fluid  is  divided  into  numerous 
subregions.  For  each  subregion,  the  instantaneous  values  of  various 
appropriate  quantities  are  averaged  over  all  the  contained  molecules. 

Then  the  fluid  is  continuous,  for  our  purposes,  provided  that  for  all 
quantities  of  interest  a  set  of  subregions  can  be  found  such  that  the 
inter- subregion  variations  of  the  averages  are,  at  every  instant  of 
time,  negligibly  small  compared  with  the  averages  themselves.  This  re¬ 
quirement  implies  the  following  features: 

1 .  Each  subregion  is  large  enough  to  contain  a  large  number  of 
molecules . 

2.  Each  subregion  is  small  compared  with  any  macroscopic  feature 
of  the  gross  configuration. 

With  the  continuity  requirement  satisfied,  the  true  fluid  can  now 
be  replaced  by  an  appropriate  continuous  model.  This  substitution  can 
be  accomplished  by  the  use  of  some  suitable  interpolation  method  applied 
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to  the  subregion  averages;  the  very  nature  of  the  continuity  requirement 
shows  that  it  does  not  matter  what  reasonable  method  is  used.  In  prac¬ 
tice,  the  replacement  is  not  actually  performed,  but  the  concept  of  it 
is  useful  in  the  interpretation  of  results. 

We  now  assume  that  there  is  no  difficulty  in  talking  about  the 
value  of  a  fluid  property  "at  a  point,"  and  that  it  is  clear  what  is 
meant,  for  example,  by  differential  elements  of  mass,  energy,  and  momen¬ 
tum,  and  gradients  of  density. 

The  derivations  to  follow  are  based  upon  two  additional  assumptions. 
First,  the  molecules  which  are  in  a  particular  subregion  at  any  given 
time  remain  forever  in  that  same  subregion;  second,  except  for  macro¬ 
scopic  slippage,  every  pair  of  contiguous  subregions  remains  forever 
contiguous.  Actually,  these  assumptions  applied  to  a  gas  include  a  re¬ 
statement  of  the  neglect  of  viscosity  and  heat  conduction. 

We  shall  not  dwell  on  these  awkward  matters,  because  it  is  possible 
to  follow  quite  a  different  and  much  more  satisfactory  approach  to  the 
formulation  of  the  macroscopic  dynamics.  (An  introduction  to  the  alter¬ 
native  is  given  in  Chapter  VIII.)  We  shall  here  proceed  on  a  "semi- 
intuitive"  basis,  using  a  liberal  interpretation  of  our  basic  assump¬ 
tions  as  required. 

Characterization  of  the  Fluid 

To  completely  specify  the  state  of  a  fluid,  it  is  necessary  to 
give  the  values  of  certain  "field  variables"  for  every  point  in  the 
fluid.  Some  of  these  are: 
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"r  =  Lagrangian  coordinate.  This  is  the  permanent  label  which  is 
attached  to  a  point  moving  with  the  fluid.  It  may  be  chosen  in  any  man¬ 
ner  as  long  as  it  is  different  for  each  point.  The  symbol  is  as  shown, 
because  frequently  the  Lagrangian  coordinate  is  chosen  to  be  the  position 
of  the  point  at  some  reference  time  (t  =  tQ) •  Its  Cartesian  components 

are  then  x  ,  y  ,  z  .  (in  many  of  our  discussions,  the  subscripts  will 
o  o  o 

be  omitted  when  it  is  clear  that  the  Lagrangian  form  is  used.) 

t =  Eulerian  coordinate.  This  is  the  position  of  a  point  rela¬ 
tive  to  some  external  reference  frame.  The  Eulerian  coordinate  is  a 

function  of  ?  and  of  the  time,  t. 
o 

Each  of  the  following  field  variables  can  be  considered  to  de¬ 
pend  upon  the  time  and  either  the  Lagrangian  or  the  Eulerian  coordinate: 
p  s  Mass  density  of  the  fluid 
V  a  Specific  volume  (V  =  i/p) 

p  a  Pressure  r  ^  1 

&&  ,t) 

u  a  Velocity  of  the  point  u  =  - ^ - J .  The  velocity  could 

also  be  called  the  specific  momentum,  just  as  l/2  u*u* is  the  specific 
kinetic  energy . 

I  a  Specific  internal  energy 
T  =  Temperature 
S  a  Specific  entropy 

E  a  Specific  total  energy  (E  =  I  +  l/2  u*u* ) 

This  notation  is  used  throughout  the  manual  except  as  otherwise  explicitly 

stated.  Additional  terms  are  introduced  and  defined  as  they  are  needed. 

♦"Specific"  quantities  are  measured  per  unit  mass* 
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One -Dimen sional  Flow 

It  is  necessary  to  establish  clearly  the  difference  between  two 
types  of  "one-dimensional"  motion  of  a  fluid.  In  pseudo  one -dimensional 
motion,  all  quantities  characterizing  the  fluid  are  functions  of  one 
space  coordinate  only  (as,  for  example,  x  in  Cartesian  coordinates  or 
the  angle  6  in  cylindrical  coordinates).  In  true  one -dimensional  motion 
there  is  the  additional  feature  that  all  velocity  components  vanish 
except  the  one  along  that  coordinate  direction. 


Conservation  of  Mass 

Consider  a  closed  surface  composed  of  points  which  move  with  the 
fluid;  according  to  the  assumptions  of  our  model,  the  mass  of  fluid 
within  that  surface  will  never  change.  This  fact  can  be  exploited  in 
the  derivation  of  one  of  the  basic  fluid-dynamic  equations. 

As  an  example  of  the  derivation,  consider  a  true  one-dimensional 
flow.  The  mass  per  unit  transverse  area  contained  between  the  Eulerian 
boundaries  x^  auid  Xg  is 


Thus 


This  must  vanish  for  any  arbitrary  interval,  so  that  the  integrand  it¬ 
self  must  vanish  at  every  point: 


dp 


+ 


dpu 
”5x  “ 


0 


(1) 
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This  is  the  one-dimensional  mass  conservation  equation  in  Eulerian 
coordinates. 

In  three  dimensions,  consider  the  closed  surface  to  be  a  cube 
at  some  instant  of  time.  Then  the  mass  within  the  cube  is 


nX2  /2  PZ2 


p  dxdydz 


/X  vjy 

1  y1  1 


dm  / 
dt  ~  Jx 


<*2  Py2  nZ  2  n  r 

,  X,  I,  I toayaz  + 1, 


'2  ^2 


f  P  || 

y-i  L  atjZi 


y-,  zt  xi  zi  xi 


Ji  i 


dzdx 


"  i,2^2^,2  ^  +  S  (pux>  +  i  (ptty)  +  S 

Again,  the  integrand  must  vanish  at  every  point,  so  that  the  equation  may 
be  written,  in  vector  form, 

+  V-(  pu)  -  0  (2) 

The  above  procedure  for  forming  the  time  derivative  of  a  quantity 
is  quite  general.  Thus,  if  X?,t)  is  any  quantity  measured  per  unit  vol¬ 
ume,  then  the  total  of  that  quantity  throughout  some  volume  moving  with 

the  fluid  is 


A  =  /  Tv  dr 
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where  dr  is  a  generalized  volume  element,  and  the  range  of  integration 
covers  the  volume.  Then  one  can  derive,  analogous  to  above 


dA 

dt 


(3) 


This  formula  can  also  be  interpreted  by  rewriting  it  in  the  form 


IS  “Si  ft  **  (u) 

where  use  has  been  made  of  the  divergence  theorem  to  transform  the 
second  term  into  an  integral  over  the  surface;  n  is  a  unit  outward 
normal  to  the  surface  and  da  is  an  element  of  its  area. 

Thus,  the  total  rate  of  change  of  A  can  be  considered  to  arise 
from  two  processes.  The  first  term  on  the  right  of  (h)  signifies  the 
change  that  would  occur  if  the  surface  remained  fixed  —  did  not  move 
with  the  fluid.  In  the  second  term,  n*u*  da  is  the  rate  at  which  volume 
is  swept  up  by  the  moving  surface;  multiplication  by  A  gives  the  rate 
at  which  the  quantity  is  swept  up.  Thus,  the  second  term  "catches  up" 
with  what  was  lost  from  the  first  term,  in  which  the  surface  was  con¬ 
sidered  fixed. 


Conservation  of  Momentum 

The  momentum  per  unit  volume  is  pu*,  so  that  the  total  momentum  in 
any  volume  moving  with  the  fluid  is 
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At  any  point,  the  force  acting  on  a  surface  element,  da,  is  (-np  da),  so 
that  also 


Transforming  this  to  a  volume  integral  and  comparing  with  the  previous 
result,  we  find 


Again,  the  integrand  must  vanish,  and  the  momentum  equation  is 


^  +  V-(pu*u^  =  -  Vp 


(5) 


Conservation  of  Energy 

The  energy  per  unit  volume  is  pE  and  the  rate  of  doing  work  on  an 
element  of  its  area  is  -n*pu*  da.  Thus,  with  no  heat  conduction  in  the 
fluid 


*The  quantity  V*  ( pu*  u)  is  obtained  by  considering  individually  the  scalar 
components  of  either  u*  and  combining  results.  Thus,  V«(pu*u)  =  ( pu*V)u> 
+  ut V* pa)  or  V* ( pu* u)  =(u>»V)pu>+  put^7*u) . 
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The  same  sort  of  manipulation  as  "before  leads  to 


+  V*(pEu^  =  - 


Lagrangian  Form  of  the  Equations 

The  three  basic  equations,  (2),  (5),  and  (6),  can  be  written  as 
the  following  equivalent  set: 


Qj:  +  =  -  pV*u* 


St  + 


u-VE  =  - 


+  =  -V.(pu) 


In  all  three  equations  there  appears  the  operator  ^  +  u*«V. 

This  is  called  the  time  differential  operator  along  the  particle  path. 
To  see  the  reason  for  this,  suppose  that  one  wishes  to  find  the  rate  of 
change  of  a  quantity,  A,  per  unit  time  along  some  space-time  path  char¬ 
acterized  by  dr  and  dt.  Then 

^  =  +  *£.VA 

dt  St  dt 


(]  r 

Along  that  special  path  for  which  ^ 


u,  then 


dA  dA 

dt  "*St 


This  special  path  is  just  that  taken  by  a  point  moving  with  the  fluid. 
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o 


where  the  subscripts  are  used  to  show  explicitly  what  is  being  held  con¬ 
stant  in  each  derivative.  The  notation  involving  the  capital  D  will 
occasionally  he  useful  in  this  manual  in  designating  an  infinitesimal 
change  along  the  motion  of  the  fluid. 

Equations  (7),  (8),  and  (9)  in  Eulerian  form  may  he  transformed  to 
the  Lagrangian  form.  The  process  is  most  conveniently  illustrated  in 
one  dimension.  As  an  example,  consider  the  mass  equation  in  which  sub¬ 
scripts  show  explicitly  what  is  being  held  constant: 


We  have  seen  from  ( 1 0)  that  this  can  be  written 


To  complete  the  transformation,  we  note  that 

p  dx  s  pdx 
o  o 

where  pQ  is  the  density  at  the  time  at  which  x  =  xq.  This  expresses 
the  fact  that  the  mass  between  any  two  adjacent  points  is  the  same, 
whether  observed  from  the  Eulerian  or  Lagrangian  viewpoint.  Thus,  with 
p  and  u  being  functions  of  Xq  and  t,  and  pQ  being  a  function  of  xq, 

(12)  becomes 
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This  is  the  Lagrangian  form  of  the  mass  equation. 

The  transformation  can  also  be  performed  in  a  somewhat  more  ele¬ 
gant  fashion.  First,  consider  the  return  from  (l4)  to  (11).  Let  the 
Eulerian  and  Lagrangian  coordinate  systems  coincide  at  time  t  =  0.  We 
then  make  the  coordinate  transformation  in  ( 1 4) : 


(15) 


where  the  integration  is  performed  with  fixed  xq.  Thus 


(16) 


where  use  has  been  made  of  (14).  Then 


This  result  is  the  same  as  (13).  Thus  (16)  becomes 


Also, 
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or 


fb\  fb\ 

fei  =u(^]  + 


08) 


If  we  put  (17)  and  (18)  into  (l4)  and  replace  z  -» t,  then  we  obtain  (ll). 

The  reverse  transformation  by  this  method  is  not  as  simple,  but 
is  worth  presenting  since  it  shows  more  about  the  choice  of  P0(x0)  an<^ 
points  out  the  procedure  required  in  the  more  general  three-dimensional 
equations.  We  start  from  the  Eulerian  equation  with  z  still  replacing  t: 

(&♦  ’&>*  09) 


Define 


0(  Vt}  =  (Jr, 


(20) 


Then,  from  ( 1 5) 


=  1  + 


Jo 

t 


cbA 


dt 


where  we  had  to  use  (16),  because  07)  is  as  yet  not  known  for  this  re¬ 
verse  transformation.  Thus,  using  ( 19),  we  obtain 


LAGRANGIAN  FORM  OF  THE  EQUATIONS 


which  becomes,  with  (l8). 


o 


The  integration  is  to  be  performed  with  xq  constant,  so  that 


This  may  be  integrated  after  division  by  0,  leading  to 


©p  =  K(xq) 


(22) 


where  K( xq)  is  an  "arbitrary"  function  of  its  argument',  which  may  be  deter¬ 
mined,  however,  by  substitution  back  into  the  integral  equation  (21 ): 


0  = 


*  K(xo)  /V 

— r  vst 


1  -  J  - S-  l  &  I  dt 


=  i  + 


^rm  * 


=  i  + 


K(xo)  K(xo) 


To  agree  with  (22),  we  see  that  K(xq)  =  pQ  so  that  combining  with  (20) 
we  finally  obtain 

dx)  fo 
' p 

and  the  transformation  can  proceed  in  a  straightforward  manner. 
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Usually  the  analogous  three-dimensional  transformation  to 
Lagrangian  coordinates  will  involve  such  a  complicated  Jacobian  that  no 
advantage  is  gained.  In  one  dimension,  however,  we  shall  have  several 
occasions  for  using  the  Lagrangian  forms,  and  it  is  useful  to  summarize 
them  here. 


0 


0 


> 


j 


(23) 


Time  Derivatives  of  a  Volume  Integral 

Let  A(?,t)  be  some  quantity  density;  then  there  are  two  main  types 
of  volume  integrals  of  A  which  will  be  of  interest: 


and 


A( Eulerian)  = 


A( Lagrangian) a 


dr 


o 


’o 

The  first  one,  which  we  have  discussed  before,  is  an  integral  over  a 
volume  V  moving  with  the  fluid.  The  second,  in  which  A  is  re-expressed 
as  a  function  of  the  Lagrangian  coordinates,  is  an  integral  of  A  over 
some  fixed  initial  volume. 

We  have  seen  before  (3)  that 
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dA(Eulerian) 

dt 


The  time  derivative  of  the  Lagrangian  integral  is  simpler  to  per¬ 
form;  since  each  element  of  volume  in  the  sum  is  constant,  the  deriva¬ 
tive  is 


o 


Now,  the  time  derivative  of  the  Eulerian  integral  can  he  derived 
in  another  illustrative  way.  We  first  transform  it  to  a  Lagrangian 
integral: 

vc 

where  pQ/ p  is  the  transformation  Jacobian  (whose  form  follows  from  the 
fact  that  p  di  s  pQ  dTQ)#  Thus 


where  all  quantities  in  the  integrand  depend  spatially  on  the  Lagrangian 
coordinates.  Now  the  mass  equation  (l4)  can  be  put  into  the  mixed  form 


r 

o 
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where  the  divergence  is  Eulerian.  Thus 


dA(Eulerian) 

dt 


dT 

o 


si 


$  ♦  *  ** 


dT 


V 

where  use  has  been  made  of  (10).  The  derivative  formula  given  before 
follows  immediately.  This  type  of  transformation  appears  at  first  cum 
bersome,  but  is  sometimes  extremely  useful.  Note  the  special  case 


3E dT  "/p  §  dt 

V  V 

which  follows  from  a  similar  derivation. 

One-Dimensional  Conservative  Form 

In  one  dimension,  the  Eulerian  equations  may  be  written 

^  +  s  (pu2  +  p)  -  ° 

( puE  +  pu)  -  0 

and  the  Lagrangian  equations  may  be  written 
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5<pc/p)  3u  . 

— St-  -  SF  ■  0 

o 


dp  u 
Ko 


dt 


> 


(25) 


Equations  (2*0  and  (25)  are  called  the  conservative  forms.  Integration 
of  any  one  of  them  over  a  fixed  space  interval  (Eulerian  or  Lagrangian, 
as  appropriate)  reveals  the  reason.  Consider  the  Eulerian  momentum 
equation  as  an  example.  With  x^  and  x^  being  fixed  Eulerian  positions, 
we  obtain 

3E  fm  =  (pu2  +  p)x»x,  ■  (p“2  +  p)x=x2 


Thus,  there  is  no  internal  contribution  to  the  timewise  variation  of 
momentum  in  the  interval;  its  momentum  changes  only  if  there  are  boun¬ 
dary  fluxes.  In  this  case,  the  boundary  flux  is  composed  of  two  terms. 
The  transport  term,  pu  ,  measures  the  rate  at  which  momentum  is  carried 
by  the  moving  fluid;  the  force  term,  p,  measures  the  acceleration  due 
to  external  pressures. 

The  conservative  equations  are  all  of  the  form 


dA  ^  SB  _ 

<5t  <5x 


0 
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If  F(x,t)  is  any  arbitrary  function  of  its  arguments,  then 
dF 


B  = 


5t 


A  =  _ 


dF 

5E 


is  the  most  general  solution  of  the  equation.  Now 
dF  =  ^  dt  +  ^  dx  a  E  dt  -  A  a* 


Since  dF  is  a  perfect  differential,  its  integral  around  any  arbitrary 
closed  path  in  the  x-t  plane  is  zero.  Thus 


j j r  (b  dt  -A  dx)  a  0 


This  equation  is  equivalent  to  the  original  one.  Thus,-  the  one¬ 
dimensional  Eulerian  equations  in  integral  form  are 


/(pudt  -  p  dx)  =  0 
/[(pu2  +  p)dt  -  pu  dx 

jT^(puE  +  pu)dt  -  pE  dx 


=  0 


(26) 


and  the  Lagrangian  equations  are 

P, 


^ (u  dt  +  dxQ)  =  0 

/(p  dt  -  pQ  u  dxQ)  =  0 
jf  (pu  dt  -  pQ  E  dxQ)  =  0 


(27) 
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These  equations  in  integral  form  will  be  useful  in  later  developments. 
Other  Coordinate  Systems 

The  equations  for  a  particular  problem  can  often  be  simplified  in 
form  if  a  coordinate  transformation  is  made.  Various  compilations  have 
been  given  of  the  equations  in  the  more  common  coordinate  systems  (see, 
for  instance,  Pai  ). 

A  convenient  starting  point  for  transforming  coordinates  is  the 
set  of  equations  in  general  vector,  Eulerian  form.  Consider,  for  example, 
the  momentum  equation: 

P  +  p(u*V)u  =  -  Vp 

If  tj  is  one  of  the  three  unit  vectors  of  some  curvilinear  coordinate  sys¬ 
tem,  then  the  dot  product  of  q  with  the  equation  will  give  the  appropriate 
component  equation  in  the  desired  system: 

P  5^  +  PR#[(u**V)£|  =  - 

R 

where  u  is  the  component  of  velocity  in  the  direction  of  n,  and  dx  is 

R  ”  T| 

the  change  in  distance  along  a  path  in  the  direction  of  ij.  Now 
R*  t(u-V)£|  =  (u«V)(iHr)  -  #  [(^.V)^] 

The  second  term  on  the  right  vanishes  only  if  ^  is  a  constant  everywhere. 
While  tj  is  always  constant  in  magnitude,  it  is  not  generally  constant  in 

S.  Pai,  "Viscous  Flow  Theory,  I  -  Laminar  Flow,"  D.  Van  Nostrand  Company, 
Inc.,  Princeton,  195 6* 
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direction,  and  the  second  term  does  not  vanish. 

As  an  example,  consider  the  transformation  to  cylindrical  coor¬ 
dinates.  The  three  unit  vectors  in  the  radial,  angular,  and  axial 
directions  are,  respectively,  r,  d,  z.  We  are  thus  concerned  with  find¬ 
ing  V£,  and  Vz.  In  general  in  cylindrical  coordinates, 


*  /n  drf  o  1  c)ti  ~  dn 

^  =  erJe+z& 

where  the  order  of  vectors  is  preserved  throughout.  Thus  we  must  find, 
for  example,  dr/d@.  A  graphical  approach  will  be  illustrative.  The 

vector  r^  -  points  approximately 
in  the  direction  of  6  (and  will  do 
so  exactly  as  6 2  -  6 ^  -+0)»  Also, 

?2  -  r1  has  limiting  magnitude 
( ©2  -  )  times  the  magnitude  of 

r2  or  r1 .  Thus  (?2  -  )/( @2  -  ©1 )  ~  © 

and  dr/de  =  6.  In  similar  fashion, 
the  other  appropriate  derivatives 
may  be  found,  and  =  (l/r)  66, 

V0  =  -  (l/r)  0r,  Vz  -  0.  Thus, 


A  A 


U 


u.[(u*V)r]  = 


u**[(irV)0]  =  - 


u  u„ 
r  6 


and  the  component  equations  of 
motion  are: 
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du 

r  /-* „\  w 

P^-  +  p(u*V)ur  -  P  —  =  - 


p  5T  +  p(u  v)ue  p  T~ 

8u  ^ 

P  +  P(U*V)U  =  «  ;£ 


Notice,  in  passing,  that  one  aspect  of  these  equations  can  he 
given  what  at  first  seems  a  puzzling  interpretation.  If  u^  =  0,  then 
(28)  becomes 


Thus 


[2  Du 
P  tvT“  <*x 


where  dx  is  the  element  of  volume  of  a  cylindrical  shell  of  length  L. 


Thus 


r2,o  chi 


3o  5T  aio 


dx  =  -  2rtL  /  r  dr 


Define  the  mean  radial  speed, u  ,  of  the  shell  contained  between  r1  and 


r2*y 


p  u  dT  =  m  u 
o  r  o  r 


where  m  is  the  shell  mass.  Then 


-35- 


I.  THE  BASIC  EQUATIONS 


dr 


=  (PA), 


inside 


-  (PA) 


outside 


+  2*L 


r 

l 


pdr 


where  A  is  the  surface  area  of  the  inside  or  outside  of  the  shell.  Thus 
the  acceleration  of  the  shell  is  produced  by  more  than  the  difference 
between  the  external  forces;  even  if  these  vanish,  the  internal  pressure 
within  the  cell  causes  its  radial  acceleration. 


Equation  of  State 

The  equations  so  far  discussed  apply  quite  generally  to  all  simple 
fluids.  Reference  to  any  set  of  them  —  (24),  for  example  —  shows  that 
the  number  of  equations  is  always  one  less  than  the  number  of  dependent 
variables,  and  thus  insufficient  for  a  complete  description  of  any  pro¬ 
cess.  The  required  additional  relation  is  one  characterizing  the  fluid 
itself.  It  is  called  the  equation  of  state;  we  often  write  it  in  the 
form 

P  =  f(p,T)  .(31) 

and  accept  its  existence  as  being  an  experimental  fact. 

It  has  been  observed  that  many  gases  closely  follow  the  "ideal 
gas"  equation  of  state: 
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THE  EQUATIONS  FOR  A  FLUID  WITH  NONLOCAL  FORCES 

where  R  is  the  universal  gas  constant  and  N  is  the  number  of  moles  per 
unit  mass  (or  the  reciprocal  of  the  molecular  weight).  Much  of  our  dis¬ 
cussion  will  relate  to  ideal  gases. 

The  Equations  for  a  Fluid  with  Nonlocal  Forces 

If  each  element  of  fluid  is  subject  to  forces  exerted  by  other 
than  its  immediate  neighboring  elements,  then  additional  terms  are  re¬ 
quired  in  the  equations.  (Such  forces  would  be  present,  for  example, 
in  a  gas  with  net  electric  charge,  or  one  acted  upon  by  an  external 
gravitational  field. )  The  mass  equation  remains  unchanged  by  such 
forces.  If  the  force  per  unit  mass  (the  acceleration)  is  g,  then  the 
momentum  equation  becomes,  in  Eulerian  coordinates, 

|+(u-V)?=-ivp  +g  (33) 

The  rate  at  which  work  is  done  by  the  external  force  on  the  element  is 
VL.g.  Thus  the  energy  equation  becomes 

||  +  (u*V)E  =  -  I  V.(p3)  +  u*g  (34) 

Suppose  now  that 

2-31<3+3,<r,t) 

where  gj  is  the  acceleration  produced  by  forces  from  outside  the  fluid, 
and  gg  is  that  produced  by  distant  action  on  elements  within  the  fluid. 
Suppose  further  that  and  gg  both  can  be  expressed  as  gradients  of 
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functions 

g ]  =  -  v  ^(r) 
^  =  -  V  cp2(r,t) 
and  finally,  that 


where  the  integration  is  over  the  entire  volume  of  the  fluid,  and 

Q(r  -  ?' )  s  -  j*).  These  assumptions  are  not  very  restrictive;  they 

are  valid  for  many  physical  situations  of  interest. 

The  momentum  equation  (33)  can  then  he  written 


^jr  +  V*(puu*+  p)  =  p  g|  -  /p(r^t)p(r,,,t)  $( r  -  r*)  dr' 

J  3? 

so  that,  to  demonstrate  momentum  conservation,  we  integrate  over  volume 


where  vanishing  of  the  integral  of  the  second  term  on  the  right  follows 
from 


Second  term  =  JJ  p(r,t)p(^*,t) 


fp(r,t)p(r',t) 
=  -  (Second  term) 


t* )  drdT» 

-  r)  dTdT1 


Thus,  the  total  momentum  of  the  fluid  is  affected  only  by  the  external 
forces. 
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The  energy  equation  ( )  can  be  treated  in  similar  fashion.  We  re¬ 
write  it  in  identically  equivalent  form 

^  pE  +  5p(91  +  92)  +  V.  pufe  +  i  put<?1  +  92)  +  pu  = 

rd(cp,  +  9p)  _» 

2  p  L — '  u'v  (cpi  +  ^ 

wherein  use  has  been  made  of  the  mass  conservation  equation  (2).  Again  we 
integrate  this  over  the  whole  fluid  and  arrive  at  the  result 

^  f  pE  +  ip(cp1  +  <P2)  dT  =  0 


Vanishing  of  the  integral  of  the  right  side  of  the  equation  is  demon¬ 


strated  as  follows: 


Integral 


r  +  %)  . 

■  *>  Har-2-  -  *** 


1  +  92)J  <*t 


o  r  ^cp  +  q>  )  ; 

i  /  P  --a— -  +  (9,  +  <P2)V-pu  dr 


where  we  have  dropped  a  term  which  integrates  to  zero.  Using  the  mass 


equation  we  get 


Integral 


!  r  r  3(<pi +  <p2) 


-*/[■ 


"  ((pi  +  *2*  5tl  dT 


=  i  *1  h  fpdx  +  ^/(p  ~d  ~  92 


-37- 


I.  THE  BASIC  EQUATIONS 


The  first  term  vanishes  since  total  mass  is  conserved.  Consider  one 
term  in  the  second  part: 


dT 


p(3?',t)  didr' 


p(?,t)  dxdr* 


so  that  the  second  integral  is  self-cancelling.  Thus  it  is  shown  that 
the  total  "potential"  energy  of  the  system  is 

i  Jp +  <Pg)  dx 

and  this  plus  the  kinetic  and  internal  energy  of  the  system  are  altogether 
conserved. 
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CHAPTER  II 

DISCUSSION  OP  THE  EQUATIONS 
The  Polytropic  Equation  of  State 

The  specific  entropy  change,  DS,  along  the  motion  of  a  point  is  de- 

* 

fined  by  the  equation 

TDS  =  DI  -  ^  Dp  (1) 

P 

The  second  law  of  thermodynamics  states  that  DS  is  a  perfect  differen¬ 
tial;  that  is,  that  the  entropy  at  any  point  in  a  fluid  is  a  function  of 
the  instantaneous  state  at  that  point.  For  an  ideal  gas,  with  equation 
of  state  (l-32), 

DS  =  1  DI  -  M  Dp 

Now,  since  DS  is  a  perfect  differential,  here  considered  a  function  of 
the  independent  variables  p  and  I,  we  have 

DS  -|f  DI  +  ||  Dp 

See,  for  example,  P.  S.  Epstein,  "Textbook  of  Thermodynamics,"  John  Wiley 
and  Sons,  Inc.,  New  York  1937* 
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II.  DISCUSSION  OF  THE  EQUATIONS 


Comparing  these  two  equations,  we  see  that  the  condition  of  exactness 
applied  to  an  ideal  gas  is 

dS  =  2  US  NR 

dl  T'  ^  ~  “  p 

I’rom  which 


so  that  T  is  a  function  of  I  alone.  In  an  ideal  gas,  the  specific  in¬ 
ternal  energy  is  a  function  of  temperature  only.  For  many  gases,  it  is 
sufficient  to  keep  only  a  few  terms  in  the  Maclaurin  expansion  of  I  as 
a  function  of  T: 

2 

I  =  a  +  bT  +  cT  +••• 

The  internal  energy  can  he  chosen  to  vanish  at  T  =  0.  (T  is  thus  the 
absolute  temperature. )  Those  special  gases  for  which  only  b  /  0  are 
called  "polytropic."  Mathematically,  they  are  relatively  simple  to 
handle  and,  under  standard  conditions,  they  often  represent  real  gases 
quite  closely. 

For  such  a  polytropic  gas  it  is  customary  to  define  the  dimension¬ 
less  constant: 


in  terms  of  which  the  polytropic  equation  of  state  can  be  written 

P  =  (y  -  l)pl  ( 2 ) 
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CHANGES  OF  STATE  IN  A  SIMPLE  GAS 


The  significance  of  the  constant,  b,  can  be  seen  by  a  considera¬ 
tion  of  the  process  of  adiabatic  (reversible)  heating  of  the  gas.  In 
such  a  process  the  heat  change  along  the  motion  of  a  point  is 
DQ  =  TDS,  so  that,  from  ( 1 ) 


DQ  =  bDT  -  ^  Dp  s  bDT  +  pDV 


Thus,  if  the  volume  is  kept  constant,  then 


DQ 


=  b 


V=constant 


so  that  b  is  the  specific  heat  at  constant  volume.  If,  on  the  other  hand, 
the  temperature  is  changed  at  constant  pressure,  then,  since  DV  =  (NR/p)DT 
[obtained  from  the  ideal-gas  equation  of  state,  (l-32)],  we  obtain 


=  b  +  NR 


^  xp=constant 

which  is  the  specific  heat  at  constant  pressure, 
these  two  specific  heats. 


Thus  y  is  the  ratio  of 


Changes  of  State  in  a  Simple  Gas 

As  long  as  the  field  variables  are  continuous,  the  differential 
equations  of  Chapter  I  are  meaningful;  we  examine  some  of  their  proper¬ 
ties  starting  from  the  Eulerian  form 


pV*u 
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XI.  DISCUSSION  OF  THE  EQUATIONS 


Du  ^ 

PDt='VP 

DE  „  -» 

P  Dt  =  “  V‘pu 


For  E,  put  I  +  \  U'U,  and  the  third  equation  becomes 


i  ^ 

DI  — >  Du  — >  — > 

p  Dt  +  pU’Dt  =  -  pV*U  “  U*^ 


Eliminate  Du/Dt  by  means  of  the  momentum  equation  and  V»u  by  means  of 
the  mass  equation,  and  the  result  is 

DI  -  ^  Dp  =  0  (3) 

P 

Comparison  with  ( 1 )  shows  that  the  specific  entropy  of  each  point  re¬ 
mains  constant  as  it  moves  with  the  fluid.  The  value  may  differ  from 
point  to  point;  if  it  does  not,  then  the  fluid  is  called  "isentropic, " 
and  the  above  result  shows  that  it  remains  isentropic. 

One  of  the  field  variables  can  be  removed  from  (3)  by  use  of  the 
equation  of  state,  and  the  result  is  an  equation  which  always  can  be  in¬ 
tegrated,  yielding  the  "adiabatic  equation  of  state."  If,  for  example, 

I  is  eliminated  by  means  of  the  polytropic  equation  (2),  the  integrated 
result  is 


p  =  A(rQ)p7  (4) 

where  AC?^)  is  an  arbitrary  function  of  the  Lagrangian  coordinate.  This 
is  the  adiabatic  equation  of  state  for  a  polytropic  gas. 

The  quantity  aC?^)  is  related  to  the  entropy;  the  relation  can  be 
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CHANGES  OF  STATE  IN  A  SIMPLE  GAS 


seen  as  follows.  Suppose  that  by  some  process  the  values  of  p  and  p  are 
changed  at  some  point  in  such  a  way  that  the  entropy  is  changed.  Eli¬ 
minate  T  and  I  as  they  occur  in  (l),  (2),  and (1-32)  to  obtain 


Dp  _  2®P  =  d 
P  P 


From  (4), 

DS  =  -^-1  D  fin  A(r  ) 
y  -  1  O 

Thus, 

S  =  111  A(?)  +  constant  a  Sq  +  ^  —  In  (-^)  (5) 

Note  that  as  T  -»  0  for  constant  density,  p  -» 0  and  S  -  ®.  Thus,  no  gas 
can  behave  as  polytropic  for  very  low  temperatures. 

For  any  simple  fluid,  the  entropy  of  any  element  remains  constant 
as  long  as  the  equations  of  motion  are  valid.  It  will  be  shown  later, 
however,  that  certain  initial -value,  boundary -value  problems  do  not 
possess  solutions  which  are  everywhere  continuous,  even  if  the  initial 
and  boundary  conditions  are  continuous.  Moving  surfaces  of  discontin¬ 
uity  can  arise,  and  at  them  the  equations  become  meaningless.  A  special 
treatment  is  required,  and  one  can  be  found  which  gives  a  useful  meaning 
to  these  discontinuities,  but  a  consequence  of  the  treatment  is  that 
there  is  a  change  in  entropy  as  a  discontinuity  passes  over  an  element 
of  the  fluid. 
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II.  DISCUSSION  OF  THE  EQUATIONS 


As  soon  as  the  real  effects  of  viscosity  and  heat  conduction  are 
added  to  the  equations,  one  no  longer  can  derive  the  result  DS  =  0.  Heat 
sources  due  to  chemical  or  nuclear  reactions  also  result  in  entropy  pro¬ 
duction.  These  properties  of  nonsimple  gases  will  he  exhibited  in  Chap¬ 
ters  VIII  and  IX. 


Angular  Velocity  and  Circulation 

Relative  to  some  axis  defined  by  the  unit  vector  q,  the  angular 
velocity  of  a  fluid  element  is 


o> 

<1 


_  q«r  X  u 


where  r  is  a  vector  from  an  origin  on  the  axis  at  the  nearest  point  to 
the  element. 

To  find  the  average  angular  velocity  about  q  in  the  neighborhood  of 
a  point  on  q,  we  proceed  as  follows.  First,  find  the  average  on  a  circle 
centered  at  the  point  and  perpendicular  to  q.  This  average  is  given  by 


*  %  ■  ^  |  =  I/S  X  V?  da 


where  the  line  integral  has  been  transformed  to  a  surface  integral  over 
the  plane  surface  bounded  by  the  circle  whose  radius  is  R.  Thus,  the 
value  is 


u> 


=  — —x  q*  /  V  X  u  da 
2nEC  J 


If  the  circle  is  small  enough  that  V  x  u  varies  little  on  it,  then 
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ANGUIAR  VELOCITY  AND  CIRCULATION 


J  V  x  u  da  ->  «R2V  x  u;  in  the  limit  as  R  -*  0  the  replacement  becomes 
correct  if  the  velocity  curl  is  continuous.  Thus,  the  average  angular 
velocity  component  over  the  circle  is 

cd  =  •§■  q*V  X  u 

q. 


independent  of  the  circle  size,  and  it  follows  that  the  angular  veloc¬ 
ity  of  a  point  in  the  fluid  is 


a?  =  k  VX  u  (6) 

q 


The  quantity  VX  u  is  called  the  vorticity. 

Actually,  the  averaging  does  'not  need  to  be  performed  over  a  cir 
cular  path.  Thus,  for  an  angular  averaging  over  any  path  in  a  plane 


perpendicular  to  q, 


CD  =  - — 

q  2« , 


2* 


q  X  r*u 


dd 


where  6  is  the  angle  between  r  and  some  fixed  initial  value  of  r.  Now 
?  =  rn  where  n  is  a  unit  vector  which  is  a  function  of  6  only;  likewise, 
r  can  be  considered  a  function  of  0. 

Expand  u  in  a  series  about  the  origin: 

u  =  u  +  (r*V)u  +  • • • 
o  o 

and  restrict  r( e)  to  be  so  small  that  the  terms  beyond  those  shown  can 
be  neglected.  Then 
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CD 


.  —  [** 
2k  J 


q  X  n* 


U_  /v  . 

—  +  (n-V)u 
r  '  o 


ae 


Integration  of  the  first  term  in  the  bracket  gives  zero,  so  that 

2* 

“q  ■  k  fo  4  *  a-(8^)“0M 

To  perform  this  integration,  we  put  R  s  nR  wherein  R  =  constant.  Then 
Rq  X  n  d@  s  dR*  and 

“4*53-  / 


(R-V)u  *d R 
'  o 


R 


Now,  to  the  same  order  of  approximation  as  before,  u*  s  u^  +  (R*V)u^. 


Thus,  . 

-  _  1  f  u-dR 

®q  2*R  J  R 

which  is  the  same  form  as  before  with  integration  around  a  circle.  Note 
that  this  result  would  not  have  been  obtained  for  lengthwise  averaging 
along  an  arbitrary  path. 

To  see  how  the  angular  velocity  of  the  field  changes  with  time 
in  a  simple  fluid,  we  start  with  the  Eulerian  equation 

du*  1 


+  (u*V)u  =  -  —  Vp 
P 


For  a  simple  isentropic  fluid,  p  is  a  function  of  p  only,  and  the  right 
side  of  the  equation  can  be  expressed  as  the  gradient  of  a  function. 
Taking  the  curl  of  the  equation,  we  obtain 

d® 


2  H  +  V  X[(u‘V)u}  =  0 
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Now 

(u.V)u  =  i  V(u.u)  -?X(VX5) 

and  so 

|=7X  (uxS) 

After  expansion  of  the  right  side,  this  can  be  transformed  to 
=  (io»V)u  -  d^V»u) 

Note  that  if  a?  =  0  at  some  instant  and  at  some  element  of  the  fluid, 
then  its  rate  of  change  along  the  motion  of  the  element  is  also  zero. 
Thus,  if  ms  0  for  the  entire  fluid  at  some  instant,  then  it  continues 
to  vanish  for  all  later  times,  regardless  of  the  motion  of  boundaries 
(provided,  of  course,  that  the  fluid  remains  isentropic). 

If,  at  some  time,  the  vorticity  is  nonzero,  then  it  will,  in 
general,  change  with  time.  There  is,  however,  for  a  simple  isentropic 
fluid  a  related  invariant  called  the  "circulation,"  F.  It  is  defined 

r  =y  u»d?  (7) 

where  the  integration  is  performed  about  some  closed  loop,  every  point 
of  which  travels  with  the  fluid.  The  relation  of  circulation  to  vor¬ 
ticity  may  be  demonstrated  by  transforming  the  integral  to  a  surface 
integral  over  any  continuous  surface  bounded  by  the  loop.  Let  n  be  a 
unit  normal  to  the  surface  and  da  the  area  element.  Then 
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=  J (n  X  V)  'u  da 


2/  n»co  da 


To  prove  the  invariance  of  T  under  the  conditions  stated,  we  differen¬ 
tiate  along  the  fluid  motion 


Dr  iYbu  ,!->  ,  Ddr 
Dt=iU  ,dr  +  U,Dt 


The  quantity  can  be  interpreted  by  reference  to  the  diagram: 


u  and  u'  are  the  velocities  at  each  end  of  dr  at  time  t.  They  carry  dr 
to  a  new  position  at  time  t  +  St,  such  that 


Thus, 


ufct  +  (dr)t+5t  =  u'St  +  (dr)t 


=  lim  (£?*  -  u) 

5t-»0 


Thus, 


!?'={?+  (d?*V)u* 


=  (d?*V)u* 
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and 

ff.SjjE  .  u-(dr-V)u  s  4  (d?-v)(?.?) 

Now  also,  for  this  fluid,  is  the  gradient  of  some  function,  say  4 
(that  is,  =  ~  VP)  ^  so 

j|-/v 

This  can  he  transformed  to  a  surface  integral  which  will  involve  the  curl 
of  the  integrand  as  a  factor,  and  thus  vanish,  proving  the  invariance  of  f. 

Note  that  it  follows  from  (8)  that  in  two-dimensional  plane  flow, 
o?V  is  a  fixed  constant  for  each  point  moving  with  the  fluid.  This  fact 
can  also  be  derived  as  follows.  In  three-dimensional  flow  of  a  simple, 
isentropic  fluid,  we  have  seen  that 

=  (a?- V)u*  -  u^V«u) 

Thus 

OT“  =  i?5J  +  vgf  =  J?vv.i?+  V[(i?-V)?-  Stv-3)]  =  v(5?-v)? 

Dt  Dt  ut 

In  the  special  two-dimensional  case,  is  normal  to  the  surface  across 
which  u  varies  and  the  right  side  vanishes. 

Angular  Momentum 

The  angular  momentum  about  the  origin  of  coordinates  of  an  element 
of  fluid  of  volume  dr  is 
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dm  r a  — >  ., 

—  =  j  nx  pr  da 

where  n  is  the  unit  outward  normal  at  the  surface  area  element  da. 

This  result  shows  that  angular  momentum  of  the  fluid  changes  with  time 
only  through  a  surface  flux;  the  angular  momentum  of  an  isolated  fluid 
is  identically  constant. 
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Local  Decomposition  of  the  Velocity  Field 

At  any  instant  of  time,  and  for  any  point  in  the  fluid  (chosen  to 
he  the  origin  of  coordinates),  the  velocity  field  can  be  expanded  in  a 
Maclaurin  series: 

at?)  a  ?  +  (?-V)£  +  a  (r2) 

where  (r2)  means  neglected  terms  of  second  or  higher  order. 

We  now  introduce  component  notation  in  which  i,  j ,  and  k  can  indi 
vidually  designate  x,  y,  or  z.  Thus  u^  can  refer  to  ux  or  uy  or 
In  addition  we  shall  employ  the  summation  convention  such  that  if  any 
index  appears  twice  in  a  term,  then  summation  over  all  three  coordi¬ 
nates  is  implied.  Thus,  for  example, 

du.  du.  duj  ^ 

uo  ST  s  ux  TST  +  uy  W  +  u*  ^ 

u 


in  which  expression  i  may  still  be  x  or  y  or  z.  With  this  notation 
the  Maclaurin  expansion  becomes 

^io  frt  2x 

"i  3  ui°  +  xj  -35EJ  i<r(r  ) 


We  now  define  the  symmetric  and  antisymmetric  tensors 
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u.  =  u.  +  x.(d.,  +  Cl..  )  +  0( r  ) 
1  10  A  lj  xjyo 


The  significance  of  D. .  can  be  seen  as  follows.  The  length  of  a  dis- 

placement,  ds,  corresponding  to  coordinate  displacements  of  dx,  dy, 

2 

and  dz  is  given  by  ds  =  dx^dx^.  Now 


(dxi}  =  §t  (^o  tXjo)  =  ^  taJo  =  toJ 


where  subscript  zero  here  means  the  Lagrangian  coordinate,  which  does 
not  change  with  time.  Thus 


D  ( ,  2 

5t  Vds 


Thus 

du. 

a2dxi 

du 

s  2  dx. 
J 

4 

where  the  two  forms  arise  from  interchanges  of  the  dummy  summation  in¬ 
dexes.  If  we  take  the  average  of  the  two  forms,  we  get  the  result 


m  vdsy  - 2  Dij 


Thus,  if  D. .  =  0  at  some  point,  then  the  length  of  any  displacement 
x  J 

near  that  point  will  remain  constant  as  the  point  moves  with  the 

fluid.  Likewise  if  D. ,  a  o  everywhere,  then  the  fluid  is  perfectly 

rigid  and  its  motions  are  limited  to  translation  and  rotation  only. 

Local  vanishing  of  D^  also  means  local  translation  and  rotation  only, 

so  that  x.Cl. .  is  thereby  interpreted  as  representing  the  rotational 
J  J 

component,  while  D„  represents  the  deformation. 
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Note  that 


XA) =  4  A  !x*V  +  * <y2) 

so  that,  with 


we  have 


ui  =  uio  + 


3d 

3*. 


+  'Yil’c,  *<’(r‘) 


A  straightforward  comparison  of  components  of  and  the  ith  compon¬ 

ent  of  i(a?x  r)  shows  that  they  differ  bytftr2)  so  that,  finally, 

u  =  u  +  (VD)  +  X  r)  +  <7(r2) 


The  velocity  field  has  been  decomposed  locally  into  a  translation,  a 
deformation,  and  a  rotation. 


Propagation  of  Small  Disturbances 

What  is  the  manner  in  which  a  very  small  disturbance  is  propagated 
through  a  region  in  which  the  field  variables  are  otherwise  constant? 
The  region  is  isentropic  so  that  everywhere  the  pressure  is  a  function 
of  density  only, 

P  =  p(p) 


We  use  the  Eulerian  equations 
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du 
p  St 


+  p(u*V)u  =  -  Vp 


+  pW+  (u*V)p  =  0 

Define  a  positive  quantity  c  (which  we  will  here  show  to  he  the  sound 
speed)  such  that 

c  =  (9) 

dp 

and  also  substitute 


u  -»u  +  u, 

o  I 


P  -*  p0  +  Pt 


c  -*  c  +  c, 
o  1 


where  p  ,  u  ,  and  c  characterize  the  undisturbed  field,  and  p, ,  u. , 

'o  o  o  .'ll' 

and  c  are  relatively  very  small  perturbations.  Then  the  equations 
become,  to  lowest  order  in  the  perturbations, 

sr  *  m  -  %  Vpi 

-3T  +  pV-Z,  +  -v)p,  »  o 

We  now  transfer  to  a  coordinate  system  moving  with  velocity  u^  relative 
to  the  original  fixed  system.  In  this  new  coordinate  system,  the  equa¬ 
tions  become 


SOUND  WAVES  IN  A  VARIABLE -DENSITY  ATMOSPHERE 


1  2 

p0  st  ■  -  co  Vp. 

\  v-» 
s r  -  -  fov-ui 


From  these  two  equations  u^  can  he  eliminated,  giving 


=  0 


This  is  the  familiar  wave  equation,  for  which  the  general  solutions  are 

waves  propagating  with  speed  cq.  Thus  small  disturbances  are  propagated 

relative  to  the  moving  fluid  with  speed  c  . 

-  ■=•  0 

For  a  polytropic  gas,  for  example. 


=y?  =  j  7(7  -  i)i 


p 

2 


(10) 


In  general,  c  is  the  constant-entropy  pressure  derivative  with 
respect  to  density: 

c2  =  &) 

^  P'S  ss  constant 

With  a  general  equation  of  state  in  the  form  I  =  I  (p,p),  the  formula  for 
sound  speed  can  thus  be  written 


2  dl 

2  p  -  p  ^ 
0  "  2  dl 

p 


Sound  Waves  in  a  Variable-Density  Atmosphere 

The  time-independent  solution  of  the  one-dimensional  Eulerian 
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equations  with  external  acceleration 

p  3?  +  pu  5E  =  "  °2  ll  +  pst(x) 

3t  +  p^+u3§=0 


gives  the  equilibrium  state  of  a  variable  density  fluid.  We  assume 
that  the  sound  speed  is  constant.  Then  the  equilibrium  solution  is 
given  as 

u  (x)  -  0 
e  ' 

Pe(x)  -  PQ  exP 

which  includes  the  boundary  condition  pg  -  pQ  at  x  =  0.  We  look  for 
solutions  differing  only  slightly  from  the  equilibrium  state 


4 yj  g(x')dx' 


u  «  c 

p  =  pe(x)e€,  €«  1 

Substitution  of  these  into  the  equations  and  elimination  of  all  but 
lowest  order  terms  lead  to 

de  du  u  dpe 

3t  +  cbc  ”  ~  p  dx 
e 


du  2  de 
+  c  = 


o 


Note  that  if  pg  is  constant,  then  these  become  the  wave  equation  with 
signal  speed  c. 
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Consider  the  special  case  in  which  g(x)  =  constant.  Then 


pe(x)  =  pQ  e' 


gx/c* 


and  the  equations,  with  e  eliminated,  become 


d2u  2  d2u  du  _  n 

72  "  c  72  ‘  g^“  ° 

du  dx 


Substitution  of  the  progressive-wave  solution 


iax  iu>t 
u  =  A  e  e 


shows  that  the  condition  for  solution  is 


2  2  2  .  _ 

05  =  a  c  -  lga 

and  that  A  is  an  arbitrary  constant,  with  dimensions  of  velocity 
(whose  magnitude,  of  course,  must  be  much  less  than  c  for  validity  of 
the  perturbation  approximation) .  If  the  wave  length  of  the  signal  is 
small  compared  with  the  distance  over  which  Pe/pD  changes  appreciably, 
then  (d  can  be  approximated  by 

cd  *»  ±  (ac  -  ||) 

Thus  a)  has  an  imaginary  part,  corresponding  to  a  change  in  amplitude 
of  the  wave;  as  the  wave  moves  in  the  direction  of  decreasing  density , 
its  amplitude  increases,  and  vice  versa.  Along  the  line  x  =  +  ct, 
u  =  A  exp(-  gx/2c2)  so  that,  along  the  motion  of  the  signal,  p  u  is 
constant. 
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CHARACTERISTIC  SOLUTIONS  AND  SIMPLE  WAVES 


Characteristics 

One  of  the  most  powerful  methods  for  obtaining  solutions  of  the  one¬ 
dimensional  hydrodynamic  equations  is  called  the  method  of  characteristics. 
As  an  example  of  its  application,  consider  the  problem  of  solving  the  one¬ 
dimensional  Eulerian  equations  for  a  simple,  isentropic  gas 
du  du  2  dp 

p  St  +  pu  S  =  ■  c  S 

(1) 

dp  dp  du 

dt +  US=  -  P^E 

The  sound  speed,  c,  is  a  function  of  the  density,  so  that  we  may  intro¬ 
duce  a  new  function,  a,  defined,  to  within  an  arbitrary  constant,  by 

db  =  (2) 

P 

Then  ( 1 )  becomes 


du 

du 

dcr 

3t  + 

U^  + 

C^  = 

d a 

do 

du 

3t  + 

us  + 

C  3£  = 
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wherein  c  is  now  considered  to  be  a  known  function  of  a.  By  summing  or 
differencing  these  equations,  one  obtains 


££-(u  +  o)  +  (u  +  c)  (u  +  a)  =  0 

^  (u  -  a)  +  (u  -  c)  (u  -  a)  =  0 


(4) 


From  the  first  equation,  one  can  see  that  along  a  line  in  the  x-t 
plane  such  that  dx/dt  =  u  +  c,  the  quantity  (u  +  ff)  is  a  perfect  dif 
ferential  and  the  equation  can  be  integrated.  With  a  similar  result 
from  the  second  equation  we  may  write 


u  +  o  =  constant,  along  —  =  u  +  c  ^ 

fly 

u  -  o  =  constant,  along  —  =  u  -  c 

at 

These  are  the  characteristic  solutions,  and  the  families  of  lines 
dx/dt  =  u  ±  c  are  called  the  characteristic  lines,  or  simply,  the 
characteristics . 

These  characteristic  solutions  are  not,  in  general,  complete, 
they  do  not  necessarily  allow  the  features  of  any  applicable  flow 
field  to  be  determined  directly.  They  are,  however,  extremely  power¬ 
ful  aids  in  obtaining  solutions  in  certain  special  cases,  or  in  cases 
when  parts  of  the  solution  can  be  obtained  by  other  means. 

As  an  example  of  the  use  of  the  characteristic  solutions,  con¬ 
sider  the  problem  of  determining  the  effect  on  a  gas,  initially 
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at  rest,  of  a  piston  being  withdrawn 

signal  from  it.  Up  to  t  =  0,  there  is  a 

piston  at  x  =  0,  and  gas  at  rest 
o 

for  x  >  0.  At  t  =  0,  the  piston 
commences  to  move  with  uniform 
velocity,  u^,  in  the  negative  x- 
direction,  and  a  sound  signal  pro¬ 
ceeds  into  the  gas.  We  now  attempt 
to  connect  the  sound  signal  and 
piston-path  lines  with  a  family  of 
characteristic  lines.  The  line  dx/dt  =  u  +  c  is  of  no  use,  since  along 
the  sound  signal,  u  =  0,  and  the  characteristic  lies  along  the  signal 
line.  The  line  dx/dt  =  u  -  c,  on  the  other  hand,  has  negative  slope  at 
the  sound  signal  line,  hence  intersects  it.  Furthermore,  this  charac¬ 
teristic  has  negative  slope  throughout  the  flow  field  between  the  sound 
signal  and  the  piston  path,  because  u  <  0  and  c  >  0.  Also,  the  slope 
is  more  negative  than  that  of  the  piston  path,  assuming  that  c  /  0  any¬ 
where  within  the  flow  field  of  interest,  so  that  each  characteristic 
dx/dt  =  u  -  c  will  also  intersect  the  piston  path.  Thus 

U  -  0  =  U  -  0 

o  o  p  p 

or,  since  u  =  0, 

*  o 

a  =  u  +  0  (6) 

p  p  o 

For  a  polytropic  gas,  for  example,  o  e  2c/(y  -  l),  so  that 
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c 

P 


u 

p 


+  c 

o 


Since  u  <0,  the  sound  speed  is  less  than  c  at  the  piston;  tne  more 
p  o 

negative  the  piston  speedy  the  smaller  would  be  the  sound  speed  there. 
This  fact  makes  plausible  the  assumption  that  c  /  0  in  the  flow  field, 
with  the  exception  being  the  case  in  which  c^  =  0.  If  the  piston  is 
withdrawn  any  faster  than  the  critical  speed  at  which  c =  0,  then  the 

x' 

gas  cannot  follow;  a  vacuum  occurs  between  the  escaping  gas  front  and 
the  piston.  This  critical  piston  speed  —  called  the  escape  speed  of 
the  gas  —  is 

2c 

u  - 2-  (7) 

escape  y  -  1 

Notice  that  c  is  independent  of  time,  simply  because  u  and  c  are 
P  P  Q 

independent  of  time.  It  does  not  matter  where  the  two  intersection 

points  are  along  the  sound-signal  and  piston  paths. 

The  method  of  solution  is  also  valid  if  u  varies  with  time,  as 

P 

long  as  a  shock  does  not  interfere  with  the  characteristic  line. 

(Note:  If  the  piston  velocity  persists  for  sufficient  time  at  values 
less  negative  than  any  which  it  has  previously  attained,  then  a  shock 
will  form.) 

An  interesting  generalization  of  this  problem  can  also  be  solved 
by  the  method  of  characteristics.  Suppose  that  the  piston  is  replaced 
by  a  wall  which  has  mass  per  unit  area  m.  Up  to  t  =  0,  the  wall  is 
held  fixed  at  x  =  0;  it  is  then  released  and  moves  away  because  of  the 
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gas  pressure  exerted  on  it.  What  is  the  history  of  its  motion? 

script  w  refer  to  conditions  at  the  movable  wall.  Then 

du 

w 


Let  sub- 


The  characteristic  equation  (6)  is  again  applicable;  in  terms  of  sound 
speed. 


Now 


cw 


c 

o 


Z± 


27 


so  that  we  now  have  three  equations  for  determination  of  the  three  un¬ 
knowns,  u^,  p^,  and  c^.  There  is  one  constant  of  integration  to  be  de¬ 
termined  in  the  solution;  for  this  we  specify  that  u  =  0  at  t  =  0. 

w 

Then  the  solution,  after  some  manipulation,  is  found  to  be 


where 


Thus  as  t  -»  ob,  u^  approaches  the  gas  escape  speed;  the  wall  mass  only 
affects  the  rate  at  which  the  final  speed  is  approached. 
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Simple  Waves 

Suppose  that  in  some  region  of  the  x — t  plane  the  flow  field  is  in 
a  constant  state:  u  =  uq,  o  =  Oq  everywhere  within  the  region.  The 
constant -state  region  can  be  bounded 
above  and/or  below  by  regions  in 
which  the  flow  field  is  not  constant. 

The  boundary  lines  will  either  be 
shocks,  or  else  they  will  be  straight 
characteristics.  If  they  are  not 
shocks,  then  the  boundaries  of  the 
nonconstant  regions  (the  disturbance 
boundaries)  will  propagate  into  the 
constant -state  region  with  sound  speed  relative  to  the  material,  and 
thus  will  have  slope  uq  ±  cq,  being,  therefore,  straight  characteristics. 

The  most  important  fact,  now  to  be  established,  is  that  the  flow 
in  the  adjacent  regions  will  always  be  of  a  particularly  simple  form. 
Consider  first  the  lower  nonconstant  region.  Through  it  will  pass  the 
family  of  characteristics  dx/dt  =  u  -  c,  which,  at  the  disturbance 
boundary,  have  slope  uq  -  cq,  and  hence  intersect  it  as  long  as  cq  /  0. 
Along  each  of  these  characteristics  u  -  o  is  constant;  indeed,  u  -  a 
will  be  the  same  constant  for  all  members  of  the  family  of  characteris¬ 
tics  which  intersect  the  disturbance  boundary,  since  along  that  boundary 

the  constant  is  u  -  a  .  Thus,  throughout  the  region  adjacent  to  the 
o  o 

region  of  constant  state,  u  -  a  will  be  one  fixed  constant.  Likewise, 
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one  can  prove  that  throughout  the  nonconstant  region  above  the  region 
of  constant  state,  u  +  a  will  be  one  fixed  constant.  Any  region  in 
which  u  +  fl  or  u  ■  ff  is  a  fixed  constant  throughout  is  referred  to  as 
the  region  of  a  "simple  wave. "  The  extent  of  the  region  of  a  simple 
wave  adjacent  to  a  region  of  constant  state  is  limited  by  the  require¬ 
ment  of  contact  with  the  disturbance  line  through  a  characteristic. 

In  the  case  that  cq  =  0,  no  infinitesimal  disturbance  can  pro¬ 
pagate  into  the  constant-state  region;  any  wave  propagating  into  the 
region  must  therefore  have  finite  amplitude  at  its  front.  Such  a  dis¬ 
continuous  disturbance  is  a  shock  and  requires  separate  treatment. 

Consider  now  the  example  of  a  simple  wave  propagating  in  the 
positive  x  direction  into  a  constant -state  region  in  which  uq  =  0. 

Then 

u  -  a  s  -  a 

o 

and  the  two  equations  in  ( 3)  become  the  same, 

do  z  *  do 

St  +  (°  -  ao  +  c)  3SE  =  0 

which  has  the  general  solution 


or 


o  =  F[x  -  (cr  -  cr  +  c)t] 


o  =  F[x  -  (u  +  c)t] 


(8) 
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wherein  F  is  an  arbitrary  function  of  its  argument. 

As  an  example  of  the  use  of  this  solution,  consider  the  problem 
of  determining  the  motion  of  a  polytropic  gas  disturbed  by  a  piston. 

It  is  assumed  that  up  to  t  =  0  the  piston  is  at  rest,  so  that  a  semi¬ 
infinite  region  of  constant  state,  with  u  =  0,  a  =  d^,  has  been  estab¬ 
lished.  Subsequent  to  t  =  0,  the  piston  moves  with  prescribed  velocity 
o(t)  to  positions  x(t)  (such  that  v>  =  dx/dt).  For  a  polytropic  gas, 
a  =  2c/ (y  -  l),  so  that  (8)  can  be  written 


or 


and 


(9) 


Thus,  the  function  F  is  to  be  determined  by  substituting  into  (9)  the 
known  conditions  at  the  piston 


Try- *-  V)  t<o 

yrrr  *  »  *  p[x  -  (c0  +  2~T1  D)tJ  *  >  0 

and  F(0)  lies  between  2c .  /(y  -  l)  and  2cQ/(y  -  l)  +  V>. 


(10) 


i' 
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Example  I 

The  problem  of  withdrawal  of  the  piston  at  constant  speed  was 
partially  solved  earlier  in  the  chapter.  Here  we  may  solve  it  com¬ 
pletely  .  (The  same  problem  is  solved  more  easily  in  Chapter  V;  we  use 
this  method  here  for  illustration,  because  of  its  power  for  more  com¬ 
plicated  problems  which  cannot  be  treated  by  the  procedure  of  Chapter  V. ) 

In  this  case,  b  is  a  constant  (b  <  0)  and  x  =  bt.  From  (lO) 

2c 

F(|)  -  i  >  0 

2c 

F(§)  =  +  u  I  <  0 


where  5  is  any  argument  of  F. 

This  solution,  put  into  (9),  becomes 


or 


o  2c  ,  /  2c  \ 

— -  =  - -  for  x  -  - - -  l  c - r  ft  >  0 

y  -  1  7  -  1  7  -  1  \  7  +  V 


+  b  for  x  -  — -j-  ( c - I^r'V  <  0 

1  7  -  1  \  7  +  V 


7-  1 

7  - 1 

2c 

2c 

0 

7~  1 

7  “  1 

c  = 

c 

0 

2c 

2c 

0 

(11) 


for  x  >  c  t 
o 


=—  =  —  -_°1  +  b  for  x  <  (cQ  +  — b)t 


(12) 


The  first  statement  says  that  the  signal  propagates  with  sound  speed. 

The  second  one  states  that  between  the  piston  and  the  path  given  by 
x  =  (c^  +  ((7  +  l)/2]b}t,  the  sound  speed  is  a  constant,  and  the  result  is 
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identical  to  that  in  (6)  for  conditions  at  the  piston.  (Note  that  if 
■q  _  _  [2c  / ( 7  -  1 ) ] ,  the  escape  speed,  then  the  thickness  of  the  constant 
state  zone  next  to  the  piston  is  zero,  since  its  boundaries  both  move 

with  piston  speed. ) 

Finally,  in  addition  to  ( 1 1 )  we  have 


From  (9)  may  be  calculated  the  instantaneous  slope  of  the  function 


where  F'  means  the  derivative  of  F  with  respect  to  its  argument.  Thus, 
if  F(x)  has  negative  slope  at  any  time,  then  the  denominator  may  even¬ 
tually  vanish;  the  result  is  a  shock. 

As  a  specific  example,  consider  the  problem  of  determining  the 
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effect  of  a  uniformly  accelerating  piston,  moving  into  the  gas.  In  this, 
with  acceleration  a. 


=  - 2-  x  >  c  t 

y  -  1  o 

(The  condition  x  <  c  t  in  the  first  equation  follows  from  the  condition 
X-(7+l)/(7-l)  (c-  [2c  /(y  +  l)]i  t  <  0.  This  can  be  verified  by 
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noting  that  at  the  point  x  =  CQt,  c  =  cq  in  the  solution  that  follows.) 
The  equations  ( 1 7)  can  solved  for  c: 


c  =  c  x  >  c  t 

o  o  J 

(As  a  check  on  this  solution,  it  may  be  noted  that  at  the  piston,  where 
x  =  1/2  at2,  ( 1 8)  gives  2c/(7  -  0  -  2c J{y  -  1 )  +  at,  being  the  result 
obtained  previously  -  see  (6),  which  is  perfectly  valid  for  a  piston 
speed  varying  with  time . ) 

The  envelope  in  the  x — t  plane,  of  values  such  that  the  square 
root  vanishes  in  (l8),  is  a  path  of  particular  interest.  Along  it, 
dc/dx  is  infinite  so  that  the  path  defines  a  shock,  whose  position 

x  (t)  is  given  by 
s 

xs(t)  =  ~  cq2  +  (7  -  l)a  cQt  +  a2t2^  (19) 

This  is  valid,  however,  only  for  xg  <  cQt,  or  for  t  >  2cq/(7  +  1  )a. 
Hence  a  shock  forms  at  time  and  position 


and  proceeds  initially  with  speed  c  .  The  solution  cannot,  however,  be 
believed  after  the  initial  formation  of  the  shock,  since  thereafter  the 
flow  pattern  is  no  longer  that  of  a  simple  wave. 
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SHOCKS 

Description  of  a  Shock 

The  compressive-piston  problem  of  Chapter  III  shows  that  there 
are  circumstances  in  which  the  differential  equations  predict  the  flow 
field  of  a  fluid  to  develop  a  discontinuity,  or  shock.  So  far,  we  have 
not  discussed  any  procedure  for  handling  such  circumstances;  the  differ 
ential  equations  become  meaningless  at  a  shock  and  offer  no  clue  for 
subsequent  treatment  of  the  discontinuity  whose  formation  they  predict. 

The  concept  of  a  discontinuous  shock  is  very  useful,  even  though 
the  field  does  not  become  really  discontinuous.  Actually,  where  the 
differential  equations  predict  formation  of  a  discontinuity,  a  very 
thin  region  is  formed  over  which  the  field  variables  change  extremely 
rapidly.  The  structure  of  these  regions  is  discussed  in  Chapter  DC, 
where  it  is  shown  that  the  inclusion  of  viscosity  effects  in  the  dif¬ 
ferential  equations  removes  the  tendency  for  the  equations  to  predict 
discontinuities.  It  is  further  shown,  for  example,  that  shocks  in 
air  are  commonly  a  small  fraction  of  a  millimeter  in  thickness,  so 
that  in  a  study  of  interaction  with  large  bodies,  the  discontinuous 
approximation  may  be  very  good. 
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There  are  two  types  of  surface  discontinuities  which  may  occur* 

One  of  these,  called  a  "contact  discontinuity,"  moves  with  the  fluid. 

It  is  present  at  the  boundary  between  two  kinds  of  fluid,  and  it  may 
also  be  generated  within  one  fluid  under  certain  circumstances.  It  is 
characterized  by  the  continuity  of  pressure  and  of  normal  velocity  com¬ 
ponent  across  it. 

The  other  type  of  discontinuity,  the  shock,  moves  relative  to  the 
fluid,  changing  the  state  of  each  element  as  it  sweeps  by.  To  see  the 
manner  in  which  this  occurs,  it  is  sufficient  to  consider  a  relatively 
simple  case:  The  shock  is  an  infinite  flat  plane  separating  two  semi¬ 
infinite  regions  in  each  of  which  the  state  is  everywhere  the  same  and 
in  which  the  material  velocities  are  normal  to  the  shock  plane.  The 
results  will  be  instantaneously  appropriate  for  curved  shocks  with  non¬ 
constant  adjacent  states  as  long  as  the  material  speed  is  everywhere 
normal  to  the  shock.  The  reason  is  that  in  considering  a  shock  to  be 
of  zero  thickness,  the  vicinity  of  a  point  on  it  can  be  examined  micro¬ 
scopically"  to  a  sufficient  degree  to  give  the  vicinity  of  the  point 
the  idealized  conditions  described. 

We  use  the  notation  in  the  diagram;  the  shock  is  moving  to  the 
right  and  the  fluid  is  more  compressed 
on  the  left.  In  general,  the  shock 
speed,  v,  is  greater  than  the  material 
speed  on  either  side.  Also  u  >  u+ 

(for  a  shock  traveling  as  pictured), 
which  follows  from  the  fact  that  the 
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shock  is  compressive.  (We  shall  show  in  this  chapter  that  all  shocks 
are  compressive. )  Four  conditions  exist  which  relate  changes  of  these 
variables  across  the  shock.  One  of  these  is  the  equation  of  state. 

The  other  three,  called  the  Rankine-Hugoniot  relations,  are  derived  as 
consequences  of  the  three  fundamental  conservation  laws.  It  will  be 
useful  to  examine  several  derivations  of  these  relations. 

Fundamental  Derivation  of  the  Shock  Relations 

1  *  Conservation  of  Mass.  The  distance  per  unit  time  that  material 

moves  relative  to  the  shock  is  (v  -  u+)  on  the  right  and  (v  -  u_)on  the 

left.  The  mass  per  unit  area  of  shock  per  unit  time,  771,  which  passes 

into  the  shock  from  the  right,  is  thus  m  =  p(v-u).  This  must,  how- 

+  + 

ever,  be  exactly  the  same  as  the  mass  per  unit  time  which  leaves  the 
shock  to  the  left: 

m  =  p+(v  -  u+)  =  p_(v  -  u^)  (l  ) 

2.  Conservation  of  Momentum.  The  momentum  per  unit  area  per  unit 
time  passing  into  the  shock  from  the  right  is  m  u+;  similarly,  that 
which  leaves  from  the  left  is  Tti  u_  so  that  the  change  in  momentum  per 
unit  time  is  ttl  (u_  -  u+).  This  must  equal  the  force  per  unit  area  on 
the  system: 

m  (U_  -  u+)  =  p_  -  p+  (2) 

3.  Conservation  of  Energy.  The  energy  change  per  unit  area  per 
unit  time  is  m  (E_  -  E+).  This  must  be  equal  to  the  net  rate  of  doing 
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work  on  the  region: 

m (E_  -  E+)  =  p_u_  -  p+u+ 


(5) 


These  simple  results  will  be  extremely  useful  for  exploring  a  num¬ 
ber  of  interesting  hydrodynamic  situations.  For  convenience  in  handling 
the  shock  relations,  we  introduce  the  notation  for  each  field  variable 

1  (4) 

*  =  i  (♦  +  *_) 

where  stands  for  any  of  the  field  variables.  The  following  identity 
occasionally  will  be  useful: 

5(t^2)  ■  ti&*2  +  (5) 

In  terms  of  these  symbols,  the  shock  relations  can  be  put  into  the  form 


TTl  =  Vp  —  pU 

(6) 

v&p  =  &(pu) 

(7) 

?n  6u  =  &p 

(8) 

m  &E  =  5(pu) 

(9) 

Integral  Derivation  of  the  Shock  Relations 

In  Chapter  I,  it  was  shown  that  the  equations  could  be  expressed 
in  integral  form  (l-2 6),  (l-27).  We  now  assume  that  the  equations  thus 
written  are  valid  even  in  the  presence  of  a  shock.  A  posteriori  justi¬ 
fication  of  this  assumption  will  follow  from  the  fact  that  the  results 
of  the  following  analysis  are  correct. 
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The  line  in  the  diagram  is  a 
plot  of  shock  position  versus  time, 
with  the  shock  moving  in  the  posi¬ 
tive  x-direction.  Draw  a  rectangle 
about  a  segment  of  the  curve  in  the 
vicinity  of  a  point  of  interest. 

The  rectangle  is  to  be  considered 
so  narrow,  in  a  direction  normal 
to  the  curve,  that  contributions 
to  an  integral  around  the  rectangle  due  to  its  ends  will  be  negligible. 
The  rectangle  cuts  the  shock  curve  at  (x^,t^)  and  ( x^ , t^ ) .  The  integral 
equation  for  a  path  around  the  rectangle  is  thus  (see  page  30) 

A_(x2  -  x^  -  B_(t2  -  t^  +  A+(x1  -  Xg)  -  B+(t1  -  t2)  =  0 

(where  higher  order  contributions  due  to  finite  rectangle  size  have 
been  neglected) .  If,  now,  the  rectangle  is  sufficiently  short,  then 
(x2  -  Xl)/(t2  -  t^ )  differs  negligibly  from  the  speed  of  the  shock,  v. 
Thus,  the  result  from  the  integral  equation  is 

v6A  -  5B  =  0  (10) 

From  the  conservative  Eulerian  equations,  (1-26)  we  obtain 
vSp  =  8(pu) 

vS(pu)  =  8(pu2  +  p)  (ll) 

v5(pE)  =  b( pEu  +  pu) 


DISCUSSION  OF  SHOCK  FORMATION 


The  first  of  these  is  identical  with  the  previous  result  (7).  It  can  he 
shown  that  all  three  of  them  are  reducible  to  the  results  obtained  from 
the  fundamental  derivation. 

From  the  conservative  Lagrangian  equations,  ( 1-27) ,  we  obtain 

v'p  &(— )  =  -  5u 
Ko  'p' 

v 1 pQ5u  =  6p 

v'PqSE  =  S(pu) 

where  v'  is  the  Lagrangian  shock  speed  (that  is,  tne  speed  of  the  shock 

relative  to  the  unshocked  material).  The  quantity  v'p  is  just  equal 

o 

tom,  the  mass  per  unit  area  per  unit  time  passing  across  the  shock. 

Thus 

7TX,5(1)  =  -  5u 

M5u  =  5p  y  (12) 

M6E  =  &(pu) 

J 

Through  the  Eulerian  definition  of  m  in  ( 6) ,  it  can  be  shown  that  ( 1 1 ) 
and  (12)  are  equivalent.  The  fundamental  derivation  produced  the  first 
equation  of  (ll)  and  the  second  two  of  (12).  This  auxiliary  integral 
method  of  derivation  produced  the  other  forms  shown,  which  will  be 
useful. 

Discussion  of  Shock  Formation 

It  was  shown  in  Chapter  III  that  the  compressive  motion  of  a  pis¬ 
ton  is  likely  to  cause  a  shock.  Here  we  examine  in  more  detail  the 


IV.  SHOCKS 


basis  for  this  peculiar  phenomenon.  There  are  two  processes  which  con¬ 
tribute  to  the  formation  of  a  shock.  Consider,  for  example,  the  com¬ 
pressive  motion  of  a  piston.  In  the  pulse  adjacent  to  the  piston,  the 
fluid  is  moving  towards  the  undisturbed  fluid  which  is  further  away. 
Signals  from  the  piston,  indicating  its  approach,  are  propagated  with 
sound  speed  relative  to  the  moving  material.  Thus,  relative  to  a  fixed 
reference  frame,  signals  approach  the  pulse  front  faster  than  it_  moves 
through  the  undisturbed  fluid.  These  signals  pile  up,  producing  the 
shock.  The  second  process  which  also  contributes  to  the  shock  forma¬ 
tion  is  the  heating  of  the  pulse  adjacent  to  the  piston.  This  heating 
also  increases  the  sound  speed  over  that  in  the  undisturbed  region. 

Both  of  these  processes  are  manifested  mathematically  by  the  non¬ 
linear  character  of  the  differential  equations.  Thus,  for  example,  in 
the  simple-wave  equation  of  Chapter  III, 

da  ,  ,  „ 

3t  +  (a  '  °o  +  o)  S'  0 

the  nonlinearity  arises  from  the  presence  in  the  parentheses  of  a  (pro¬ 
ducing  the  first  shock-forming  process)  and  of  c (a)  (producing  the 
second  shock-forming  process).  If  we  linearize  the  equation  by  putting 
a  =  and  c  =  cq,  the  general  solution  is  that  of  a  propagating  pulse 
whose  shape  does  not  change  with  time 

O  =  F(x  -  CQt) 

However,  the  more  general  solution,  as  shown  in  Eq.  ( III-1 4) ,  can  lead 
to  shock  formation. 
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Usually,  the  formation  of  a  shock  makes  the  theoretical  treatment 
of  a  problem  much  more  difficult.  The  shock  is  a  moving  boundary,  along 
which  boundary  conditions  can  be  supplied,  but  whose  motion  with  time  is 
not  a  priori  known.  In  addition,  the  fluid  may  no  longer  be  isentropic 
(as  will  be  shown  later).  As  a  result,  analytical  solutions  in  the  pres¬ 
ence  of  variable  shocks  are  veiy  difficult  to  obtain  except  in  certain 
special  cases  (see,  for  example.  Chapter  Vi). 


Shock  Relations  for  Special  Cases 

Probably  the  three  most  useful  of  the  general  forms  of  the  shock 
relations  (obtainable  by  putting  E  =  I  +  l/2  u  and  eliminating  nvfrom 
the  shock  relations)  are 

v&p  =  &(pu) 

(&p)(&  £)  =  -  (6u)2  03) 

51  =  -  p5(£) 

In  particular,  if  the  coordinate  system  is  such  that  u+  =  0,  then 
v(p  -  p  )  =  p  u  0*0 


2 

(p_  -  p+)(p_  -  p+)  =  p_p+u. 


p  +  ph 

J-  -  J+  ^1“ 


(p  -  p  0 


05) 


06) 


In  the  special  case  that  the  gas  is  polytropic,  these  equations 
can  be  rewritten  in  a  number  of  convenient  forms.  Let  c_  and  c+  be  the 
sound  speeds  behind  and  ahead  of  the  shock,  respectively,  and  define 
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P_ 


Relative  to  the  gas  ahead  of  the  shock,  M  and  U  are,  respectively,  the 
Mach  numbers  of  the  shock  and  of  a  piston  producing  it.  M'  is  the  Mach 
number  for  the  flow  behind  the  shock.  It  follows  from  the  shock  rela¬ 
tions  that  U,  Z,  P,  and  M'  are  all  functions  of  M  alone.  Thus 


or,  conversely. 
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If  the  shock  is  very  strong  (that  is,  the  shock  speed  is  large 
compared  with  the  sound  speed  ahead)  then  M  -» °°  and 

m  r  +  1 


M'  -» 

Thus,  for  an  infinite-strength  shock, 

1 .  The  shock  speed  is  determined  by  the  gas  speed  behind.  If 
the  shock  is  formed  by  a  piston  moving  with  unifomn  speed,  then  the  ma 
terial  speed  behind  the  shock  equals  the  piston  speed  and 


2.  The  compression  produced  by  the  shock  is  independent  of  the 
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shock  speed  and  depends  only  on  the  nature  of  the  gas.  Notice  from  ( 1 8 ) 
that  the  greatest  possible  compression  that  can  be  produced  by  a  shock 
in  a  polytropic  gas  is  (7  +  1  )/{y  -  l).  For  air,  for  example,  7  =  1.4 
so  that  the  greatest  possible  shock  compression  in  air  is  6. 

5.  The  maximum  Mach  number  for  the  flow  behind  a  shock  moving  into 
a  polytropic  gas  at  rest  is  a  function  of  the  value  of  7  only.  For  air, 
this  limit  is  1.89* 

The  relations  for  an  infinite  strength  shock  can  be  summarized  in 
more  convenient  forms 


(25) 


Note  from  the  third  of  these  that  there  is  equal  partition  between  inter¬ 
nal  and  kinetic  energy  behind  an  infinite  strength  shock.  Note  also  that 
for  most  real  gases,  which  often  can  be  considered  polytropic  with  7  <  2, 
the  material  speed  behind  the  infinite  stre  ;h  shock  is  greater  than  the 
sound  speed. 


-z  =  Z  +  1 
p+  7-1 


v  = 


I  = 


-  y  t  1 


u 


i  u2 


2  ( 7  —  1  \  2/2^ 


u 


2 

7(7  -  1 ) 


-80- 


entropy  change  across  a  shock 


Entropy  Change  Across  a  Shock 

The  entropy  change  across  a  shock  in  a  polytropic  gas  is,  according 


to  (II-5), 


where  h  is  the  specific  heat  at  constant  volume.  Combining  with  the 
general  shock  relations,  we  may  put  this  in  the  form 


-  S  =  b 
+ 


^0  "  7  ftf)  -ln  (j  +  V i  +  r) 


-y  In  R 


(26) 


If  R  =  1  (no  density  change  across  the  shock),  then  S_  =  S+.  For  R  only 
slightly  different  from  1 ,  then  (26)  may  be  expanded  to  give 


S.  -  S+  ~  V-d  (R  -  1  ?  (27) 

where  higher  power  terms  in  (R  -  1 )  have  been  dropped. 

The  Second  Law  of  Thermodynamics  states  that  entropy  increases  or 
is  constant  in  any  closed  macroscropic  system.  Thus,  when  a  shock  is 
formed,  R  must  be  greater  than  1 .  Furthermore,  for  any  shock,  R  must 
always  be  greater  than  1  since,  to  change  to  R  ^  1,  the  region  of  valid¬ 
ity  of  (27)  must  be  approached.  This  means  that  shocks  must  always  be 
compressive  in  the  polytropic  gas  discussed  here.  A  more  general  argu¬ 
ment  along  similar  lines  can  be  used  to  prove  that  all  shocks  are 
compressive. 
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Hugoniot  Relation 

There  are  two  extreme  ways  (with  no  heat  flow)  in  which  a  fluid  can 
be  compressed  monotonically:  impulsively  or  very  slowly.  In  the  latter 
case,  the  fluid  remains  isentropic,  and  the  adiabatic  equation  of  state 
holds.  Alternately,  in  the  case  of  shocks,  we  have  seen  (l 6)  that 


P 


This,  together  with  the  equation  of  state  in  the  form  I  =  l(p,p),  gives 
a  relation,  called  the  Hugoniot  relation,  between  pressure  and  density 
for  shock  compression. 


A  plot  of  the  adiabatic  and 
Hugoniot  pressure-density  relations 
for  a  material  shows  that  the  "adia- 
bat"  lies  below  the  "Hugoniot."  This 
is  because,  in  general,  a  shock  com¬ 
pression  increases  the  internal  energy 
nore  than  an  adiabatic  compression, 
and  further,  for  a  fixed  compression, 
pressure  increases  with  internal 
energy. 

For  any  monotonic  compression 
of  a  fluid,  its  state  will  always  be 
such  that  its  pressure  and  density 
determine  a  point  which  lies  between 
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the  A  and  H  curves.  /  Any  subsequent  expansion  of  an  element  of  the  fluid 
will  be  adiabatic.  In  general,  the  slope  of  the  adiabat  at  any  point 
is  less  than  the  slope  of  the  Hugoniot.  As  a  result,  it  is  impossible, 
in  the  absence  of  heat  flow,  for  a  fluid  to  be  changed  in  such  a  way 
that  its  (p,p)  point  will  lie  directly  to  the  right  of  any  point  where 
it  has  previously  been.  This  irreversible  trend  of  migration  of  (p,p) 
points  to  the  left  is  identically  equivalent  to  the  trend  of  irreversibly 
increasing  entropy. 
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DIMENSIONALITY  AND  SIMILARITY 

Scaling 

Consider  the  problem  of  determining  the  flow  field  about  an  air¬ 
plane.  We  assume  for  now  that  the  air  can  be  treated  as  a  simple  gas 
(i.e.,  the  effects  of  viscosity  and  heat  conduction  can  be  neglected.) 
Even  so,  the  analytical  calculation  would  be  extremely  difficult,  and 
an  experimental  approach  is  required.  A  model  is  built  and  tested  in  a 
wind  tunnel.  The  question  then  arises:  What  bearing  do  the  test  re¬ 
sults  have  on  the  nature  of  flow  around  the  real  airplane?  Another 
question  of  interest  is:  Can  the  results  of  a  test  at  the  extreme  of 
wind-tunnel  capabilities  be  extended  to  give  information  in  the  region 
beyond  those  capabilities?  Considerable  effort  has  been  expended  to 
answer  these  questions,  and  much  of  many  aerodynamic  textbooks  is  de¬ 
voted  to  related  discussions.  We  shall  discuss  the  matter  only  slightly. 

Consider  the  Eulerian  equations  in  the  form 

P  +  p(u*V)u  -  -  Vp 
+  V.(pu)  =  0 
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p  ^  +  p(u-v)l  =  -  pW 

Now  let  the  entire  system  be  rescaled  such  that 
P  -PQP 


SCALING 


P  »PGP 

I  -» I  I 
o 

~4  — 5 

u  —*  u  u 
o 


t  -» t  t 
o 


r  -♦  x  r 
o 

where  the  subscript  o  here  indicates  a  constant,  dimensionless  scaling 
factor.  Then  the  equations  become 

2 


©*>£  *  (r©> 


These  equations  will  be  identically  the  same  as  the  original  ones,  pro¬ 
vided  that 


P  x 
*o  o 


2  u  t  ;2 
d  u  o  o  u 
'o  o  o 


=  1 
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p 

In  addition,  the  sound  speed  is  defined  by  c  =  (dp/dp)  so  that,  with 

s 

c  -*  c  c,  we  also  have 
o 


and  the  scaling  laws  can  be  rewritten 


o 


Thus,  in  the  case  of  the  airplane  and  its  model,  results  of  tests  can  be 
applied  to  the  full  scale  provided  that  the  ratio  of  flow  speed  to  sound 
speed  is  the  same  in  both  cases.  This  important  ratio  is  called  the  Mach 

p 

number  of  the  flow.  The  requirements  Iq  =  Pq/po  =  c^  are  certainly  sat¬ 
isfied  for  a  polytropic  gas.  The  relation  x  /t  =  c  gives  the  time 

o  o  o 

scale  for  occurrences  in  tenons  of  the  scaling  of  linear  dimensions.  Thus, 
if  the  sound  speed  is  the  same  for  the  airplane  and  its  model  and  if  the 
airplane  is,  say,  5  times  as  long,  then  it  takes  5  times  as  much  time  for 
a  similar  change  in  flow  pattern  to  take  place. 

When  the  gas  is  viscous  and  conducting  of  heat,  additional  param¬ 
eters  result  from  scaling  arguments.  These  are  discussed  in  Chapter  VIII. 


Similarity 


In  some  cases  of  interest,  a  flow  field  may  retain  geometric  simi¬ 
larity  to  itself  as  it  develops  in  time.  This  is  the  case  in  one  dimen¬ 


sion,  for  instance,  if  all  features  of  the  flow  field  are  functions  of 

x/t  alone.  The  example  of  piston  withdrawal  in  Chapter  III  [see  (ill-l^)] 
is  an  illustration. 
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Consider  the  one-dimensional  equations  in  the  form 


du  du 

p  St  +  pu  SE  = 


c)A  .  dA  a 

3t  +  us;  =  0 


where  A  is  any  function  of  entropy  alone;  for  a  polytropic  gas,  we  take 


We  now  assume  that  there  exists  a  function  y(x,t)  such  that  p,  u,  and 
p  are  functions  of  y  only.  Then 


and  similarly  for  p  and  p.  Here  a  prime  denotes  differentiation  with 
respect  to  y.  Note  that  this  is  really  a  part  of  the  more  general 
transformation 

y  =  y(x,t) 
z  =  t 

after  which  we  assume  independence  from  z,  so  that  all  features  of  the 


-87- 


V.  DIMENSIONALITY  AND  SIMILARITY 


flow  field  are  independent  of  time  in  the  transformed  coordinate  system. 
Note  that  the  transformed  coordinate  system  moves  with  respect  to  the 
original  one  with  a  velocity 


v 


(dy/dt) 

(dy/dx) 


(2) 


Putting  these  assumptions  into  the  original  equations,  we  get  expressions 
which  reduce  to 

pu'(u  -  v)  =  -  p' 

p'  (u  -  v)  =  -  pu'  (3) 

A*  (u  -  v)  =  0 


Thus,  to  be  consistent  with  the  original  assumptions,  v  must  be  a  func¬ 
tion  of  y  only.  This  is  equivalent  to  the  condition 

y  =  f[x  -  v(y)t]  (4) 

where  f  is  an  arbitrary  function  of  its  argument.  Thus,  in  this  general¬ 
ized  similarity  transformation,  f  and  v  can  be  chosen  as  arbitrary  func¬ 
tions  and  the  equation  for  y(x,t)  cam  be  obtained. 

The  special  case  in  which  v  =  u  is  trivial.  It  is  an  equilibrium 
solution  in  which  all  characteristics  of  an  element  of  fluid  are  con¬ 
stant.  It  can  be  written  in  the  form 

A(y)  s  arbitrary  function 

u(y)  =  constant 

p(y)  s  constant 

p(y)  =  arbitrary  function 
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and  I(.y)  is  related  to  p(y)  in  such  a  way  that  the  pressure  is  everywhere 


the  same. 


If  v  £  u,  then 


A  s  constant 


(pu')2  =  p'p' 

Thus  p  is  a  function  of  p  only,  and  we  may  use  the  definition  (III-2) 


to  obtain 


du  =  ±  dcr 


u  =  u  ±  a 
o 


where  u  is  a  fixed  constant.  Note  that  this  is  the  same  result  as  ob- 
o 

tained  from  the  theory  of  characteristics  except  that  now  there  is  not 
yet  any  specification  of  the  path  along  which  u  ±  or  is  constant.  We  put 
(5)  into  (3)  to  obtain 

±  pa’(uQ  ±  a  -  v)  =  -  cp<J * 

The  case  o'  =  0  is  again  trivial;  in  order  that  o'  o,  it  is  necessary 


u  ±a-v±c  =  0 
o 


v  =  u  £  c 
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Thus  the  transformed  coordinate  system  must  everywhere  move  with  the 
velocity  of  the  characteristic  lines. 

The  problem  of  the  constant- speed  retracting  piston  of  Chapter  III 
is  now  solved  with  ease.  Put  v  =  a(x/t)  (where  a  =  constant)  and,  for  a 
polytropic  gas,  cr  =  2c/(y  -  l).  Then 

Q£  ^  =  u  ±  c 
X/ 


Application  of  appropriate  boundary  conditions  then  yields  the  solution. 

Various  other  types  of  similarity  solutions  are  applicable  to 
special  cases.  Sometimes  considerable  ingenuity  is  required  to  find  one 
which  satisfies  the  desired  initial,  and  boundary  conditions.  In  many 
cases,  no  such  solution  has  yet  been  found.  Chapter  VI  shows  the  appli¬ 
cation  of  another  type  of  similarity  assumption  to  the  problem  of  a  de¬ 
caying  shock. 
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DECAY  OF  A  SHOCK  PULSE 


General  Discussion 

We  have  seen  in  previous  chapters  how  the  compressive  motion  of  a 
piston  can  produce  a  shock.  In  most  practical  situations,  however,  the 
piston  pushes  for  only  a  finite  time.  When  it  stops,  the  fluid  at  first 
continues  to  rush  away  and  a  rarefaction  is  formed.  In  general,  the 
rarefaction  will  overtake  the  shock  and  modify  it;  the  shock  will  decay. 
The  sequence  is  shown  qualitatively  in  the  following  successive  plots  of 
pressure  as  a  function  of  distance. 


1 


/ 


i 


=1 
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Most  problems  concerning  shock  decay  are  very  difficult  to  solve. 
We  shall  consider  two  special  cases.  In  the  first  one,  the  shock  is 
formed  by  the  very  short  impulsive  motion  of  a  piston  adjacent  to  a 
cold  polytropic  gas.  In  the  second  one,  a  weak  shock  decays  under  the 
action  of  an  overtaking  simple  rarefaction. 

Decay  of  a  Very  Strong,  Impulsively-Formed  Shock  Pulse 

In  this  section  we  discuss  a  special  case  of  a  more  general  result 
as  summarized  from  the  Russian  literature  by  Hayes  and  Probstein.  Their 
calculations  were  performed  in  Eulerian  coordinates;  we  have  used 
Lagrangian  coordinates  and,  in  addition,  have  tried  to  show  how  the  var¬ 
ious  similarity  assumptions  arise  from  one  basic,  initial  assumption. 

Because  we  shall  refer  to  both  Eulerian  and  Lagrangian  coordinates, 

we  shall  carefully  keep  track  of  them  by  the  designations  x  and  x  ,  re- 

o 

spectively  (at  t  =  0,  x  s  xq).  The  problem  is  described  as  follows. 
Before  t  =  0,  the  polytropic  gas  is  cold  and  at  rest  (thus  it  has  zero 
sound  speed).  At  t  =  0,  an  instantaneous  impulse  is  applied  by  a  pis¬ 
ton  at  x  =  0;  thereafter  the  piston  remains  fixed  at  x  =  0  so  that  the 
amounts  of  energy  and  mass  in  the  region  for  x  >  0  remain  constant. 

A  shock  propagates  away  from  the  piston;  its  position  is  x  which  is  a 

s 

function  of  time,  and  its  velocity  is 


W.  D.  Hayes  and  R.  F.  Probstein,  "Hypersonic  Flow  Theory,"  Academic 
Press,  Inc.,  New  York,  1 959 >  PP*  52  ff. 
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The  boundary  conditions  at  the  shock  are 


=  2. 


+  1 


7  ~  1 


u 

s 


2 

7  +  1 


v 


P 


s 


2 

7  +  1 


Thus,  the  density  is  forever  constant  at  the  shock,  while  the  pressure 
and  material  speed  vary  with  the  shock  speed. 

We  work  with  the  Lagrangian  equations  of  motion 


The  basic  assumption,  which  will  be  given  a  posteriori  justifica¬ 
tion,  is  that  there  is  a  generalized  similarity  solution  in  which  the 
density  is  a  function  of  x^/x^  only.  That  is, 

p(xQ,t)  =  Pgf(y)  0) 


where 


x 

o 


From  this  assumption  only,  there  result  the  other  features  of  the 


(2) 


-93- 


VI.  DECAY  OF  A  SHOCK  PULSE 


DECAY  OF  A  VERY  STRONG,  IMPULSIVELY-FORMED  SHOCK  PULSE 


so  that  g( 1 )  =  1 ,  and  the  mass  equation  becomes 


2  2 

yf’  =  f  s' 


(5) 


7  -  1 

Next,  consider  the  momentum  equation.  By  a  similar  process,  it  becomes 


2 _  (  dv  2  A 

—  p0  Vxs  dt s  - v  y 6  7 =  " 


Thus,  p,  considered  as  a  function  of  y  and  t,  must  be  written  in  the 
form 

p(y,t)  =  x  ^(function  of  y)  +  v  (function  of  y) 

S  QX 


But 

Pd.t)  S7TT  V2 

2 

Thus  x  (dv/dt)  must  be  identically  proportional  to  v  ,  and  it  follows 
that  we  can  write 

x  =  Ktn  (4) 

s 


where  K  and  n  are  arbitrary  constants.  (This  comes  from  integrating 
the  equation 


and  applying  the  boundary  condition  that  at  t  =  0,  x  =0.)  Thus,  we 

s 

put 

P  =  — P  h(y) 

*  y  +  1  *o 
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so  that  h(  1 )  =  1 .  As  a  result,  the  momentum  equation  has  become 


,  ,  .  1  -  n 

h'  =  yg*  +  1  -n-  ■  ■  £ 


Finally,  there  is  the  entropy  equation  which  states  that 

0 

Thus  the  quantity 

v\(y) 

f7(y) 

must  be  independent  of  time.  Because  v^  is  proportional  to  t  ,  we 


see  that 


h  (2n-2)/n 

-r  ~  y 


Since  h(  1 )  =  f( 1 )  =  1 ,  the  constant  of  proportionality  is  1 ,  and 


h  -  f 


The  over-all  conservation  laws  give  further  information  concerning 
the  problem.  Thus,  for  over-all  conservation  of  mass,  we  require 


J  pdx  =  p  x 
o  os 


This,  however,  is  satisfied  identically,  since 


!q  pd x=JQ  pQdxo 
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The  momentum  equation  gives  more  valuable  results.  We  require 


dt  J q  pu  dx  =  p(O) 


where  p(O)  is  the  pressure  at  the  point  x  =  0,  that  is,  at  the  piston. 


4  r» ax  .  EM 
dt  Jo  °  po 


L  8dy " h(0) 


From  this,  we  may  distinguish  several  classes  of  solutions.  If  n  =  l/2, 
then  h(o)  =0.  If  n  /  l/2,  then  h(o)  /  0,  or  else  if  h(0)  =  0,  then  g 
must  change  sign  in  the  interval  or  identically  vanish.  We  shall  see 
below  that  n  =  2/5  and  thus  conclude  that  h(o)  /  0  is  required. 

The  energy  conservation  equation  states  that 


f  (y*r 


+  5  pu  j  dx  =  Ej, 


where  ET  is  the  constant  total  energy  of  the  system.  This  can  be  put 
into  the  form 


Jo  (J  +  <?) 


(7  +  1 

2 

2x  v  p 
s  Ho 


Thus  x  v  must  be  constant  in  time.  Since 
s 
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2  v-3  2  3n-2 

x  v  s  K^n  t 
s 


it  follows  that  n  =  2/3. 

In  summary,  with  n  =  2/3, 


we  now  have  the  three  equations 


h'  =  yg'  +  \  g 
yf 1  =  y-§T  f2g’  ” 

yh  =  f 7 


together  with  the  conditions  h( 1 )  =  f ( 1 )  =  g( 1 )  =  1 ,  and  the  expectation 
that  h(0)  /  0,  g(0)  =  0,  and  that  in  the  interval  0  <  y  <  1  the  depen¬ 
dent  variables  are  all  positive.  If  these  conditions  can  be  fulfilled, 
then  the  similarity  assumption  may  be  reasonable.  It  may  be  noted  that 
f(o)  =  0,  so  that  the  temperature  at  the  origin  must  be  infinite  in  or¬ 
der  that  the  pressure  not  vanish.  This,  then,  adds  another  physical  con¬ 
dition  to  the  limit  of  applicability. 

Once  a  solution  has  been  obtained,  then  the  energy  integral  can  be 
evaluated  and  K  determined  as  a  function  of 


dy  = 


9(7  +  1 )% 


8k\ 


The  left  side  is  a  function  of  y  only,  so  that 


(10) 


l/3 


(function  of  y) 


(11) 


One  useful  relationship  which  may  be  obtained  without  a  complete 
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solution  is  that  between  the  Eulerian  and  Lagrangian  position  variables. 
We  use  the  relation 

pdx  =  p  dx 
o  o 


from  which 


o 


Now 


yf _  2g* 

2  ~ 

f  7-1 


so  that 


/ 


dy  =  Mil 

o  /  7  7-1 


Integrate  this  by  parts  and  combine  with  the  equation  for  x  to  obtain 


__  ,  ag(y) 

7  +  1  s|_7  -  1 


7  -  1 

X  =  - - r  X 


f&TJ 


(12) 


This  shows  that  x/xg  is  a  function  of  xQ/xg  only,  so  that  h(y),  f(y), 

and  g(y)  could  just  as  well  have  been  written  as  functions  of  x/x  . 

'  s 

If  the  solution  is  obtained  in  one  coordinate  system,  it  can  thus  be 
transformed  to  the  other. 
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The  process  of  solving  the  equations  can  be  simplified  by  using 
the  fact  that  there  must  be  no  energy  flux  across  any  similarity  line. 
Thus 


where  us^m  is  the  velocity  of  a  similarity  line  at  fixed  y.  From  the 
relationship  between  Eulerian  and  Lagrangian  coordinates,  we  may  see 
that 


so  that  the  desired  relation  is 

(h  +  fg2)y  =  2hgf  (14) 

This  may  be  used  in  place  of  any  of  the  differential  relations. 

Hayes  and  Probstein  display  the  complicated  solution  to  the  prob¬ 
lem  and  show  graphs  of  the  results  (as  a  function  of  the  Eulerian  simi¬ 
larity  variable)  for  y  =  1.4.  We  here  derive  instead  a  simple  approxi¬ 
mate  solution  which  shows  the  qualitative  features  and  which  satisfies 
the  boundary  conditions.  We  start  by  trying  as  a  solution 

h(y)  =  h(0)  +  >iyn 

s(y)  =  ym 

where  m  and  n  are  to  be  determined.  To  satisfy  h(l)  =1,  A  =  1  -  h(0). 
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These  satisfy  the  equation  which  does  not  contain  f  provided  that 

n  =  2(1  -  A) 

2A  -  1 
m  =  2(i~n^ 


The  equation  between  g  and  f  can  also  be  satisfied  provided  that 


1-m 


1-m  2m  /•  t 

f  +  (7  -  l)(l  -"m)U  “  7  ' 


Finally,  however,  there  is  a  relation  between  h  and  f  which  cannot  be 
satisfied,  even  to  lowest  order  in  y  unless  7  =  1  (in  which  case  it  can 
be  satisfied  identically) .  ’  In  order  that  the  lowest  order  exponents  be 
the  same,  we  choose  7(1  -  m)  =  1.  As  a  result,  there  follows  the  approxi- 
mate  solution 


g(y)  ■  ^)h 

h(y)  =  2T2r±rT)Lr  +  (5r  '  2)yl 

11  h 


(27-1 )/7 


f(y)  = 


yh(y) 


05) 


This  predicts  that 

h(0)  =  2(277-  1) 

which  is  surprisingly  accurate.  The  exact  result  is 


(16) 
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g 


A  comparison  is  shown  of  the 
approximate  and  exact  solutions  for 
g(y)  for  the  case  y  =  5/3 •  For 
smaller  values  of  y ,  the  agreement 
is  even  better,  and  in  the  limit 
y  ~¥  1  the  approximate  solution  be¬ 
comes  exactly  correct. 


It  is  an  important  feature  of  the  solution  that  the  shock  pressure 


decreases  in  proportion  to  t  '  J  in  this  constant -energy  case. 


Interaction  of  a  Rarefaction  Overtaking  a  Weak  Shock 

Consider,  first,  a  shock  moving  into  a  polytropic  gas  which  is  at 
rest  but  which  has  sound  speed  c+.  We  refer  to  equations  (lV-17)  and 
(lV-l8),  put  into  the  form 


u-  2(M^  -  1) 
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^  =  — — — r*  /[ 7  -  1  +  (2/M2)]  [27M2  -7+1] 
c+  7  +  1  V 

We  shall  need  these  in  the  weak-shock  limit;  let  M  =  1  +  e,  where  e  «  1; 
then 


c+  7  +  1 


(18) 


=>  1  +  2€ 


<M) 


This  shock  is  followed  hy  a  simple  rarefaction  which  arises  at 


x  =  0,  t  =  0.  Using  the  method  of  Chapter  III,  we  find  that  the  rarefac¬ 
tion  can  he  characterized  by 


within  the  rarefaction  region.  The  constant  K  is  related  to  the  shock 
conditions  at  the  instant  of  collision  by 

•^C-  7-1 

—  -  r— -  u 

at  collision 


K 


■  (^rr  -  fTT  0 


The  question  now  arises:  Can  the  rarefaction  solution  be  joined  identi¬ 
cally  to  that  of  the  shock?  In  general,  the  answer  is  no.  Only  in  the 
special  case  that  the  shock  is  weak  can  the  Joining  be  made: 


VI.  DECAY  OF  A  SHOCK  PULSE 


where  xg  is  the  shock  position.  These  two  relations  must  he  identically 
the  same,  if  the  joining  is  correct;  with  the  general  shock  relations 
this  cannot  he  accomplished;  in  this  special  case  it  can,  as  follows. 

The  two  equations  are  written 


7  +  1 


-r  -  2e  c  = 
t  + 


K 


x 


s 


t 


-  2e  c 

+ 


7  +  1 
7  -  1 


These  are  the  same  provided  that 
2c 


K  = 


7  +  1 


In  the  weak-shock  limit  this  requirement  is  satisfied.  Thus  the  two 
equations  "become 


f  -  2€  C+  =  C+ 


Now  €  =  M  -  1  so  that  c+e  a  v  -  c+.  Also, 
dx 


dt 


so  that 


x  /dx 

_  2  ( — S  -  c 
t  Vdt 


or 
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This  has  the  solution 

x  =  c  t  +to-/t  (^9) 

s  + 

in  which  ^  is  a  constant.  Thus,  the  shock  speed  asymptotically  approaches 
sound  speed  ahead  of  it: 

v=c  +JL  (20) 

+  2  Vt 

Note  that  this  solution  predicts  infinite  shock  speed  at  t  =  0,  in  con¬ 
tradiction  to  the  weak-shock  assumption.  Thus  the  solution  is  valid  only 
for  late  times.  Note  also  the  prediction  of  shock -pres sure  decay  as 
p  ^ /2  compared  to  t  in  the  preceding  section. 


-105- 


CHAPTER  VII 


INVERTED  HYDRODYNAMIC  EQUATIONS 


The  Transformation 

Even  for  a  simple,  isentropic  gas,  the  hydrodynamic  equations  are 
difficult  to  handle  in  complete  generality.  The  complications  arise 
mainly  from  nonlinearity  of  the  equations.  One  means  of  circumventing 
the  difficulty  is  by  the  method  discussed  in  this  chapter  (sometimes 
called  the  "speedgraph"  method) .  The  hodograph  transformation  for 


steady  flow  is  part  of  a  very  similar  method. 

We  start  with  the  equations  in  the  form  used  in  Chapter  III  (see 
III-3);  the  gas  is  assumed  to  be  simple  and  isentropic. 


du 

du  A 

da 

5t  + 

U£E  + 

it 

a 

do 

do 

du 

+ 

(D 


To  remove  the  nonlinearity  these  are  transformed  to  a  set  of  equations 
in  which  the  dependent  variables  are  x  and  t  and  the  independent  var¬ 
iables  are  u  and  a;  the  resulting  equations  are  then  linear  and  homo¬ 
geneous,  and  thus  amenable  to  treatment  by  more  familiar  methods. 

We  introduce  the  shorthand  notation  xff  =  ( dx/da )u=cons^an-(: >  etc.. 
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Det  s  x  t  -  x  t 
o  u  u  of 


(3) 


Validity  of  the  transformation  requires  Det  4  0* 

With  these  transformation  equations,  the  equations  in  (l)  become 


ax 

dt 

at 

^  ' 

u  35  = 

• c  3^ 

ax 

dt 

at 

5u  ’ 

- 

w 

(5) 


Since  c  is  a  function  of  a  only,  these  equations  are  linear  and 
homogeneous  in  their  dependent  variables.  From  these  two  equations, 
x  can  be  eliminated 


d^t  8^t  _  1  /  dc ,  dt 

Su2  ‘  So2  =  “  +  do  35 


(6) 


A  similar  equation  for  x  can  also  be  obtained,  but  it  is  not  as  concise. 
The  amount  of  difficulty  involved  in  solving  (6)  depends  upon  the  nature 
of  the  function  (l/c)[  1  +  (dc/dff)],  which,  in  turn,  depends  upon  the 
form  of  the  equation  of  state.  For  a  polytropic  gas,  for  example, 


a2t  a2t  (y  +  i\  i  at 
au2  aa2  ^  -  v  ff  35 


(7) 


Even  simpler  is  the  equation  resulting  from  the  equation  of  state 

p  =  a-£  (8) 

where  a  and  (3  are  constants.  [Such  an  equation  of  state  may  be  useful 
if  a  very  small  range  of  densities  is  involved;  then  (8)  may  fit  the 
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true  equation  of  state  sufficiently  well  over  the  range  of  interest.] 

The  form  of  (6)  with  the  relation  (.8)  is 

d2t  =  o 

7~2  ",  2 
ou  8a 

the  simple  wave  equation,  which  has  the  general  solution 
t  =  i  (u  +  o)  +  l2(u  -  a) 

where  l 1  and  t 2  are  arbitrary  functions  of  their  arguments  (arbitrary, 
that  is,  except  for  the  restrictions  imposed  by  continuity  and  nonvanish¬ 
ing  of  the  transformation  determinant).  The  corresponding  solution  for 

x  ■'■s  u+a  u-a 

x  =  J  (6  -  K)  l'(|)d6  +  f  (I  +  K)  !£U)dS 

where  K  and  the  lower  limits  of  the  integrations  are  arbitrary  constants. 
With  this  solution  (3)  becomes 

Pet  =  -4c 

which  must  not  vanish  if  the  transformation  is  to  be  valid. 

Returning  to  ( 7 ) ,  we  make  the  substitution 

2n  +  1 
7  "  2n  -  1 


Then,  for  integer  n  >  0,  the  most  general  solution  is 


(9) 
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where  ^  and  are  arbitrary  functions  of  their  arguments.  The  special 
cases  of  a  monatomic  gas  (7=5/3,  n=  2)  and  a  diatomic  gas  (7  =  7/3, 
n  =  3)  are  covered  by  this  solution.  The  fictitious  case,  7  =  3  (n  =  1) 
has  a  particularly  simple  solution 

t  =  §  +  a)  +  i2(u  -  n)3  (10) 

In  this  special  case,  c  =  0. 

The  solution  for  x  follows  from  (4)  and  (5),  either  directly  or 
through  the  following  transformation: 

y  a  x  -  ut 
Z  s  (It 

Then,  for  a  polytropic  gas, 

by  _  7  -  1  ciz 

do  ~  2  3u 

_  1  ~  1  z  I  7  -  3] 

ou  2  00  0  \  2  J 

For  7=3,  these  become  particularly  simple,  and  the  solution  is  (10) 
together  with 

y  =  -i^u  +  0)  +  *2(u  -  0) 
or 

x=  ut  -  i^u  +  0)  +  i2(u  -  0)  .  (11) 

A  useful  form  of  these  general  solutions,  which  follows  directly  from 
(lO)  and  (ll),  is 
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u  +  o  =  F[x  -  (u  +  a)t]  02) 

u  -  o  =  G[x  -  (u  -  o)t]  (13) 

Here  F  and  G  are  arbitrary  functions  of  their  arguments.  [It  should  be 
noted  that  (12)  and  (13)  can  be  derived  simply  and  directly  from  the  origi¬ 
nal  equations  as  a  direct  consequence  of  the  fact  that  for  7  =  3,  a  and  c 
are  identically  equal.*  A  directly  analogous  derivation  for  any  other  value 

of  7  is  not  possible.] 

Example:  Expansion  into  a,  Vacuum 

Initially  the  gas  is  at  rest;  for  x  <  0,  there  is  vacuum,  for  x  >  0 
there  is  gas.  Define 


0 

x  <  0 

H(x)  3- 

1 

V. 

x  >  0 

► 

W 

0  <  h(x)  <  1 

x  =  0 

J 

and  let  a  be  a  fixed  constant.  Then,  at  t  =  0, 


cr  =  a  H(x) 
u  =  a[H(x)  -  1] 


(15) 


The  form  of  u  follows  from  the  fact  that  the  motion  is  that  of  a  simple 
wave  in  which  u  -  cr  has  everywhere  the  same  value.  Combining  with  the 


*  For  7  =  3,  the  one-dimensional  Eulerian  equations  are 

— —  ®  „£(.±. >♦<»*.)£(«**-<>■ 

The  general  solution  of  this  pair  is  seen  to  be  u  ±  cr  =  F±[x  -  (u  ±  cr)t]  • 
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solution  (12)  and  (13),  we  get 

u  +  0  =  a  {2H[x  -  (u  +  a)t]  -  1 j 
u  -  0  =  -a 


or 


u  +  0>  * 
t 

U  +  0  <  £ 
"0 

u  +  0  =  £ 
t 


for  u  +  0  =  -a 
for  u  +  0  =  a 
for  -a  <  u  +  0  <  a 


which  seems  like  a  strange  way  to  express  the  solution,  but  this  can  be 
transformed  to  the  equivalent,  more  familiar  form 


«  -  I  <f  -  «) 

0  *  1  (f +  °> 


for  -a  <  ^  <  a 


for  the  region  of  the  rarefaction  wave.  In  particular,  at  x/t  =  -a, 
u  =  -a,  the  "escape  speed." 

It  is  also  strange  that  the  final  solution  depends  upon  the  nature 
of  the  velocity  profile  in  the  vacuum  (15).  Had  we  taken  u  =  0  at  t  =  0, 
the  final  result  would  have  been  meaningless. 

Considerable  additional  discussion  of  this  method  has  been  given  by 

t-s  w ,  *  _  , 

R.  von  Mises  ,  by  Courant  and  Friedrichs,  and  by  Landau  and  Lifshitz. 

R.  von  Mises,  Mathematical  Theory  of  Compressible  Fluid  Flow,"  Academic 
Press,  Inc.,  New  York,  1958. 

**Richard  Courant  and  K.  0.  Friedrichs,  "Supersonic  Flow  and  Shock  Waves," 
Interscience  Publishing  Company,  New  York,  1948,  pp.  88-91;  pp.  191-197. 

***L.  D.  Landau  and  E.  M.  Lifshitz,  "Fluid  Mechanics"  Pergamon  Press  Ltd., 
London;  Addison-Wesley  Publishing  Company,  Inc.,  Reading,  Mass.,  1959, 
pp.  586-592. 
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THEORY  OF  VISCOSITY 

General  Description 

Viscosity  in  fluids  may  arise  from  two  main  sources  which  are  quite 
different  hut  may  he  approximated-  hy  the  same  mathematical  form: 

1 .  Fluids  which  are  "thick"  and  "sticky"  derive  their  viscosity 
from  their  strong  cohesion,  or  intermolecular  attraction.  Shear  motion  is 
accomplished  hy  doing  work  to  break  these  attractive  bonds.  Re-forming 
of  the  bonds  leaves  the  molecules  vibrationally  excited,  equivalent  to 
conversion  of  energy  to  heat.  The  greater  the  heating,  the  more  easily 
the  bonds  are  broken,  so  that  such  viscosity  decreases  with  temperature. 

2.  In  gases,  the  intermolecular  forces  are  of  short  range  compared 
to  separation,  and  the  individual  molecules  move  with  much  greater  veloci¬ 
ties  than  the  local  average  of  many  of  them  (this  average  being  the  local 
fluid  velocity").  The  molecules  diffuse  about  among  each  other  (in  viola¬ 
tion  of  our  early  assumption  that  elements  remain  adjacent).  Thus,  for 
example,  a  boundary  between  helium  and  hydrogen  which  is  initially  sharp 
will  gradually  be  replaced  by  a  widening  zone  of  mixture.  Likewise,  the 
boundary  between  a  cold  region  (average  particle  speed  small)  and  a  hot 
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one  (average  particle  speed  large)  will  be  replaced  by  a  widening  warm 
zone  as  the  fast  moving  molecules  diffuse  into  the  slower  ones.  In  like 
manner,  if  there  is  a  slip  plane  in  the  gas  (that  is,  a  discontinuity 
in  mean  molecular  velocity),  then  the  effect  of  diffusion  will  produce 
an  apparent  drag.  This  is  because  the  molecules  moving  each  way  across 
the  slip  plane  carry,  on  the  average,  the  mean  velocity  of  the  side  of 
origin,  and  mix  this,  through  collisions,  with  the  mean  velocity  from 
the  other  side.  This  is  gaseous  viscosity.  Since  the  random  molecular 
velocities  (and  thus  the  diffusion  rates)  increase  with  temperature, 
the  viscosity  and  heat  conduction  will  increase  with  temperature  for 
this  type  of  process. 

The  two  viscosities  are  thus  considerably  different  in  their  de¬ 
pendence  on  temperature.  They  also  differ  in  that  heat  conduction  is 
directly  associated  with  the  second  type  but  arises  (if  even  signifi¬ 
cant)  from  quite  different  processes  in  the  first. 

The  usual  concept  of  viscosity  is  that  it  produces  a  drag  along  a 
shear  plane,  but  it  is  easily  seen  for  the  gaseous  type  how  viscosity 
can  affect  one-dimensional  motion.  If  there  is  a  velocity  gradiert 
parallel  to  the  velocity,  then  molecular  diffusion  can  tend  to  decrease 
the  gradient  in  just  the  same  way  as  when  the  gradient  is  perpendicular 
to  the  velocity.  The  fact  that  the  velocity  vector  can  diffuse  both 
parallel  and  perpendicular  to  itself  leads  us  to  suspect  that  a  precise 
mathematical  formulation  will  involve  the  use  of  tensors. 
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Formulation  from  the  Diffusion  Viewpoint 

To  accomplish  the  formulation  for  gaseous  viscosity,  we  consider 
first  the  momentum  equation  in  conservative  form.  In  true  one-dimensional 
motion,  the  equation  is  ( 1-24) 


The  term  in  parentheses  is  the  momentum  flux  past  a  point  fixed  in  space; 
the  first  part  is  the  transport  due  to  fluid  motion,  the  second  is  that 
produced  by  inter-element  forces.  If  there  is  a  gradient  in  fluid  veloc¬ 
ity  (i.e.,  in  mean  local  molecular  velocity),  then  we  expect  an  additional 
flux  term,  F,  which  depends  upon  the  size  of  the  gradient.  We  expand 
the  flux  in  a  Maclaurin  series: 

F  =  F  (^)  =  F(0)  +  ^  F'(0)  +  ••• 

The  flux  depends  also  upon  such  quantities  as  temperature  (and  so  there¬ 
fore  do  the  expansion  coefficients).  There  certainly  will  be  no  viscous 
flux  if  du/dx  =0;  also,  for  many  situations  only  the  lowest  contribut¬ 
ing  terms  in  the  expansion  are  required.  Thus,  we  write 


and  the  momentum  equation  becomes 


(D 


The  coefficient  of  viscosity  is  (iQ.  Note  that  the  viscous  term  can  be 
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considered  a  pressure  modification. 

In  three  dimensions,  the  momentum  equation  is  (l-5) 


+  V*(p  u*u}  +  Vp  =  0 

To  avoid  the  awkward  double  vector,  this  is  written  in  index  notation 

sr  *  <puiV  + =  0  (2) 

The  summation  convention  is  here  employed;  any  term  in  which  an  index 
appears  twice  is  to  be  summed  over  that  index: 


a 


3x7  (puiV  means  L  7ST  (puiuj} 


°  j=i  u 

The  three  components  are  Cartesian,  but  our  final  results  can  be  trans¬ 
formed  to  any  desired  coordinate  system.  It  is  convenient  to  introduce 
the  symbol  6^  defined  by  6^  =  0  for  i  £  J;  6  =  1  for  i  =  J.  Then 

(2).  can  be  put  into  the  form 


sr +  3^  Ieuiuj +  i>V  ’ 0  (3) 

The  bracket  term  is  again  a  flux  of  momentum.  The  first  part,  arising 
from  fluid  transport  past  the  reference  point,  signifies  the  rate  of 
transport  of  j -direction  momentum  in  the  i  direction.  Likewise,  the 
viscous  flux,  F^,  must  be  a  tensor  whose  components  denote  the  rate  of 
diffusion  of  j-component  momentum  in  direction  i. 

The  flux  tensor  must  be  symmetric.  This  can  be  demonstrated  by 
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the  following  simple  example.  Consider  a  small  cube  of  fluid  with  sides 

parallel  to  the  coordinate  axes.  Then  F^  is  the  x-direction  momentum 

per  unit  area  per  unit  time  which  enters  the  face  normal  to  the  y  axis 

closest  to  the  origin.  This  contributes  to  the  cube  a  rate  of  increase 

of  z-directed  angular  momentum  about  the  cube  center  of  (s  F^)(s/2) 

where  s  is  the  cube -side  length.  Through  the  other  face  normal  to  the 

y  axis  there  likewise  escapes  x-direction  momentum  at  the  rate  per  unit 

area  F  ,  and  this  likewise  adds  angular  momentum  in  the  same  direction 
xy 

as  that  through  the  other  face,  the  total  angular  momentum  per  unit  time 

3 

passing  through  the  y-direction  faces  thus  being  s  F  .  In  similar 

manner,  the  angular  momentum  rate  in  this  same  angular  direction  passing 

3 

through  the  x-direction  faces  is  -  s  r  .  Now  the  moment  of  inertia  of 
0  yx 

the  cube  is  l/6  (ps5)s2  so  that  the  rate  of  change  of  the  angular  ve- 

* 

locity  a)  is  given  by 


1  5  • 
£  p  S  CD 


Since  must  remain  finite  as  s  -*  0,  it  follows  that  F  =  F^.  The 
same  can  be  demonstrated  in  like  manner  for  the  other  component  pairs. 

The  viscous  flux  of  momentum  is  again  expected  to  vanish  when  there 
are  no  velocity  gradients;  we  keep  only  first  order  terms  in  the  expan¬ 
sion  in  powers  of  the  gradients.  Further  F^  must  depend  only  upon,  sym¬ 
metric  combinations  of  these  gradients.  These  symmetric  combinations 


(called  the  "rates  of  strain")  are 

du.  du 

eij  s  3^  +  5^ 


(4) 
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Besides  being  symmetric,  the  form  of  F. .  is  further  restricted  by 
the  requirement  that  the  equations  must  have  the  same  form  no  matter 
what  may  be  the  assignment  of  x,  y  and  z  to  i  and  j .  (More  precisely, 
the  viscous  flux  must  be  a  tensor.)  As  a  result,  F  must  be  related 

X  J 

to  e^  by  a  scalar  function  of  proportionality.  In  addition,  it  may  also  be 

related  to  5^  by  a  scalar  function  of  proportionality  provided  one  can  be 

found  which  vanishes  with  the  rates  of  strain.  Such  can  be  found  to  be  e, , 

kk 

and,  in  conventional  nomenclature,  we  finally  write  the  flux  tensor 


“  2  A  ekk  8ij 


ue 


i  J 


(5) 


Here  |i  and  A  are  scalar  functions  of  position,  termed,  respectively,  the 

* 

first  and  second  coefficients  of  viscosity.  The  generalized  momentum 
equation  has  thus  become 


p  i  d  1 

SC  +  3^  [puiu,)  +  (p  •  2*  ekk)  61J  -  -  0  <6> 

Notice  with  regard  to  scaling  of  this  equation  that  if  two  flows 
are  to  be  related  by  scaling,  then  X  and  u  must  be  proportionally 
changed,  and  the  "Reynolds  number"  scaling  factor (p^u^x^ )/mq  must  be 
unity  (i.e.,  the  Reynolds  number,  Rg  =  pux/u,  must  be  the  same).  These 
conditions  are  in  addition  to  those  of  Chapter  V. 

From  a  mathematical  point  of  view,  the  combined  effect  of  the 


*See,  however,  W.  H.  Liepmann  and  A.  Roshko,  "Elements  of  Gasdynami cs, " 
John  Wiley  and  Sons,  Inc.,  New  York,  p.  337  in  which  }X  +  2p  is  termed 
the  second  coefficient  of  viscosity.  The  reason  for  this  choice  is 
that  in  many  common  situations  3X  +  2u  s  0. 
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scalar  pressure  and  the  diffusion  flux  tensor  is  the  same  as  that  of  a 
generalized  stress  tensor 

8iJ  +  “eiJ  (?) 

which  is  interpreted  as  composed  of  the  components  of  normal  and  tangen¬ 
tial  force  per  unit  area  on  Cartesian  planes  within  the  fluid.  The  in¬ 
terpretation  is  clearly  artificial,  but  nevertheless  useful.  Most  dis¬ 
cussions  of  viscosity  use  this  stress-tensor  concept  almost  exclusively, 

obscuring  the  true  origin  of  viscosity. 

The  negative  average  of  the  three  normal  pressures  at  a  point  is 


called  the  total  pressure,  p^: 


-  "  3  pii 


-p-O^fiOsJ 


The  conservation  of  mass  equation,  which  is  not  affected  by  the  general¬ 
izations  of  this  chapter,  states  that  du./dx.  =  -  O/pJ  Up/dt),  so  that 
the  contribution  to  total  pressure  due  to  the  fluid  motion  is  propor¬ 
tional  to  the  compression  rate  of  the  fluid.  It  has  been  argued  that 
this  additional  pressure  is  fictitious  and  should  not  be  present,  so 
that  7  +  2/3  p  s  0.  Many  authors  have  treated  flow  problems  using  this 
assumption,  but  we  shall  keep  the  two  coefficients  independent  in  the 
discussions  that  follow.  The  second  term  in  (8)  is  called  the  "dilata¬ 
tion  pressure." 

Note  that  for  an  incompressible  fluid  (for  which  V*u  =  0)  second 
viscosity"  plays  no  role,  since  the  coefficient  of  *  vanishes  identically. 
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From  (6)  many  specialized  forms  may  be  derived  which  are  useful. 
Some  of  these  are  discussed  in  Chapter  DC,  in  examples  of  simple  viscous 
flows.  Equation  (6)  in  full  generality  can  be  transformed  to  vector 
notation 

P  «  Pg  -  Vp  +  V(AV* 3)  +  2(  V*pV)  ?+Vx(nVXu) 

where  g* is  the  body  acceleration.  Most  of  the  terms  can  be  verified  by 
inspection  of  components;  the  last  two  require  some  straightforward  re¬ 
arrangement  for  verification.  Transformation  to  various  curvilinear  co¬ 
ordinate  systems  follows  easily  from  the  vector  form.  If  A  and  p  are 
both  constants  (a  useful  approximation  for  many  problems),  then 

P  ^  =  Pg^  -  Vp  +  (A  +  p)  V  (V»u^  +  p(V.V)  n 


The  Energy  Equation 

It  is  instructive  to  examine  closely  the  derivation  of  the  energy 

equation  when  one  has  a  generalized  stress  tensor,  p  .  We  use  a  method 

>  ij 

which  is  much  like  that  of  Chapter  I,  but  carried  out  with  slightly  more 
* 

elegance.  Consider  a  closed  surface  moving  with  the  fluid 


Let 

W1  s  rate  at  which  work  is  being  done  on  the  fluid 
within  by  forces  at  the  surface 

Wg  =  rate  at  which  work  is  being  done  on  the  fluid 
within  by  body  forces 


*The  derivation  is  essentially  an  elaboration  of  that  presented  in 
L.  Howarth,  ed.,  "Modern  Developments  in  Fluid  Dynamics  —  High  Speed 
Flow,"  pp.  52-55,  Clarendon  Press,  .Oxford,  1953* 
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Q  =  rate  at  which  heat  is  conducted  across  the  surface 
into  the  fluid  within 

K  =  the  total  kinetic  energy  within 

dr  a  element  of  volume 

The  rest  of  the  symbols  are  as  used  before. 

Conservation  of  energy  then  requires 

i +  £  /p  1  dT  -  "i +  Ma + «  <*» 

where  the  integral  is  over  the  volume;  surface  integrals  will  be  over 
the  surface. 

Consider  now  an  element  of  area,  da,  on  the  surface.  Let  its 
outward  unit  normal  vector  be  n^.  Then  njP^  da  is  the  total  force  on 
the  area  in  the  i  direction.  This,  multiplied  by  the  velocity  compon¬ 
ent  in  the  i  direction,  gives  the  total  flux  in  that  direction,  and  the 
sum  of  these  gives  the  total  work  flux  across  the  element  of  area: 
n^p^u^da.  Thus 

wi  ‘/vuuida 

which,  by  the  divergence  theorem,  becomes 

Mi  ’  / 3^  <Vi)dT 

Tnia  may  be  transformed  as  follows: 

"l  *JU l 
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Now  this  equation  must  hold  for  any  surface;  therefore,  it  must  be  true 
that  the  integrands  form  the  equation 


DIIATATIVE  DISSIPATION 


Usually,  it  is  more  convenient  to  work  with  the  total  energy  per  unit 
mass,  E  =  I  +  Ju^  s  I  +  |u*u.  The  energy  and  momentum  equations  can 

be  combined  in  vector  form  to  give 

o  —  =  p  u.g  +  V‘[kVT  -  p  u*  +  Xutv*u *)  +  5  hV(u-u)  +  p(u*V)u)  05) 

K  dt 

Dilatative  Dissipation 

Any  element  of  a  compressible  fluid  may  be  expanding  or  contract¬ 
ing  as  it  moves.  Suppose,  in  particular,  that  an  element  is  being  com¬ 
pressed  the  same  from  all  sides  so  that  it  remains  geometrically  similar 
to  itself  in  shape.  If  we  transform  our  equations  to  a  uniformly  moving 
coordinate  system  which  instantaneously  has  the  average  velocity  of  the 
element,  then  the  velocity  in  the  vicinity  of  the  element  will  be  u*=  icr 
(where  k  is  a  constant)  to  lowest  order  in  l?|.  This  corresponds  to  an 
expansion  of  the  element  if  k  >  0  or  contraction  if  K  <  0.  (if  the 
velocity  to  lowest  order  were  u  =  ar  r,  then  this  would  correspond  to 
an  infinite  rate  of  change  of  density  for  n  <  0  and  no  rate  of  change 
of  density  for  n  >  0.)  Inserting  u*  =  icr  into  (l4)  reduces  it  to 

p  g  =  V-(k*T)  -  [p  -  3/C  (A  +  |  n)] 

The  quantity  in  brackets  is  the  total  pressure  defined  in  (8),  and  the 
significance  of  this  quantity  now  becomes  apparent.  Unless  the  dilata¬ 
tion  pressure  vanishes  (by  having  A  +  2/3  d  =  0),  then  there  is  a 
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viscous  dissipation  in  spherically  symmetric  dilatation. 

There  are  expansions  or  contractions  which  do  not  maintain  the 
geometric  similarity  of  an  element.  An  example  is  the  contraction  due 
to  a  plane  piston  pushing  against  a  fluid.  For  these,  vanishing  of  the 
dilatation  pressure  does  not  remove  the  dissipation. 

There  is  considerable  evidence  that  many  fluids  behave  as  though 
their  first  and  second  coefficients  of  viscosity  were  very  different > 
the  relation  between  them  being  quite  variable  from  fluid  to  fluid.  For 
such  fluids,  the  dilatation  pressure  is  a  real  and  important  quantity. 

Discussion  of  the  Dilatation  Pressure 

In  the  next  section  of  this  chapter  it  is  shown  how  the  total 
kinetic  energy  of  an  aggregate  of  monatomic  molecules  can  be  divided, 
in  macroscopic  viewpoints,  into  kinetic  and  internal  parts.  The  former 
results  from  the  mean  velocity  of  all  the  molecules,  while  the  latter, 
which  really  is  also  a  kinetic  energy,  results  from  the  individual  fluc¬ 
tuations  from  the  mean  velocity.  For  more  complicated  molecules,  there 
is  an  additional  contribution  to  the  internal  energy,  arising  from 
vibrations  and  rotations  within  the  molecules  themselves.  The  two  atoms 
of  a  diatomic  molecule,  for  example,  may  vibrate  along  the  line  between 
them,  or  may  rotate  about  either  of  two  axes  perpendicular  to  that  line. 
Thus  there  are  three  internal  "degrees  of  freedom"  in  addition  to  the 
three  translational  degrees  (corresponding  to  motion  in  each  direction 
of  three-dimensional  space) .  Only  at  high  temperatures,  however,  is 
the  vibrational  degree  of  freedom  excited. 
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Statistical  mechanical  arguments  show  that  in  equilibrium  of  the 
gas,  all  "active"  degrees  of  freedom  have,  on  the  average,  the  same  amount 
of  energy.  (Some  degrees,  such  as  the  diatomic  vibrational  one  mentioned 
above,  are  excluded  under  some  conditions  by  quantum  mechanical  effects.) 
If,  however,  some  influence  changes  the  state  of  the  gas  rapidly,  then 
there  may  be  a  delay  before  new  equilibrium  is  reached.  The  internal 
motions  may  be  excited  or  relaxed  only  after  numerous  collisions,  and 
so  lag  behind  the  changes  in  translational  internal  energy.  The  dilata¬ 
tion  pressure  arises  as  a  result  of  this  lag. 

Let  I.j  arid  Ig  be  the  specific  internal  energies  associated  with 
the  translational  and  internal  degrees  of  freedom,  respectively.  We  as¬ 
sume  that  the  translational  degrees  all  have  the  same  specific  energy, 
e1 ,  and  the  internal  degrees  all  have  the  same  specific  energy,  eg. 

Then  with  a  total  of  n  degrees  of  freedom  1^  =  Ig  =  (n  -  3)eg, 

the  total  specific  internal  energy  is  I  =  1^  +  Ig  -  +  (n  -  3)6g. 

It  is  also  shown  in  the  next  section  that  for  a  monatomic  gas  the 
pressure  is  given  by 

2  _ 

P  -  3  Pi 

In  that  special  case,  there  are  no  internal  degrees  of  freedom,  so 
that  I  arises  completely  from  translational  energy.  Likewise  for  more 
complicated  molecules,  the  pressure  is  still  associated  with  the  trans¬ 
lational  internal  energy  only,  so  that,  in  general. 
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-iI  +  t(n-3>  <ei  -e2> 

so  that 

p  =  —  pi  +  —  (n  -  3)p(e,  -  e-) 

(Note  that  in  equilibrium,  when  e1  =  e^,  this  equation  of  state  is  that 
of  a  polytropic  gas,  for  which  y  -  1  a  2/n,  or 

n  +  2 


The  validity  of  this  result  is  well  demonstrated  by  the  fact  that  under 
most  circumstances  such  diatomic  gases  as  oxygen  and  hydrogen,  which 
possess  five  active  degrees  of  freedom,  are  well  represented  by  y  =  7/5.) 

Next,  as  an  approximate  model  of  nonequilibrium  processes,  we  as¬ 
sume  that  the  energy  e^  is  always  tending  towards  e^  at  a  rate  which  is 
proportional  to  their  difference 


where  a  is  a  constant  with  dimensions  of  time,  the  "relaxation  time." 
Thus  the  pressure  is 

p  -  (7  -  1)pl  (n  -  5)p  5^= 


Finally,  we  assume  that  the  processes  are  nearly  in  equilibrium,  so  that 
a  is  small;  the  relaxation  time  is  short  compared  with  times  for  signifi¬ 
cant  change  in  the  gross  fluid  configuration.  Then,  since 
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1  2* 
n  Dt 


+  &  (a) 


we  get 

P  -  (7  -  1)pI  +  %  (n  -  3)p 
n 

p 

plus  terms  proportional  to  cc  which  are  here  neglected#  Finally,  we  use 
the  hydrodynamic  equation 


to  obtain 

p  =  (7  -  i)pi  - 

This  is  the  total  pressure,  as  defined  in  (8),  and  comparison  with  that 
equation  shows  that 

X  +  |  u  ■  — ( n  -  3) 

and  thus  vanishes  for  monatomic  gases  (n  =  3)  or  very  fast  relaxa¬ 
tion  (a  =  0) . 


L  n 


n  -  3) 


V*u 


Introduction  to  Rigorous  Viscosity  Theory 

The  diffusion  approach  to  derivation  of  the  viscous  equations  was 
filled  with  assumptions  which  may  he  accepted  as  reasonable,  but  which 
have  no  firm  basis  for  argument.  A  really  satisfactory  mathematical 
description  of  hydrodynamic  processes  can  arise  only  from  a  careful 
examination  of  the  detailed  molecular  dynamics.  We  here  present  a 
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highly  simplified  version  of  such  an  examination,  which,  however,  intro- 

* 

duces  many  of  the  features  of  the  more  elegant  treatments.  We  shall  see 
that  this  new  viewpoint  will  lead  to  some  surprising  reinterpretations 
of  previously  studied  hydrodynamic  processes. 

We  start  by  reviewing  the  concept  of  a  gas  (to  which  the  present 
discussion  is  limited) .  A  gas  consists  of  numerous  molecules  moving  very 
rapidly  with  infrequent  encounters.  We  give  more  precision  to  this  state¬ 
ment  as  follows: 

1 •  "Numerous"  means  that  within  any  region  of  macroscopic  space 
(i.e.,  any  region  over  which  the  mean  molecular  characteristics  do  not 
change  much),  the  number  of  molecules  is  large  compared  with  one. 

2.  "Rapidly"  means  that  the  velocity  of  an  individual  molecule  is 
large  compared  with  the  average  velocity  of  its  numerous  neighbors. 

3.  "infrequent"  means  that  a  molecule  spends  nearly  100$  of  its 
time  traveling  in  "free  space"  (i.e.,  not  in  interaction  with  any  other 
molecule) . 

Consider  the  following  idealized  situation.  At  t  =  0,  a  uniform 
distribution  of  molecules  is  moving  in  a  direction  normal  to  a  flat 
rigid  wall,  all  molecules  having  precisely  the  same  velocity.  Each 
molecule  hitting  the  wall  reflects  and  moves  back  into  the  approaching 


*  See,  for  example:  G.  N.  Patterson,  "Molecular  Flow  of  Gases,"  John 
.  Wiley  and  Sons,  Inc.,  New  York,  195^,  for  a  much  more  extensive 
introduction;  or  the  classical  treatment  in  detail:  Sydney  Chapman 
and  T.  G.  Cowling,  "The  Mathematical  Theory  of  Non-Uniform  Gases," 
Cambridge  University  Press,  1952. 
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cloud  with  the  negative  of  its  original 
velocity*  Thus  there  arises  a  widening 
zone  against  the  wall  defined  by  the 
motion  of  the  first  group  of  molecules 
to  reflect  from  the  wall.  Assume  at 
first  that  there  are  no  molecular  col¬ 
lisions.  Then  in  the  zone  next  to  the 
wall  there  are  just  as  many  molecules 
moving  to  the  left  as  to  the  right,  so 
that  the  mean  velocity  of  the  gas  is 
zero;  the  gas  has  "come  to  rest." 

The  leftward  moving  boundary  of  the 
zone  is  the  "front  of  a  shock." 

In  the  zone  there  is  no  macro¬ 
scopic  kinetic  energy,  because  the 
mean  gas  velocity  is  zero.  The  energy  of  the  molecules  is  now  called 
internal  energy,  and  it  is  easy  to  see  that  the  internal  energy  per 
unit  mass  is  just  equal  to  the  kinetic  energy  per  unit  mass  of  the  gas 
outside  of  the  shock  zone.  The  shock  speed  relative  to  the  approaching 
fluid  is  twice  the  speed  of  the  fluid  itself.  (Note  that  this  corres¬ 
ponds  to  y  -  3  in  our  previous  treatment  of  an  infinite  shock  in  a  poly¬ 
tropic  gas.  The  correspondence  is  not  accidental,  as  we  shall  see  below.) 

The  assumption  of  no  molecular  collisions  is  like  the  assumption 
of  no  viscosity  in  the  macroscopic  viewpoint.  The  effect  of  molecular 
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encounters  is  to  decrease  the  mean  leftward  speed  of  the  reflected  par¬ 
ticles,  particularly  of  the  earliest  ones  reflected,  thus  "smearing  out" 
the  otherwise- sharp  shock. 

The  individual  velocity  of  molecule  j,  u  ,  can  be  written  in  terms 

of  the  mean  velocity  in  its  neighborhood,  <  u*  >,  as 

J 

u .  =  <  u .  >  +  5u, 

J  J  J 


Thus  <  8u,  >  a  0.  Also,  since 
J 


— ►  2  -  -4  .2  2  _  ^  ^ 

u.  =  <  u.  >  +  5u.  +  2  <  u.  >*Su. 

J  i  0  <5 


we  obtain 


2<Uj  >  =  2<Uj>  +  2<  > 


06) 


That  is,  the  average  kinetic  energy  per  unit  mass  of  the  molecules 
(which  is  the  time  total  energy  per  unit  mass)  is  equal  to  the  kinetic 
energy  per  unit  mass  of  the  mean  motion  plus  an  additional  term  which 
can  only  be  identified  as  I,  the  internal  energy  per  unit  mass,  as  viewed 
macro scopi cally . 

The  fundamental  entity  in  gas  dynamics  is  the  distribution  func¬ 
tion,  f,  of  molecular  velocities.  It  is  defined  in  such  a  way  that 


f(u*,?,t)  d?  du* 


is  the  number  of  molecules  at  time  t  which  have  velocities  in  the  inter¬ 
val  du*  about  u*  and  which  lie  in  the  space  interval  d?  about  ru  We  shall 
suppose  throughout  that  f  ->0  as  |u*|  -»  ±  o»,  sufficiently  rapidly  that 
all  desired  integrals  converge.  Then,  with  M  =  molecule  mass. 
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00 

I  tKS) 


t)  du* s  p(?,t)  <  Q  > 


where  Q(u)  is  any  function  of  u  for  which  an  average  may  be  desired,  and 
p(r,t)  is  the  macroscopic  density,  mass  per  unit  volume. 

The  distribution  function  is  the  subject  of  Liouville's  Theorem, 
which  is  the  basic  starting  place  for  the  derivations  to  follow.  It  can 
be  stated:  Along  the  motion  of  any  particle,  f(u*,?,t)  is  a  constant; 


that  is 


f(u,r,t)  =  0 


To  demonstrate  the  validity  of  this  theorem,  we  consider  the 
integral  of  f  over  some  region  of  six  dimensional  coordinate-velocity 


space 


du  dr 


This  integral  represents  the  total  number  of  particles  within  the  region 
of  integration;  that  is,  the  total  number  of  particles  whose  coordinates 
lie  within  the  spatial  part  of  the  region  and  whose  velocities  lie  within 
the  velocity  part.  Now  suppose  that  the  region  of  integration  moves  in 
such  a  way  that  a  particle  on  its  surface  remains  always  on  its  surface. 
Then  no  particle  away  from  the  surface  can  ever  come  to  lie  thereon,  and 
the  number  of  interior  particles  is  forever  constant. 


—  JJf(u,r,t)  du  dr  =  0 

The  rule  for  such  differentiation,  (i-j),  can  be  generalized  to  this 
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six-dimensional  case;  with  the  summation  convention  for  the  integrand. 


du  dr  =  0 


or 


Jim 


du  dr  + 


JJ& 


‘du.  d(Du  /Dt)“| 
+  - 


-.  du. 

J  J 


du  dr  -  0 


Now  the  partial  derivatives  in  the  bracket  are  with  respect  to  the  inde¬ 
pendent  variables  x.,  u.,  and  t,  so  that  du./dx  =  0.  Also,  we  shall 

J  J  J  J 

assume  the  forces  on  the  particles  are  independent  of  velocity,  so  that 


Du. 

mBt= 


as  a  result  of  which  d(Du ./Dt)/du  =  0.  (The  assumption  of  velocity- 

J  J 

independent  forces  is  not  restrictive  for  our  present  purposes.)  Thus 

must  hold  for  any  arbitrary  region  of  integration  and  the  integrand  must 
vanish  identically. 

We  now  carry  out  the  differentiation  indicated  in  (l8)  and  write 
the  result  in  component  form 


or 


df  df  ^  df  ^ 

3t  +  5x7  Dt  +  du7  Dt 
J  J 


0 


df  df  df 

dt  +  U j  STT  +  m  SdT 
0  J 


0 


(19) 
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This  is  called  the  Liouville  Equation. 

As  a  simple  example  of  the  derivation  of  hydrodynamic  equations, 

we  consider  the  special  case  in  which  the  molecules  are  confined  to 
move  in  one  direction  only  (chosen  parallel  to  the  x-axis),  and  all 
quantities  of  interest  are  functions  of  that  coordinate  only.  Then 
(19)  becomes 


df  bf  F  df  _  n 

St  +  U^  +  m^'  0 


Multiply  this  by  Q(u)  and  integrate  over  all  values  of  u.  Then 

H s  *  ®  I 


du  =  0 


or 


00 


Qf  du  + 


Qu  f  du  + 


du 


0 


from  which 

dp  <Q>  dp<Qu>  _  pF  <^>  =  0 
St  5x  m  OU 


(20) 


If,  for  example,  Q  =  1,  then 

do  ,  dp<u>  _  0  (21) 

5t  5x  0 

which  is  the  familiar  conservation-of-mass  equation  with  a  rigorous  inter 

pretation  of  the  velocity  therein. 

The  other  two  conservation  equations  (1-24)  can  also  be  derived. 
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In  (20)  put  u  =  <  u>  +  6u.  Then 

dp<Q>  dp<Q>  <u>  dp<Q6u>  pF  ^  dQ 

~^r  +  — & —  +  &  i<5i 

Let  Q  =  u.  Then  this  becomes  the  momenttmi  equation 


> 


(22) 


O 

dp<u>  dp<u>  .  dp<u8u>_  pF 

+  it  +  ST"  *  m 

The  quantity  F/m  is  the  force  per  unit  mass  exerted  on  the  fluid.  Also, 
<u  8u>=  <(<u>  +  8u)8u>  m  <(8u)^>.  In  this  enforced  true  one -dimensional 
case,  there  can  be  no  viscosity.  A  molecular  encounter  can  only  result 
in  a  direct  back  scatter  which,  to  conserve  energy  and  momentum,  is  equi¬ 
valent  to  no  encounter  at  all.  This  means  that  we  must  identify  p<(8u)^> 
with  the  pressure  in  the  gas: 

p  =  p<(Su)2>  (23) 

1  2 

Next  we  let  Q  3=  —  u  ,  so  that 

<Q>  -  £  <(<u>+  8u)2>  =  ^  <u2>  +  i  <(&u)2> 

=  K  +  I 


the  sum  of  the  specific  kinetic  and  internal  energies  —  see  (16)  and 
the  discussion  following.  Thus  (22)  becomes 

+  4c  tP<u>(K  +  I)1  +  4c  [p<l"  u28u>]  = 

Now 

<  u28u>  =  <(<u>2  +  2<u>5u  +  5u2)8u>=  2<u>  <(5u)S»  +  <5u)^> 

=  2  <u>p  +  <  (  5u)^> 

P 
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so  that 

+  |  [p<u>  (X  ♦  I)]+  -  P  ^  -  |  (f  W>] 

The  terms  on  the  right  show  the  generalization  of  (l-24) .  The  first 
one  is  the  work  done  by  the  force  on  the  gas;  the  second  one  is  a  heat 

conduction  term. 

Note  that  the  identifications 

p  “  p  <(5u)2> 

1  =  1  <(6u)2> 

lead  to  the  equation  of  state 
P  -  2p  I 

which  is  the  polytropic  equation  of  state  with  7=3* 

Much  more  realistic  is  this  same  treatment  applied  to  a  three- 
dimensional  gas.  Multiply  (19)  by  and  integrate  over  all 

velocity  space.  Then,  in  the  same  way  as  above, 

F 

+  ar  [p<y <Q>  +  p<<i8uj>]  =  p<^,>  (2k) 

j  j 

With  ft  s  1 ,  we  get  the  mass  conservation  equation 


dp<u  > 


(25) 


With  ft  a  u±. 


we  get  the  momentum  conservation  equation 
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dp<ui>  d 

“St —  +  ST  [P<ui><uj>  +  P  <&ui6uj>]  =  — 

Thi3  is  to  "be  compared,  with  (6),  (in  which  equation  F^  s  o).  With 
Qs?u2,  so  that  <Q>  =  i<ui>2  +  i<(8ujL)2>=  K  +  I,  we  get  the 
energy  conservation  equation 

^%-"I)+  a|"  Ip  <uj>  (K  +  I)  +  p<uA>  <  SUjSu^] 

<3 

pFiui  a  o 

=  ‘  ST  ^  p  <(5ui}  5uj>] 

0 

The  generalized  stress  tensor,  corresponding  to  (7) ,  is  thus 


pu  ■  - p  <8ui  ‘V 


so  that  the  total  pressure,  analogous  to  (8),  is 


pT  =  \  p  <(5ui)2> 


1  2 

Because  I  =  ^  <  (6u^ )  >,  we  are  led  to  the  equation  of  state 


which  is  the  polytropic  equation  for  a  gas  with  y  =  5/ 3.  This  supports 
the  assertion  in  earlier  chapters  that  a  monatomic  gas  has  this  form  of 
state  equation. 

This  discussion  serves  only  as  an  introduction  to  the  rigorous 
study  of  gas  dynamics.  The  next  step,  which  we  shall  not  here  follow. 
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is  to  determine  the  distribution  function  for  situations  of  interest. 
With  it,  one  can  then  perform  the  averaging  integrations.  The  formula¬ 
tion  requires  the  use  of  additional  physical  principles,  and  leads  to 
the  famous  Boltzmann  equation.  The  subsequent  manipulations  are  long 
and  tedious,  but  generally  rewarding  in  their  results.  They  lead  to  the 
stress  tensor  (7)  as  the  lowest  order  approximation  to  reality,  and  pre¬ 
dict  the  dependence  of  A  and  p  on  the  details  of  molecular  structure. 
They  also  lead  to  higher  order  corrections  to  the  stress  tensor  —  in 
principle  the  results  can  be  as  accurate  as  desired.  The  reader  is 
directed  to  the  references  cited  earlier  in  this  chapter  for  the  de¬ 
tailed  discussion  and  for  additional  bibliography. 
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EXAMPLES  OF  VISCOUS  AMD  CONDUCTING  FLOWS 

The  Equations  Specialized 

To  illustrate  some  of  the  properties  of  viscous  and  conducting 

flows,  we  present  a  few  specialized  solutions.  Many  more  of  these  are 

* 

available  in  numerous  books  and  papers. 

We  specialize  to  plane,  two-dimensional  flow,  in  which  case  the 
appropriate  equations  are 


*See.  for  example,  S.  Pai, "Viscous  Flow  Theory,  I-Laminar  Flow,"  D.  Van 
No strand  Company,  Inc.,  Princeton,  1956. 
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Here  u  and  v  are  the  velocities  in  the  x  and  y  directions,  respectively, 

and  g  ,g  are  the  corresponding  accelerations  due  to  exterior  forces. 

^  12  2 
The  quantity,  k,  is  the  heat  conduction  coefficient.  E  =  I  +  ^(u  +  v  )• 

We  shall  limit  our  discussion  to  polytropic  gases  for  which 

p  =  (7  -  1 )p  I.  We  shall  further  simplify  the  equations  by  additional 

assumptions: 

1 .  The  two  coefficients  of  viscosity  are  proportional  to  each 
other:  A  =  Aj-i  where  A  is  a  dimensionless  constant*  (The  special  case 
in  which  A  +  2/5  p  3  0  is  included.)  In  general,  A  >  -  2/3* 

2.  The  simplest  model  of  a  gas  predicts  that  the  heat  conduction 
coefficient  and  viscosity  coefficient  are  related  to  the  specific  heat 


(here  a  constant)  in  such  a  way  that 


dl 


where  B  is  a  dimensionless  constant  related  to  y  by 


We  shall  assume  this  to  be  true  for  all  cases  discussed. 

The  equations  can  be  rearranged  and  condensed  considerably, 
especially  with  introduction  of  the  abbreviations 


Then 
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du  ci  chi\  eiQ 

p  dt  "  pgx  +  SE  (p  +  *  &  +  ^ 

dv  .  ^  dvv  dft 

P  ^  =  Pgy  +  ^  (P  +  2H  ^)  + 

P  f  =  p(ugx  +  Vgy} 

+  Jc  ^  +  ^  +  ^  Jc  ^BI  +  u2^ 

+  ^  [Pv  +  Qu  +  p  ^  (BI  +  V2)] 


(7) 

(8) 

(9) 


We  shall  not  try  to  work  with  these  equations  even  in  this  relatively 
simplified  form,  hut  shall  specialize  even  further.  Suppose,  for  example, 
that  all  quantities  are  independent  of  x,  and  g  =0.  Then 

y 


‘(Ie  +  T  If)  '  psx +  ^  %s>) 

4?  - 1)  =  £  [- p  +  “<A  +  2>  l' 

{- 


/8E  dE^ 

\5F  +  v  3^  ‘  pus*  +  ^  1‘  pv  +  M  ^ 


[ 


/  -  Ax  2  12 

( 1  +  g)  v  +  2  u 


+  BI 


(10) 

(11) 

(12) 


Diffusion  of  the  Viscous  Effect 

Consider  the  problem  of  determining  the  motion  of  a  semi -infinite 

gas  adjacent  to  a  flat  plate  at  y  =  0.  The  plate  commences  impulsively 

to  move  at  t  *  0  with  velocity  u  parallel  to  itself.  We  assume  that  u 

o  o 

is  small  enough  that  viscous  heating  of  the  gas  is  negligible.  Then  p 
and  p  can  be  considered  constants,  and  v  is  effectively  zero.  Under 
these  circumstances,  (10)  alone  is  sufficient;  it  becomes,  with  p/p  a  v, 
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du 

St 


(13) 


We  note,  from  arguments  of  dimensionality,  that  the  solution  must  be  a 
function  of  y/ -J  vt  only.  Substituting  this  into  ( 1 3)#  one  obtains  an 
easily  solved,  ordinary  differential  equation.  Application  of  the  boun¬ 
dary  conditions  yields 


u 


u=(' 


£ 


y/2/vt" 


-Z 


(1*0 


What  is  the  thickness  of  the  disturbed  area  after  the  plate  has 

moved  a  distance  L?  The  thickness  is  defined  as  the  distance  to  the 

point  at  which  u  =  1/2  uq.  At  that  point,  y^  =  0.95^  J  ^  V  vt". 

But  t  =  L/u  so  that 
o 


p 

For  air,  for  example,  v  ~  0.1 56  cm  /sec.  Thus,  if  the  plate  moves  at, 
say,  1000  cm/sec  (well  below  the  speed  of  sound  at  ordinary  conditions), 
then,  when  L  =  100  cm,  the  thickness  of  the  "boundary  layer"  is  only 
0.125  cm.  The  extreme  thinness  of  this  layer  for  objects  moving 
through  air  has  led  to  a  large,  highly  successful  field  of  approximate 
viscous  hydrodynamics  called  boundary- layer  theory.  Outside  of  this 
layer,  nonviscous  flow  theory  is  appropriate.  Within  it,  special  tech¬ 
niques  are  required. 

Problems  concerning  vehicles  or  projectiles  traveling  through 


-141- 


IX.  EXAMPLES  OF  VISCOUS  AND  CONDUCTING  FLOWS 


air  can  often  be  solved  without  requiring  the  use  of  general  viscous 
theory.  This  is  true  when  the  boundary  layer  is  always  in  thickness  a 
negligible  fraction  of  the  size  of  the  important  components  of  the  sys¬ 
tem.  A  vast  and  important  amount  of  literature  treats  many  interesting 
phases  of  these  boundary- layer  problems  —  stability,  resistance,  heat 
transfer,  etc.  We  shall  not  be  concerned  further  with  these  problems 
here. 


The  quantity  on  the  right  of  (15)  is  of  special  significance. 

If  L  is  the  length  of  a  projectile,  then  \j v/Luq  is  a  measure  of  the 
relative  influence  of  the  boundary  layer  on  the  external  flow  character¬ 
istics.  The  quantity  R  s  Lu  /v  is  called  the  Reynolds  number  for  the 
flow.  It  arises  frequently  in  problems  of  viscous  flow.  Such  quantities 
as  drag  and  lift  coefficients  are  often  expressed  as  functions  of  it. 


It  is  not,  perhaps,  inappropriate  to  notice  that  the  formula  for 

force  per  unit  areas  on  the  plate  F  =  u  (du/dy)^_Q  follows  easily  from 

(13).  The  momentum  per  unit  plate  area  of  the  fluid  above  the  plate  is 
00 

IQ  P  udy.  Thus,  from  ( 1 3),  its  time  rate  of  change  is 


This  is  the  force  on  the  fluid,  which  is  equal  and  opposite  to  that  on 
the  plate.  In  the  problem  at  hand. 
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Suppose  that,  in  this  nonsteady  problem,  the  plate  remained  fixed 
but  that  there  was  initially  some  velocity  distribution.  How  fast  does 
this  smooth  itself  out  due  to  viscosity?  It  is  still  appropriate  to  use 
(13).  We  expand  the  velocity  function  into  Fourier  integrals 

u(y,t)  =2 eikXe1CUt  dkdco 

00  -00 

Substitution  of  this  into  (13)  shows  that  it  is  a  solution  provided 

(io>  +  vk2)  =  0 


or 

V,.  \ 5(131  +  vk2) 

where  5(x)  is  the  Dirac  delta  function 

5(x)  m  0,  x  ^  0 
+b 

/  &(x)  dx  =  1 ,  a,b  >  0 

J  -a 

Thus,  the  general  solution  of  (13)  is: 

U(y,t)  .  ±  f\  elk*  e-^  dk  06) 

—  00 

This  shows  that  a  component  of  the  solution  with  wavelength  A  =  2x/k  de- 
cays  exponentially,  dropping  to  l/e  of  its  value  in  time  A  /4*  v.  In  air, 
for  example,  if  the  disturbance  is  of  length  100  cm,  its  decay  time  is 
~  400  sec. 
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Two  Steady-State  Solutions 

These  two  problems  concern  the  flow  between  parallel  plates. 

In  the  first  one,  there  is  no  external  force;  the  lower  plate  is  moving 
with  constant  velocity  uq.  In  the  second  one,  the  lower  plate  is  at 
rest  and  a  constant  acceleration  is  felt  by  the  gas,  parallel  to  the 
channel.  In  each  case,  the  v  component  of  velocity  is  zero,  so  that 
(10),  (11),  (12)  become 


d  ,  dux  - 
pg  +  ^  (H  ^p)  =  ° 

tr 0 

pug  +  [p  ^  (i  U2  +  BI)]  =0 


(17) 


In  problem  1,  g  =  0,  and  the  equations  can  be  integrated  [with 


p  =  (y  -  l)pl]  to  give 


P 


du 


a 


a 


(18) 


where  a,  £  and  a  are  constants. 

It  is  not  yet  necessary  to  specify  the  manner  in  which  p  varies. 

Let 


du  du  dl 

dy  ~  "5T  3y 


This  is  valid  provided  that  y  can  be  eliminated  from  the  equations. 
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We  substitute  this  into  (18)  and  eliminate  the  terms  with  y,  obtaining 

du  cr  du 
B  +  u  = 

for  which  n  is  also  eliminated.  This  can  be  integrated  to  give 

IB  +  \  u^  =  —  u  +  constant  09) 

As  an  example,  suppose  that  both  wall  temperatures  are  held  fixed  so 

that  I  =  I  at  top  and  bottom.  Then 
o 

p 

BI  +  i'  u  =  h  uu  +  BI 
*  ^  o  o 


At  mid-channel,  we  thus  have 

2 


I  = 


I 

o 


+ 


u 

o 


fcffi 


To  proceed  with  the  complete  solution  we  must  specify  the  variation  of 
p.  The  strongest  variation  of  p  in  most  gases  is  with  temperature;  we 
assume  that 


=  M- 


I 


n 


o 


where  p  and  n  are  constants.  From  U8) 
o 


2 

o 


From  this  may  be  obtained  the  first  integral 


n  dl  /  2  2  a2I 

*7  '7^ 


C2U) 
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p 

wherein  C  is  a  constant  of  integration.  A  second  integral  is  easily  ob¬ 
tained  for  certain  values  of  n  (such  as  0,  1);  for  others,  the  integra¬ 

tion  can  be  performed  numerically.  Qualitatively,  the  solutions  appear 
as  shown  in  the  figure,  in  which  the  ordinate  is  i/l^  or  u/uq. 

If  n  =  0,  then  u(y)  is  linear  and 
l(y)  is  quadratic. 

For  the  case  in  which  I  varies 

little  from  its  maximum  value,  I  ,  it 

m 

is  possible  to  obtain  an  approximate 
solution 

I  =  I _ 1  (y  -  y)^ 

x  m  2B  VSy/  Vym  y' 

0  y  top 

where  &u  is  the  total  change  in  velocity  over  the  channel  width,  &y,  and 

y  is  the  height  of  the  channel  center, 
m 

In  problem  2,  we  take  p  =  constant .  Then  ( 1 7)  becomes 


It  is  straightforward  to  obtain  the  lowest  order  solution  for  g  small: 
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where  I  =  I  at  top  and  bottom,  and  p  is  the  average  density.  Note  that 
o  o 


Im  =  Xo  +  3B 


so  that  the  central  material  speed  can  be  a  fairly  large  fraction  of  the 
initial  sound  speed  before  the  approximate  solution  is  no  longer  valid* 

We  may  proceed  to  obtain  an  approximation  to  the  build-up  of  this 
steady- state  solution.  We  begin  with 


du 


g  + 


(21) 


where  v  s  p/p.  This  equation  comes  from  (10)  in  which  we  put  v  =  0  and 

p  s  constant.  Using  the  result  of  (l6),  we  could  write  the  solution 

immediately  ^ 

.  1  P  ./,  \  ikx  -vk^t 

U  =  gt  +  2x  /  A'k'  e  e  ^ 

■■*00 

and  put  into  this  the  boundary  conditions  that  u  =  0  at  y  =  0,  y  =  2  y^. 

We  proceed  instead  as  follows.  Assume  that  there  exists  a  solu¬ 


tion 


Substitution  of  this  into  (21)  gives  an  equation  for  ¥ 


The  assumed  solution  is  clearly  wrong;  the  signal  from  the  edges  affects 
gradually  the  otherwise  uniformly  accelerating  gas.  Nevertheless,  the 
following  analysis  gives  a  good  representation  of  the  gross  gas  behavior. 
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Since  the  equation  for  ¥  must  not  depend  upon  y,  we  replace  the 
equation  by  its  average  over  the  channel 


1  2  d¥ 

5  ym  dt  =  6  ym 


or 


«  +  a, .  5g 

dt  2  1  y 
y  m 

•'m 


This  has  the  solution  (in  which  V  =  0  at  t  =  0) 


¥  = 


e  y. 


m 


1  -  expl 


'm 


This  leads  to  an  expression  for  u  which  can  be  integrated  over  the  chan¬ 
nel  to  give  the  time  dependence  of  the  momentum  per  unit  length 


2  pg  y. 


m 


3v 


1  -  expl 


'm 


Comparison  of  this  result  with  that  obtained  by  numerical  integration  of 
the  exact  equations  shows  surprisingly  good  agreement  for  the  approxi¬ 
mation. 


Finite  Thickness  of  Shocks 

The  inclusion  of  viscosity  and  heat- conduction  effects  in  the 
hydrodynamic  equations  can  remove  the  tendency  for  solutions  to  become 
discontinuous,  shocks  being  replaced  by  narrow  zones  over  which  the 
fluid  characteristics  vary  rapidly  but  continuously.  To  illustrate 
this,  we  examine  in  detail  the  structure  of  a  one-dimensional  com¬ 
pression  wave  moving  unchanged  through  a  polytropic  gas,  which  is  at 
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rest  at  x  =  +<».  The  appropriate  equations  are  derived  from  (7)  and  (9) > 
together  with  the  mass  equation: 


dp  .  dpu 

<5t  ox 


0 


du  ,  du  d 

p  St  +  pu  SE  =  55 


It  (I  +  \  U2)  +  PU  Ie  (I  +  \  U2)  =  4e  (-pu  + 


du 
du  t— 
3  dx 


where  p  a  s  A  +  2p  and  p  5  s  Bp,  and  po  is  the  initial  density.  We  look 
o  o  Q 

for  a  traveling  wave  solution,  in  which  all  quantities  are  functions  of 
x  -  vt  only.  The  constant  velocity,  v,  will  then  be  the  rate  of  propaga 
tion  of  the  wave.  Thus,  with  y  =  x  -  vt 


(-vp  +  pu)  =  0 

du  _  d 
dy  ~  dy 


(-vp  +  pu) 


(■ 


du 
P  o  j- 
o  dy 


\ 


(-vp  +  pu)  ^  (I  +  \  U2)  =  ±  (-pu  +  PodU  g  +  P. 


dy 


f  dl 
o^  dy 


These  may  be  integrated  once;  we  apply  the  boundary  conditions  that  at 
x  »  +«,  u  =  0,  p  =  pQ,  P  =  PQ,  I  =  IQ»  and  du/dy  =  dl/dy  =  0.  Then 


-vp  +  pu  -  -  VpQ 
du 

-vp  u=pcj  —  -  p+p 
Ho  o  dy  o 

,  12s  du  *  dl  *r 

-vpo(l  +  2  u  )  -  -  pu  +  PQPU  —  +  pc£^  -  vpoIo 

We  put  p  =  (y  -  1)  pi,  and  p  is  readily  eliminated  to  give 
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-vu 


(y  -  1)  Iv 

V  -  u 


+  (y  -  1) 


I 

o 


(22) 


v(l  -  I)  - 


V  u 


=  -  (y  -  i) 


uvl 

V  -  u 


„  du  ,  , 
cru  —  +  £ 
dy 


dl 

dy 


(23) 


Next  we  eliminate  I  (but  not  dl/ dy ) : 


„  dl  ,  v  -  u  „  du  _  u  f  2 
dy  7-1  dy  7  -  1  V  o 


2 


v 


+ 


1  +  1 
2 


(2!f) 


2 

where  c  B  y(y  -  1)  I  is  the  square  of  the  far-right  sound  speed.  Note 

2  2 

that  if  £  and  a  are  both  zero,  then  either  u  =  0  or  2( cq  -  v  )  +  ( y  +  1 )  uv 
=  0.  These  two  end  conditions  exactly  characterize  a  shock  [see  IV- 
(17)],  and  thus  justify  association  of  the  results  with  the  structure  of 
a  shock. 

We  shall  now  consider  two  extreme  cases  in  which  0  ss  0  or  £  =  0. 

Many  more  details  than  we  shall  present,  as  well  as  discussions  of  in- 

* 

between  cases,  have  been  given  in  articles  such  as  that  by  Hayes.  Even 
though  in  all  real  gases  the  effects  of  viscosity  and  heat  conduction 
are  of  nearly  equal  importance,  there  is  still  interest  in  these  two  ex¬ 
treme  cases.  With  viscosity  only,  the  results  are  applicable  to  analysis 
of  certain  finite-difference-approximation  methods  which  employ  "artifi¬ 
cial  viscosity"  to  improve  their  results.  With  conduction  only,  the 
analysis  refers  to  a  gas  in  which  radiation  is  transported.  The  true 
shock  structure  in  a  real  gas  does  not  differ  qualitatively  from  that 
in  the  viscosity-only  case.  The  results  of  the  conduction-only  case 

*W.  D.  Hayes,  in  "Fundamentals  of  Gas  Dynamics,"  H.  W.  Emmons,  ed., 

Vol.  Ill,  Section  D,  Princeton  University  Press,  1958* 
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show  that  a  discontinuity  may  still  be  present  even  if  £  is  very  large. 

Case  I,  Viscosity  Only.  With  £  a  0,  Eg..  (24)  is  easily  integra- 
ble  if  a  is  constant;  the  qualitative  features  can  be  seen  even  with 
variable  0.  The  equation  becomes 


du  u  ( 2 

0  —  =  - l  c 

dy  v  -  u  \  o 


Since  0  is  always  positive,  and  the  right  side  is  always  negative  and 
finite  for  u  between  its  extreme  values,  there  is  always  a  smooth  transi¬ 
tion  through  the  "shock  region."  The  velocity  profile  appears  as  shown 
in  the  figure 


The  width  of  the  shock  can  be  defined  as 

u  t 

w  -  du/dx|0 

where  subzero  refers  to  the  point  at  which  u  has  decreased  to  half  its 
value  at  x  =  Thus 
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and 

w  =  2a(yv2  +  <£) 

(7  +  l)v(v2  -  c^) 


(25) 


For  a  weak  shock  (v  «  c  )  the  shock  width  may  become  large.  For  a  shock 
of  infinite  strength  (c  =  0)  the  width  becomes 

-  207 

(7  +  1  )v 

2 

For  normal  air,  for  example,  a  «  0.1 5  cm  /sec  and  7  =  1.4,  so  that  for  a 
strong  shock 


/*+ 
sec.  Normal  air  sound  speed  is  about  3.3  X  10 

cm/sec;  the  speed  of  a  strong  shock  would  be  much  greater  so  that  w 
would  be  much  less  than  5  X  10~^  cm  —  a  very  short  distance  in  most  sit¬ 
uations  of  interest. 

Even  for  a  weak  shock,  the  thickness  may  be  quite  small.  In  the 
weak- shock  limit, 

v  w  _  >  2* _ (jL) 

(7  +  1)  cQ  \W 
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where  p  is  the  average  pressure  and  5p  is  the  change  in  pressure  across 
the  shock.  For  normal  air,  this  becomes 


w  *  1  X  10 


-5 


so  that  even  if  Bp  were  1 $  of  p,  the  shock  thickness  would  only  be  about 
0.001  cm. 

Case  II,  Conduction  Only.  With  o  s  0,  (22)  becomes  a  relation  be¬ 
tween  u  and  I  from  which  we  may  observe  directly  some  of  the  features  of 

* 

interest.  The  results  have  been  given  in  detail  by  Landau  and  Lifshitz. 
With  finite  heat  conductivity,  I  must  be  continuous  and  monotonically  in¬ 
creasing  towards  the  left.  Consider,  however,  the  root  of  (22)  which 
starts  at  the  far  right  with  u  =  0,  I  =  IQ.  As  u  increases,  so  also 
does  I,  and  if  the  shock  is  strong  enough,  I  will  reach  its  far-left 
value  before  u  does.  At  this  point  a  discontinuity  must  occur  in  which 
u  changes  to  the  other  root  of  (22)  for  which  I  equals  its  far-left 
value.  The  result  is  an  "isothermal  discontinuity"  in  velocity  and  den¬ 
sity.  The  velocity  profile  is  shown  qualitatively  in  the  following 
figure . 


*L.  D.  Landau  and  E.  M.  Lifshitz,  "Fluid  Mechanics,"  Pergamon  Press, 
Ltd.,  London;  Addison-Wesley  Publishing  Company,  Inc.,  Reading,  Mass.; 

1959,  p.  5^2. 
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Consider  a  gas  whose  state  is  perturbed  only  slightly  from  the  uni¬ 
form  state  in  which  p  =  pQ,  p  =  pQ,  I  =  Iq,  and  u  =  0.  The  equations  of 
the  previous  section  can  be  linearized  to  a  good  degree  of  approximation. 


Put  p  a  (7  -  1  Jpl  and  define  <p  a  p/p  ,  a  s  i/i  . 

o  o 

Then 


<*P  .  du  _ 
St  0 


lt  +  (7  -  >)  1 


d(q  +  qp' 
.  ox 
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£*o-oi-t§-o 


We  look  for  progressive  wave  solutions  of  the  form 


a  =  A  exp[i(kx  +  rnt)] 
cp  =  B  exp[i(kx  +  uit)] 
u  -  D  exp[i(kx  +  cot)] 


and  find  that  the  condition  for  existence  of  such  solutions  is 


0  ci) 

ik  ik 

(in)  +  £k2)  0 

or 

yok2  -  iCk4-  +  J — ^  ly--."  (ixii  +  dk2)(iCD  +  £k2)  =  0 

v  7  *"  /  o 

Consider  now  three  special  cases.  If  cr  =  £  =  0,  then 
O)2  «  7(7  -  1)  k2  IQ 

The  sound  speed  is  given  by 

Re  a) 
c  =  _ 

where  Re  means  the  real  part.  Thus  the  nonviscous,  nonconducting  limit 
gives  the  familiar  formula  of  previous  chapters. 

In  the  case  o  =  0,  £-><»,  the  sound  speed  (the  "isothermal  sound 


in)  +  ok 


17  -  U  Ic 
(7  -  l)ik 


=  0 
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speed")  is  thus  found  to  be 

a  —  °i  sot  henna  1  -  *J  ^  “  "*  ^o 

in  confirmation  of  the  speculation  of  the  last  section. 

If  a  -4 eo,  then  <o  is  purely  imaginary  and  the  sound  speed  is  zero. 
In  the  limit  that  o  and  £  are  small,  the  equation  for  co  can  be 
solved  for  the  lowest  order  viscosity  and  conduction  effects.  The  result 
is 

®  -  k  JA~7  -  i)io  +  [(7  -  D  £  +  yo] 

neglecting  quadratic  and  higher  terms  in  £  and  cf.  Thus  to  this  order  of 
approximation,  the  change  in  cu  due  to  viscosity  and  heat  conduction  is 
purely  imaginary;  the  sound  speed  is  not  affected.  The  correction  term 
results  in  a  damping  of  the  sound  waves  by  the  factor 

exp  {“  C(x  -  1)  S  +  7®]} 

For  a  sound  wave  in  normal  air  of  wave  length  100  cm  (so  that 
k  =  2x/l00  cm  ^ )  the  damping  time  will  thus  be  many  hundreds  of  seconds, 
during  which  time  the  signal  will  have  gone  many  miles.  The  viscous 
damping  is  in  that  case  negligible  compared  with  that  from  inverse-square 
attenuation. 

Neglect  of  the  higher-order  terms  (which  affect  sound  speed)  is 
justified  provided  that 

k  a  «  c 

o 
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where  is  the  unperturbed  sound  speed.  For  normal  air,  only  very 
short  wave  length  signals  are  affected,  with  negligible  change  noticed 
for  wave  lengths  significantly  greater  than  10  cm. 

Very  Large  Conductivity 

If  a  fluid  is  very  hot,  then  its  dynamics  may  be  affected  appre¬ 
ciably  by  the  presence  of  radiation.  At  sufficiently  high  temperatures, 
the  fluid  may  be  effectively  completely  transparent  and  in  continuous 
local  thermodynamic  equilibrium  with  the  radiation  field.  The  tempera¬ 
ture  is  then  the  same  everywhere  and  radiation  pressure  gradients  vanish. 
The  equations  of  mass  and  momentum  conservation  are  as  before,  but  the 
equation  of  state  becomes  a  relation  between  pressure  and  density  only. 
The  energy  equation  is  therefore  no  longer  needed  in  a  study  of  the 
fluid  dynamics  and  is  useful  only  for  determining  the  radiation  energy 
flux  necessary  to  maintain  the  isothermal  state.  We  may  thus  write 


||  ♦  v.(fiu)  =  o 

p  ||  +  p(u-V)u  +  a2  Vp  =  0 


where 


in  which  derivative  the  temperature  is  held  constant.  The  quantity  a, 
which  has  the  dimensions  of  velocity,  is  the  isothermal  sound  speed. 

It  is  a  function  of  density,  and  may  depend  explicitly  on  time,  through 


-157- 


IX.  EXAMPLES  OF  VISCOUS  AND  CONDUCTING  FLOWS 


variations  of  temperature  with  time. 

For  a  truly-one-dimensional  polytropic  gas, 


a  =J{ 

r 

- 1) 

I 

du 

Su 

bo 

St  + 

u 

SEE  + 

a 

bo 

So 

bu 

St  + 

u 

Sx  + 

a 

5x 

where  [see  III-(2)] 


a  s  a  In 


(28) 


(29) 


and  pQ  is  an  arbitrary  constant  with  the  dimension  of  density. 

Linearization  of  (28)  leads  to  the  wave  equation,  with  a  propaga¬ 
tion  speed  of  a  for  low-amplitude  signals,  thus  justifying  again  the 
equation  for  isothermal  sound  speed.  The  linearization  can  be  accom¬ 
plished  by  dropping  the  terms  with  u  undifferentiated,  or  by  the  follow¬ 
ing  procedure.  In  Lagrangian  coordinates,  (28)  can  be  written 


Su  °/a  bo 

St  +  a  e  SEE 

o 


=  0 


bo  a/a  Su 

St  +  a  e  5 £ 

o 


These  can  be  rearranged  to  give  the  following  exact  equation  for  o: 


a2  e~g/a 

St2 


2  82eff/a 


+  a 


=  0 


(30) 
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To  linearize,  expand  the  exponentials  in  power  series,  keeping  only  up 
through  the  first  order  terms;  the  result  is  again  the  simple  wave  equa¬ 
tion.  An  approximation  method  can  be  developed  based  on  retention  of 
successively  higher  terms  in  the  expansion.  Alternately,  with 
<p  =  e  3  p/  p  ,  the  exact  equation  can  be  written 


We  shall  now  demonstrate  a  few  simple  properties  of  the  isother¬ 
mal  polytropic  gas,  particularly  where  these  properties  differ  qualita¬ 
tively  from  those  of  the  corresponding  nonconducting  gas.  We  first  ob¬ 
serve  that  there  is  a  formal  correspondence  between  the  isothermal 
equations  and  the  isentropic  equations.  If  in  the  latter  we  let  7  -* 1 > 
then  p  is  proportional  to  p  and.  the  mass  and  momentum  equations  have  the 
same  essential  form  as  the  isothermal  equations.  Thus  it  is  possible  to 
derive  some  of  the  isothermal  relations  from  previous  results.  In  many 
cases,  however,  it  is  simpler  to  work  with  the  isothermal  equations  them' 

selves . 

Structure  of  the  isothermal  rarefaction  is  determined  from  the 
similarity  solution  in  which  u  and  a  depend  on  y  a  x/t  only.  The  equa¬ 
tions  become 


/  \  du  ^  «  da  -  n 

(u  -  v)  T““  +  a  U 

v  J'  dy  dy 


,  \  do  du  n 

(u-y)^-+a  —  =  0 

v  dy  dy 
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IX.  EXAMPLES  OF  VISCOUS  AND  CONDUCTING  FIflWS 


The  boundary  conditions  we  choose  are 

at  x  =  p  =  PQ  (so  that  d  =  0),  and  u  =  0 

at  x  =  +«>,  p  <  p  (so  that  a  <  0),  and  u  =  u  >  0 

J  O  00 

The  corresponding  solution  within  the  rarefaction  zone  is 


u  =  a  +  y 
a  =  -a  -  y 


Thus  the  density  profile  in  the  rarefaction  is 

(-1-x/at) 

e 

K  Ho 

At  x  =  -at,  p  =  p  .  No  matter  how  large  may  be  u  ,  the  density  will  not 
o 

drop  to  zero*  Thus  cavitation  cannot  occur  in  such  an  isothermal  rare¬ 
faction. 

The  shock  relations  can  be  derived  in  the  same  manner  as  in 
Chapter  IV.  For  a  steady  shock  moving  with  velocity  v  into  an  isothermal 
gas  at  rest,  the  following  relations  are  summarized 


(3D 


VERY  LARGE  CONDUCTIVITY 


Note  that  for  u  »  a,  then 


and  the  compression  across  the  shock  can  be  arbitrarily  high.  [This  is 
in  contrast  to  the  nonconducting  case  where  the  density  ratio  is  bounded 
by  the  quantity  [ y  +  1 )/( T’  -  1)]. 
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CHAPTER  X 


SOME  SPECIAL  PROBLEMS 

There  axe  certain  fluid-dynamic  problems  for  which  complete, 
exact  solutions  have  been  obtained.  Several  of  these  are  discussed  in 
this  chapter,  because  they  illustrate  well  some  of  the  methods  of  solu¬ 
tion  and  because  they  pertain  to  common  situations  of  interest.  They  are 
also  useful  for  comparison  with  the  results  of  approximation  methods,  for 
evaluation  of  their  accuracy. 

Problem  1 .  Flow  Past  a  Wedge 

Gas  flowing  past  an  infinite  two-dimensional  wedge  of  half  angle 
a  approaches  a  steady  state  configuration  as  time  passes.  If  the  incom¬ 
ing  gas  is  moving  sufficiently  fast, 
then  an  attached  shock  (dotted  line) 
is  formed.  Since  there  is  no  signi¬ 
ficant  length  to  the  system  (the 
appearance  of  the  configuration  is 
independent  of  magnification),  then 
the  shock  line  must  be  straight. 
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Furthermore,  the  trajectory  of  each  element  of  fluid  will  be  a  pair  of 
straight  lines  as  shown;  beyond  the  shock,  the  trajectory  must  be 
parallel  to  the  side  of  the  wedge.  In  the  diagram,  the  solid  line  is 
the  trajectory  of  a  fluid  element.  / 


It  is  convenient  to  resolve 
the  velocity  on  each  side  of  the 
shock  into  components  parallel  and 
perpendicular  to  the  shock.  Then 
the  relations  across  the  shock  can 


/ 


be  formed  in  just  the  same  way  as 

they  were  for  the  simple  one -dimensional  shock  in  Chapter  IV.  Let  m 
be  the  mass  per  unit  area  per  unit  time  crossing  the  shock. 


m  =  p  u  sin  0  =  pu  .gin  ( 6  -a) 
o  o 


(D 


To  conserve  tangential  momentum, 

mu  cos  0  -  mu  cos  ( 0  —  cat) 
o 

To  conserve  normal  momentum, 

p  +  mu  sin  0  =  p  +  mu  sin  (0  -  a) 
■*0  o 


(2) 

(3) 


To  conserve  energy, 

mE  +  p  u  sin  0  =  mE  +  pu  sin  (0  -  a)  W 

o  rO  o 

These  four  relations,  together  with  the  equation  of  state,  are  suffi¬ 
cient  to  determine  the  shock  angle,  as  well  as  the  field  variables  be¬ 
hind  the  shock,  all  in  terms  of  the  known  input  field  variables. 
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Consider,  for  example,  the  case  of  a  polytropic  gas  for  which 
2 

E=p/[(7-l)p]+l/2u.  In  terms  of  the  incoming  Mach  number 
(M  a  ratio  of  uq  to  initial  sound  speed),  one  can  show  by  straightfor¬ 
ward  manipulation  that 

tan  (e  -  a)  =  gig  (5) 

( 7  +  1 )  M  sin  6  cos  6 


For  given  values  of  M  and  a,  6  can  be  found,  and  the  other  field  varia¬ 
bles  obtained  from 


tan  6 

p  ~  po  tan(e  -  a) 


(6) 


cos  e 

o  cos(e  -  a) 


2  2 
p  =  p  +  p  u  sin  6  I  1 
o  o  o  1 


tan( 6  -  q) 
tan  6 


(7) 

(8) 


Note,  from  (5),  that  for  infinite  Mach  number  (incoming  gas  is  cold), 

T  -  1 

tan  (@  -  a)  =  - .  tan  6 

7+1 


so  that 


P  = 


7+1 

7  -  1 


in  agreement  with  (lV-25).  In  this  infinite  Mach  number  case,  one  can 
solve  explicitly  for  6 


tan  6  -  T7~ lltin-  ['  - J ’  -  (r2  -  1)  tan2  a] 


This  solution  is  real  only  if  tan  a 


<JVi 


7  -  1 ) .  For  greater  angles. 
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PROBLEM  2.  THE  SHOCK  TUBE 


the  shock  is  no  longer  attached  to  the  wedge,  and  a  different  procedure 

is  required  to  solve  the  problem. 

If  a  =  0  (there  is  just  a  point  disturbance  to  the  flow  field), 

then  from  (5)  sin  6  =  l/M. 

Convenient  graphs  and  tables  have  been  prepared  by  various 

* 

authors  to  show  numerical  results  for  many  situations. 


Problem  2.  The  Shock  Tube 

A  one-dimensional  hydrodynamic  system  is  divided  into  two  semi¬ 
infinite  sections  by  a  diaphragm. 


Initially,  there  is  gas  at  rest  on 
both  sides,  all  at  the  same  tem¬ 
perature.  To  the  left  of  the  dia¬ 
phragm,  the  gas  is  initially  at  a 
higher  density  and  pressure,  p_ 
and  p  ,  than  on  the  right,  p+  and 
p+.  On  both  sides  the  specific  in¬ 
ternal  energy  is  IQ. 

At  t  =  0,  the  diaphragm  is 
removed,  and  at  any  later  time  there 
is  observed  a  shock,  s,  moving  to 
the  right;  a  contact  discontinuity. 


*See,  for  references,  A.  H.  Shapiro,  "The  Dynamics  and  Thermodynamics  of 
Compressible  Fluid  Flow,"  Volume  1,  Art.  16.2,  The  Ronald  Press  Company, 
New  York,  1955* 


-165- 


X.  SOME  SPECIAL  PROBLEMS 


C ,  moving  to  the  right  and  a  rarefaction  wave,  hounded  by  points  a  and 
b,  moving  to  the  left. 

There  is  no  significant  length  to  the  system  (the  appearance  of 
the  configuration  at  a  later  time  is  a  magnification  of  an  earlier  ap¬ 
pearance),  so  that  each  of  the  points  moves  with  constant  speed.  The 
appearance  to  the  left  has  already  been  discussed  in  Chapter  III.  The 
contact  discontinuity  in  density  is  tentatively  allowed,  since  similar¬ 
ity  arguments  cannot  remove  it.  We  shall  see  that  if  the  density  is 
assumed  continuous  at  C,  then  the  problem  is  overdetermined.  The  point 
at  C  behaves  as  if  it  were  a  piston  pushing  with  uniform  speed,  and 
thus  produces  the  hind  of  shock  profile  already  discussed. 

There  are  eight  unknown  quantities  to  be  determined  in  the  sys¬ 
tem:  pL,  Uj,  pL,  IL  and  Pp,  pR,  IR.  Through  the  equation  of  state, 
two  of  these,  1^  and  IR,  can  be  eliminated.  Since  no  gas  passes  over 
the  contact  discontinuity,  we  must  have  u^  =  Up  (we  call  them  both  uj . 
Also  the  pressure  must  be  continuous  across  the  contact  discontinuity 
(otherwise,  there  would  be  an  infinite  acceleration);  we  put  pT  =  pD  =  p 

b  X\  £ 

We  are  thus  left  with  four  unknown  quantities,  p^,  Pp,  u^ ,  p^  ,  for 
which  four  relations  are  needed.  These  are  obtained  as  follows  (we  as¬ 
sume  a  polytropic  gas):  Across  the  shock,  we  use  the  relations  of 
Chapter  IV  in  the  forms 


(p+-pc 


2 

C 


(9) 


and 
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p+  ~  %  p+  ~  % 
p+  +  Pe  "  7  P+  +  Pr 

Across  the  rarefaction,  entropy  is  conserved  —  see  (II-4), 


(10) 


01) 


Finally,  a  characteristic  line  dx/dt  =  u  +  c,  goes  across  the  rarefaction 
—  see  Chapter  III  —  so  that 


(12) 


We  thus  have  the  required  four  relations  among  the  four  unknowns.  It  is 
convenient  to  define 


(13) 


so  that  X  is  known  from  the  initial  configuration.  Straightforward  mani¬ 
pulation  of  (9)  through  (12)  results  in  an  equation  determining  P  in 
terms  of  X  and  7: 


(1  -  P) 

1  +  P)  -1  +  P 


(H) 


Some  values  of  P  for  various  values  of  X  and  y  are  shown  in  the  following 
table . 
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X*4  7 

i-h  \ 

4/2 

5/5 

6/4 

7/5 

8/6 

3 

1.66859 

1.68018 

1 .68789 

1 .69339 

1 .69751 

4 

1.88474 

1.90514 

1 .91875 

1 .92846 

1 .93574 

5 

2.06428 

2.09391 

2.11371 

2.12787 

2.13849 

8 

2.47558 

2.55295 

2.57150 

2.59917 

2.61997 

16 

5.15841 

3.28224 

3.38650 

3.42740 

3.47343 

With  P  determined,  the  rest  of  the  unknown  quantities  are  easily- 
calculated: 


The  shock  moves  with  absolute  speed 


Uc  PR 
^  "  p+ 


(20) 


The  point  c  moves,  of  course,  with  absolute  speed  u^  ,  while  the  points 
a  and  bmove  with  sound  speed  relative  to  the  gas,  or  with  absolute  speeds 
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PROBLEM  J.  SHOCK  HITTING  A  DENSITY  DISCONTINUITY 

v  "c  -y^7  ■ 15  jl  (2'} 

v  -  y^K  (22) 

Problem  3.  Shock  Hitting  a  Density  Discontinuity 

In  the  simplest  form  of  this  problem,  a  one-dimensional,  steady- 
state  shock  passes  through  a  uniform  polytropic  gas  which  is  cold  and  at 
rest.  At  some  point,  it  strikes  a  discontinuity  in  material,  beyond 
which  there  is  another  (different)  uniform  polytropic  gas  which  also 
initially  is  cold  and  at  rest.  We  shall  tentatively  assume  that  a  rare¬ 
faction  is  reflected  back.  The  condition  under  which  this  is  the  case 
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will  follow  from  the  analysis.  In  any  case  a  shock  will  he  transmitted 
into  the  material  to  the  right. 

The  incoming  shock  is  characterized  by  the  infinite  shock  rela¬ 
tions  (IV-2-5): 


The  analysis  follows  closely  that  for  the  shock  tube  problem.  The  same 
quantities  (with  the  same  symbols)  are  unknown,  and  the  four  equations 
determining  them  are  [arising  as  in  (9)  to  (12)]: 
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Define,  for  convenience 


A  ■ 


0-0  Pi  +  'J 

o  ^2  +  1  ) 


(25) 

(26) 


where  A  is  known  from  the  input  conditions  and  P  is  to  he  determined. 


Some  manipulation  shows  that  P  is  determined  in  terms  of  A  and  7-j  hy 


2A 

7i  -  1 


V  -r  VP  -  ( V  P  J  +  y  -~T 

?i  -  1  L  v  /  J  V  7i 


(27) 


If,  now,  a  rarefaction  is  indeed  reflected,  then  P  >  1  is  required. 
This  corresponds  to  A  >  1 ,  which,  therefore,  is  the  required  condition 
for  a  reflected  rarefaction.  Some  values  of  P  for  various  values  of  A 
and  71  are  given  helow. 


4/2 

5/3 

6/4 

Hi _ 

8/6 

2 

1 .44587 

1 .47088 

1 .49165 

1.50933 

1 .52477 

4 

2.18104 

2.24690 

2.30289 

2.35154 

2.39453 

8 

3.42041 

3.55215 

3.66725 

3.76930 

3.86089 

12 

4.52404 

4.71103 

4.87742 

5.02684 

5.16226 

With  P  known,  the  other  unknown  quantities  in  the  system  can 
he  calculated? 
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The  points  o,  and  Jb  move  to  the  left  with  sound  speed  relative  to  the 

material.  One  may  notice,  however,  that  the  point  jb  moves  to  the 

right  relative  to  the  rest  frame  if  7^  <2. 

Notice  that  if  p  =0  (.the  shock  hits  a  free  surface),  then  one 
+>° 

can  derive  the  fact  that  the  free  surface  will  move  with  the  sum  of  the 
material  speed  in  the  shock  plus  the  escape  speed  of  the  shocked  mater¬ 
ial  (as  derived  in  Chapter  III).  This  is  proved  as  follows.  As 
p+  q  -»0,  A  -»  oo  and  P  -» <»,  From  (27),  we  see  that 


so  that,  from  (28), 


With  (IV-25),  this  becomes 


2c. 

UA  -»  U  +  - =—zr 

C  -  71  -  1 

with  the  second  term  being  just  the  escape  speed  ( III-7 ) . 
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A  second  case  must  also  be  considered,  "that  in  which  a  shock  is  re¬ 
flected  back  from  the  discontinuity.  It  is  expected  that  this  will  occur 
if  A  <  1 .  We  shall  indeed  see  the  general  result 


Since  P  <  1  is  required  for  the  shock  reflection,  we  see  that  this  means 
A  <  1,  completing  the  proof  of  (29).  Some  values  of  P  for  various  values 
of  A  and  7  are  given  below. 
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4/2 

5/5 

6/4 

7/5 

8/6 

0.2 

0.505768 

0.480057 

0.460067 

0.443920 

0.430509 

0.4 

0.658085 

0.639622 

0.625000 

0.613003 

0.602905 

0.6 

0.704182 

0.772321 

0.762822 

0.754954 

0.748275 

0.8 

0.896522 

0.890783 

0.886148 

0.882281 

0.878982 

With  P  known,  the  other  quantities  in  the  system  can  he  calculated 


"7,(1  +  P)  +  (1  -  P)~ 
PL  =  P-  17,(1  +  P)  -  (l  -  P). 


PR  =  p+,o 


,  %  (P_u_)  -  (p2uc) 

v  (absolute)  =  - UL- 

SL  p  -  PL 


> 


J 


(35) 


In  the  special  case  that  P+,Q  -  00  (the  shock  strikes  a  rigid  wall  and  re 
fleets  back),  then  A  -  0  and  P  -  (7]  -  1)^  .  In  this  case. 


v  (absolute)®  -(7.  -  1)  u 

ST  I  * 


(36) 
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Another  question  of  interest  is:  How  can  one  use  such  a  reflective  sys¬ 
tem  to  produce  a  stationary  shock  for  study?  This  can  be  accomplished 
if  the  reflected  shock  at  s^  moves  with  zero  speed.  According  to  (35), 
the  condition  for  this  is  p  u  =  p  u  .  The  result  of  some  algebraic 
manipulation  shows  that  this  can  be  realized  provided 


A  ■ 


(r-,  +  i) 
Ttprry* 


(7,  -  D5 

5-  h 


(37) 
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Problem  4.  The  Reactive  Shock  or  Detonation 

We  assume  that  the  conditions  on  each  side  of  the  detonation  front 
are  uniform,  and  that  the  detonation  itself  takes  place  instantaneously 
along  a  moving  discontinuity. 


Actually,  the  total  internal  energy  per  unit  mass  on  the  right 
side  is  not  zero,  but  K,  the  chemical  energy  per  unit  mass.  With  this 
change,  the  shock  relations  [Equations  (l4),  (15)>  06)  of  Chapter  IV] 
are  directly  applicable 


It  is  assumed  that  u  is  known;  it  is  the  velocity  of  the  driving  pis- 

s 

ton.  We  assume  further  that  the  detonation  products  can  be  represented 
by  a  polytropic  state  equation;  experiments  show  that  this  is  a  reason¬ 
able  assumption  for  many  explosives . 

These  equations  can  be  solved  as  follows  [we  have  put 
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Ps  =  (7  -  1)  PSIS1: 


"(2K/u^)  +  (7  +  1 )jy  -  1) 
2K/u  +  1  . 


(38) 


Is  =  K  +  iu^  (39) 

v  -  (7  -  1)  ug  (^0) 

s 

In  the  limit  that  K  -»  0,  these  reduce  to  the  familiar  relations  for  an 
infinite  shock.  With  K  >  0,  the  detonation  front  moves  faster  than  the 
infinite  shock,  and  the  density  "behind  the  front  is  correspondingly  less. 

We  note  that  the  results  predict  that  as  u  -*  0  ( for  K  fixed) 
the  detonation  speed  "becomes  infinite.  We  thus  expect  that  the  results 
are  in  error  for  small  piston  speeds.  Indeed,  there  is  another  reason 
for  believing  this.  We  note  that  the  sound  speed  just  behind  the  detona¬ 
tion  is 

Cg  =  J  7(7  -  1)(K  +  i  uj) 

For  sufficiently  high  piston  velocities,  u  +  c  is  thus  greater  than  vj 
that  is,  signals  from  the  piston  can  overtake  the  detonation  front  from 
behind  and  thus  influence  it.  For  small  piston  velocities  the  above 
model  would  predict  u  +  c  <  v  and  no  signal  could  overtake  the  front 
from  behind.  Thus,  for  piston  speeds  less  than  v  -  cg,  no  signal  from 
behind  the  front  can  catch  up  with  it,  and  the  detonation  process  must 
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proceed  in  a  manner  independent  of  the  motion  of  the  piston.  We  thus 

conclude  that  for  an  underdriven  detonation,  the  front  must  move  exactly 

as  for  a  critically  driven  detonation;  that  is,  one  for  which  c  +  u  =  v. 

s  s 

(This  condition,  which  determines*  the  motion  of  an  under-driven  detona- 
tion,  is  called  the  Chapman -Jouguet  condition.)  Substitution  of  the  con¬ 
ditions  for  an  over-driven  detonation  into  the  Chapman-Jouguet  relation 
allows  derivation  of  the  following 


(41) 

(42) 


Ps  -  2(r  -  1)  P„  K 


I  --52 


7  +  1 


7  +  1 

P  =  1 -  P 

s  7  o 


m 

m 


Just  as  for  an  infinite  shock,  the  compression  depends  on  7  only. 

Determination  of  the  profile  behind  an  underdriven  detonation 
can  be  made  by  means  of  the  similarity  method  of  Chapter  V.  We  outline 
the  procedure  briefly.  The  equations 

du  .  du  dp 

p3t  +  '  S 

P  =  Apr 


are  applicable  wherein  A  is  a  constant  determined  by  the  value  of  K. 
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The  substitution  y  =  x/t  and  the  assumption  of  dependence  upon  y  only 


lead  to 


„3£  +  te=  O 

J  dy  dy 


,  N  du  dp 

p(u  -  y)  55  -  -  aj 

A  7 

p  =  Ap 

If  these  are  to  represent  a  detonation,  then  it  is  necessary  that  the 
boundary  conditions  be  satisfied: 


It  thus  is  necessary  that  the  value  of  A  (which  determines  the  amount 
of  entropy  behind  the  detonation)  be  given  by 

A  -  2 (r  -  i)  (^rfr)  o\'y  k 


The  solution  can  be  carried  through,  and  the  boundary  conditions  indeed 
satisfied;  the  result  is 
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1  (ox.  ^ 
U  y  +  1  \t  "  V 


p  =  21 


+  1  * 
7  Po- 


T  7  -  1 )  r  +  ▼* 


TV 


2 
7- 1 


The  solution  continues  back  to  the  point  where  u  =  u^,  the  piston  ve¬ 
locity.  If  p  =  0  at  the  piston  (so  that  the  detonation  takes  place 
against  a  vacuum)  then  the  free  surface  of  the  gas  moves  such  that 

v 


X 

t 


1 


If  the  detonation  takes  place  against  a  rigid  wall,  so  that 

^  =  0,  then  the  velocity  profile  at  some  time  after  t  =  0  is 

V 


the  similarity  solution  that  a  =  \  vt.  What  is  the  total  displacement 
of  material  caused  by  the  passing  over  of  such  a  detonation  wave?  This 
is  determined  by  solving  the  differential  equation: 


dx 

dt 


u(x,t) 


r2x 


+  1 


v) 


=  0 


for  —  <  x  <  vt 
otherwise 
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This  has  the  solution 

2 


vt 

7  -  1 


A  graph  of  this  result  is  shown  as  follows: 


A  nonviscous,  nonconducting,  semi-infinite  polytropic  gas  flowing 


parallel  to  a  flat  surface  approaches  a  comer  in  the  surface.  If  the 
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surface  "bends  into  the  gas  flow,  then  the  wedge  theory  (Problem  1 )  may 
be  applicable.  If  the  surface  bends  away,  then  a  steady  expanding  flow 
may  occur,  (in  both  cases,  the  gas  speed  must  be  at  least  that  of  sound 
for  a  steady  state  to  be  present.)  The  latter  flow,  which  concerns  us 
here,  is  often  called  Prandtl-Meyer  flow. 


We  choose  cylindrical  coordinates  with  origin  at  the  corner;  u 
and  v  are  the  components  of  velocity  in  the  r  and  6  directions,  respec¬ 
tively,  and  the  equations  are 

du  (  du  v  duN  pv2  2  dp  _ 

p  dt  +  <\u  d?  +  ?  da/  "  V  +  c  d?  “  0 

dv  f  dv  v  dv\  puv  c2  dp  _ 

p  dt  +  p\  dr  +  r  dp/  +  r  +  r  df? 

c2  ■  y  Ap^-1 

where  A  is  a  fixed  constsmt  related  to  the  entropy  of  the  gas  —  see 
(II-4).  In  steady-state  flow,  the  time  derivatives  vanish.  Also,  since 
there  is  no  length  scale  to  the  system,  its  appearance  must  be 
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independent  of  magnification;  that  is,  u,  v  and  p  must  be  independent  of  r. 
With  these  conditions,  the  equations  simplify  considerably: 


From  the  first  and  third  of  these, 

(v2-=2)i- 

Either  v  s  c,  or  p  is  identically  constant.  The  latter  alternative  leads 
to  the  trivial  solution  of  flow  with  no  corner  at  all,  parallel  to  the 
initial  boundary.  It  is  applicable  up  to  the  point  at  which  v  =  c.  In 
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or 


=  COS 


(47) 


This  is  the  angle  of  the  "Mach  line"  from  the  corner.  Beyond  that  angle, 
the  gas  begins  to  turn. 

With  v  3  c,  the  equations  become 


pu 


de 

du 

de 


=  0 


Now 


3  b2p  ^ 

de  p  de 


where  we  have  introduced  the  abbreviation 


b2 


Thus 


,2  dc 

u  +  b  —  =  0 

de 

du  ~ 

C  -  da’  0 


These  have  the  solution 


0  B 

u  =  bA  sin  —  +  bB  cos  — 

b  b 

6  Q 

c  =  A  cos  r-  -  B  sin  — 

b  b 


m 


where  A  and  B  are  constants  of  integration.  The  general  boundary  condi¬ 
tions  are 


At  e  =  v 


V  =  C  =  co. 


u  =  vo  sin  ^ 


(49) 
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We  look  first,  however,  at  the  case  0,  =  0  (i.e.,  v  =c).  For  this 
’  *  M  ’  o  o 

simpler  case,  the  solution  becomes 


u 

v 


b  c 

o 


sin 


0 

b 


c 


0 

c  cos  r- 
o  b 


(50) 


How  does  the  radial  distance  to  a  stream  line  vary  with  angle?  From 


dr 

dt 


d0 

dt 


u 

v 

r 


we  derive 


dr  ur 
d0  _  v 


For  the  solution  (50),  this  can  be  integrated  to  give 


(51) 


Note  that  the  sound  speed,  and  thus  the  density,  drops  to  zero  at 
0  =  bx/2;  the  gas  will  turn  through  no  greater  angle  than  that.  One  can 
show  easily  that  the  angular  deflection  of  the  stream  line  from  horizon¬ 
tal,  a,  is  given  by 


a 


0  -  ^  +  tan 


(52) 


For  the  more  general  boundary  conditions,  (49),  analogous  formulas  may 
be  derived.  Alternately,  the  supersonic  input  flow  can  be  assumed  to 
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have  turned  from  sonic,  through  some  angle  6q,  and  subsequent  turning 
measured  from  the  radius  vector  at  -0q.  To  accomplish  this  it  is  useful 
to  express  6  as  a  function  of  local  Mach  number  M^  which  is  defined  as 
the  ratio  of  the  gas  speed  to  the  sound  speed: 


Combination  of  this  with  (50)  leads  to 


6  -  b  tan 


Thus,  if  the  incoming  Mach  number  is  Mq,  then 


6 


o 


■  b 


(53) 


In  turn,  one  can  find  the  fictitious  sound  speed  from  which  this  turn 
started,  and  the  fictitious  initial  radius  for  any  stream  line,  and  con¬ 
tinue  the  solution  as  if  from  the  initial  sonic  flow.  This  alternative 
procedure  is  especially  useful  for  situations  for  which  tables  have  been 
provided  for  Prandtl-Meyer  solutions  turning  from  sonic. 

Finally,  for  reference,  we  include  the  formula  for  a  as  a  function 
of  M,  derived  from  (52)  and  (53): 


(5^) 
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Problem  6.  Instability  of  Interface  Between  Fluids 

There  are  various  examples  of  fluid  motions  in  which  the  question 
of  interface  stability  arises.  If,  for  example,  a  heavy  fluid  is  sus¬ 
pended  over  a  lighter  one  in  a  downward  gravitational  field,  any  irregu¬ 
larity  of  the  interface  will  increase  in  amplitude.  The  upper  fluid 
will  fall  into  the  lower  one  in  a  set  of  narrow  penetrating  spikes, 
while  the  lower  fluid  will  float  up  in  round-topped  bubbles.  This  is 
an  example  of  "Taylor  instability." 

Another  type,  known  as  "Helmholtz  instability"  occurs  along  a 
slip  plane  between  two  fluids  (or  within  one  fluid).  Any  slight  irregu¬ 
larity  is  amplified,  resulting  in  mixing  if  there  is  no  counter  process. 
An  example  is  seen  in  the  formation  of  water  waves  from  the  wind;  another 
is  in  the  flapping  of  a  flag. 

These  two  types  of  instability  are  best  known  for  incompressible 
fluids  (that  is,  for  fluids  whose  motions  proceed  at  velocities  very 
small  compared  with  their  sound  speeds).  When  accelerations  are  great 
and  Velocities  large,  then  effects  of  compressibility  can  become  impor¬ 
tant.  Taylor  instability  problems  then  refer  to  the  effects  of  a  com¬ 
pression  wave,  or  shock,  sweeping  across  an  irregular  interface. 

We  shall  here  discuss  several  types  of  surface  instability,  in¬ 
cluding  mainly  those  for  incompressible  fluids  because  of  their  ease  in 
solving  and  because  of  their  qualitative  applicability  to  some  compresr 
sible-fluid  situations. 
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Taylor  Instability,  Fluid  Incompressible 

We  consider  the  case  of  low  amplitude  interface  motions.  For  an 
incompressible  fluid,  in  which  the  density  of  an  element  remains  forever 
constant,  the  mass  equation  in  two-dimensional  flow  becomes  —  see  I-(2) 


du  dv 


o 


(55) 


in  which  u  and  v  are  the  velocity  components  in  the  x  and  y  directions* 
respectively.  In  addition,  the  momentum  equation,  with  vertical  acceler¬ 
ation  only  (positive  upwards)  breaks  into  the  two  components  —  see  I-(33). 


8u  dp 

p3t  +  ^  = 


o 


(56) 


dv 

p3t 


+ 


pg 


0 


(57) 


in  which  we  have  dropped  the  transport  term  (u *V)u*  because  of  the  smallness 
Of  the  velocities.  (The  various  approximations  employed  in  this  section 
can  be  verified  as  applicable  by  carefully  keeping  higher  order  terms  in 
the  velocity  components  or  in  the  perturbation  amplitudes,  and  noting  the 
smallness  of  their  contributions  in  the  ranges  of  magnitudes  we  are  con¬ 
sidering.  ) 

Thus,  we  have  three  equations  in  three  unknowns ,  u,  v,  and  p.  To 
solve  them,  we  first  note  that  (55)  is  satisfied  if  we  can  find  a  function, 
cp,  such  that 


(58) 
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and 

^2  +  ^ |  =  0  (59) 

bx  by 

The  function,  9,  is  called  a  potential  function.  Furthermore,  with 
these  substitutions,  (56)  and  (57)  become 


which  two  equations  are  consistent  and  lead  to 

cKp 

P  =  p0  +  p  3t  +  pgy 

Now  this  is  a  two- fluid  problem  with  a  horizontal  surface  of 
separation.  Let  the  surface  of  separation  be  denoted  by  the  equation 

y  =  A(t)  cos  kx  (60) 


We  label  the  upper  fluid  with  subscript  1 ,  and  the  lower  with  subscript 
2,  and  suppose  that  there  is  a  potential  function  for  each  fluid,  9^ 
and  9^ •  Likewise  there  is  a  pressure  solution  for  each  fluid 


Pi  =  po,l  +  Pi67  +  pi  "3tT 


P2  =  Po,2  +  p2&  +  p2  “St 


It  may  be  verified  that  the  following  potential  functions  are  solutions 
of  (59): 
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<P1  =  e  ky  f(t)  cos  kX 
92  =  "  eky  f(t)  cos  kx 


(62) 


(We  have  not  specified  the  boundary  conditions  used  in  finding  the  solu¬ 
tion  of  (59);  they  are  such  that  the  fluid  is  at  rest  at  y  =  ±  00  and 
that  all  features  have  the  same  periodicity  in  x  as  the  initial  inter¬ 
face  between  fluids.) 

We  may  now  obtain  the  solution  in  complete  form  by  applying  the 
matching  conditions  at  the  interface.  These  are: 

1.  The  interface  moves  with  the  fluid  velocity. 

2.  The  pressure  is  continuous  across  the  interface. 

The  first  of  these  can  be  approximated  in  this  low  amplitude  study  by 

v(y  -  0)  -  (£) 

N  y=o 

which  reduces  to 

i-fw 

The  second  interface  condition  expresses  the  equality  of  the  two  pres¬ 
sures  in  (6l)  at  y  =  0.  Thus 

Po,1  +  P1  A  cos  kx  +  ^  cos  kx^ 

/  df  \ 

=  Po,2  +  p2  Vg  A  cos  “  at  cos  **) 

Now  the  low  amplitude  stage  we  are  considering  is  supposed  to  be  but  a 


-190- 


PROBLEM  6.  INSTABILITY  OF  INTERFACE  BETWEEN  FHJIDS 

slight  perturbation  from  complete  equilibrium  at  zero  amplitude,  so 
that  the  zero-amplitude  pressures  must  Just  balance,  and  the  results  of 
applying  the  two  interface  conditions  all  reduce  to 


(65) 


Thus,  with  g  <  0  (i.e.,  the  acceleration  pointing  downward),  then  the  co¬ 
efficient  of  A  on  the  right  is  positive  when  the  upper  fluid  is  the  more 
dense  __  leading  to  exponentially  increasing  amplitude  —  and  negative 
when  the  upper  fluid  is  the  less  dense  —  leading  to  a  time-wise  oscil¬ 
lation  of  amplitude.  The  first  case  is  that  of  Taylor  instability. 

Combined  Helmholtz  and  Taylor  Instabilities 

Again  we  consider  the  case  of  low  amplitude  perturbations.  The 
situation  is  as  in  (a)  except  that  now  the  upper  fluid  moves  to  the 
right  with  velocity  u  =  U  relative  to  the  lower  one  (which  can  be  con¬ 
sidered  at  rest  without  loss  of  generality).  The  equations  must  be 
generalized  somewhat,  and  we  shall  employ  a  slightly  different  technique 

for  solving  them. 

The  mass  equation  remains  the  same  as  before,  (55),  and  we  satisfy 
it  again  with  upper  and  lower  potential  functions  and  q>2»  The  momen¬ 
tum  equations  for  the  upper  fluid,  however,  must  retain  the  one  compon¬ 
ent  of  transport  relative  to  the  rightward  motion,  so  that  the  equations 

now  are  written 

p||+pUli  +  M  =  0  1  (a) 

P^+pUH  +  ^'pg  =  °J 
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Thus  the  two  pressure  integrals,  analogous  to  (6l),  become 


/  3cp1  <*PA 

p1  =  Po,1  +  plVF  +  ~$t~  +  U  ~SxJ 


We  now  write  the  equation  of  the  interface  in  the  form 


cut  -  ikx 

y  -  yQ  e 


(66) 


and  take  for  the  potential  functions  (which  satisfy  (59)  and  proper  boun¬ 
dary  conditions) 


=  a  yQ  e 


-ky  +  d>t  -  ikx 


Ux 


q>2  =  P  yo  e' 


ky  +  cut  -  ikx 


(67) 


where  a,  (3,  and  co  are  constants  to  be  determined,  and  k  is  the  wave  num¬ 
ber  of  the  interface  disturbance.  Again  the  same  two  interface  condi¬ 
tions  must  apply,  and  there  must  also  in  this  manner  of  treatment  be  an 
expression  for  continuity  of  vertical  velocity  across  the  interface. 


Thus,  for  the  upper  fluid 


v(y  =  0) 


while  for  the  lower  fluid 

v(y  -  0)  =  (|) 

x  "y=o 

These  two  conditions  lead  to  the  determination  of  a  and  p 
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The  continuity-of -pressure  condition  shows  that  in  complete  equilibrium, 

*0,1  -  *1  ^  *  *0,2 

and  that  in  lowest  order  perturbation, 

Pl(g  +  ao>  -  i  kUa)  =  p2(g  +  Pa>) 
which  reduces,  finally,  to 

_  ilC[,P1„  *  V  I  &  fell +  (68) 

“  p,  +  p2  y  k  'P 2  +  PV  (o2  +  P, ) 

With  0  =  0,  the  result  is  the  same  as  before.  With  U  4  0,  there  are  two 
additional  contributions.  The  first  on  the  right  of  (68)  expresses  the 
mean  translation  velocity  of  the  waves.  The  second,  under  the  square 
root,  expresses  the  Helmholtz  instability  effect.  It  is  always  positive, 
hence  always  contributes  to  giving  to  a  real  part  corresponding  to  ex¬ 
ponential  growth  of  the  instability.  Even  in  the  case  that  g  =  0  and 
p  =  p2,  the  instability  remains,  in  which  case 

to  =  \  kU(i  ±1) 

The  results  can  be  extended  even  further  to  include  the  effects 
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of  surface  tension;  rather  than  repeat  the  analysis,  we  simply  present 
the  final  equation  for  the  exponential  growth  factor 


CD 


lkUP1  .  ,  /g  ^  ,  "Vi 

pi  +  p2  y  k  \p2  +  ( p2  +  p1  r 


kT 

p2  + p! 


(69) 


in  which  T  is  the  surface  tension.  It  is  seen,  as  could  have  been  ex¬ 
pected,  that  surface  tension  has  a  stabilizing  effect.  Likewise  it  can 
be  seen  that  there  is  a  most -unstable  wave  length,  obtained  by  maximiz¬ 
ing  the  second  term  on  the  right  of  (69)  —  provided  the  second  term  has 
a  real  maximum. 

The  effects  of  viscosity  are  much  more  difficult  to  include  in  full 
generality;  however,  in  some  cases  it  is  sufficiently  accurate  to  write 


CD 


ikUPl 


p1+p2 


U2P2P1  ^  kT 
(p2  +  Pl  )2  P2  +  Pl 


(70) 


where 


^  +  ^2 

Pi  +  p2 


and  1^,  are  the  coefficients  of  viscosity.  Discussions  of  the  vali¬ 
dity  of  various  approximation  procedures  have  been  given  in  the  refer¬ 
ences  mentioned  above.  Chandrasekhar  and  Hide  both  give  in  addition 


*  See,  for  example,  Horace  Lamb,  "Hydrodynamics,"  Sixth  Edition,  Dover 
Publications,  New  York,  1932,  Section  268,  p.  46l  . 

**Various  treatments  have  been  given  by  Bellman  and  Pennington,  Quart. 
Appl.  Math.,  _12,  151,  1954;  by  S.  Chandrasekhar,  Proc.  Cambridge  Phil. 
Soc.,  51_,  162,  1955;  and  by  R.  Hide,  Proc.  Cambridge  Phil.  Soc.,  51  > 
179,  1955. 
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discussions  of  the  effects  of  gradual  density  gradients,  and  later 
authors  have  solved  various  problems  concerning  stability  of  variable 
atmospheres. 

At  present  there  appears  to  exist  no  comprehensive  discussion  of 
the  effects  of  diffusion  on  interface  instability,  although  a  qualita¬ 
tive  argument  can  be  given  to  show  that  the  effects  are  likely  to  in¬ 
crease  stability. 


-195- 


INDEX 


Accelerated  viscous  flow,  146 
Acceleration,  35 
Active  degrees  of  freedom,  125 
Adiabat,  82 

Adiabatic  equation  of  state,  42 
Adiabatic  heating,  4l 
Aerodynamics,  7*84 
Angular  momentum,  49 
conservation  of,  50 
Angular  velocity,  44 
rate  of  change  of,  46 
Atmosphere,  equilibrium  state  of,  ^6 
sound  in  variable  density  of,  55 

Boltzmann  equation,  1  37 
Boundary  layer,  thickness  of,  l4l 

Chapman- Jouguet  condition,  1 78 
Characteristic  line,  59,63 
Characteristic  solution,  58,89 
Characterization  of  fluid,  9,15 
Chemical  energy,  1 76 
Chemical  reactions,  14 
Circulation,  47 

Coefficient,  of  heat  conduction,  122 
of  viscosity,  115,118 
Cohesive  fluids,  113 
Compression,  maximum  from  shock,  80 
monotonic,  82 
by  shock,  71 

Compression  wave  in  viscous  fluid,  148 
Conduction,  coefficient  of,  122 
of  heat,  10,12,44 
rate  of,  1 20 


Conservation,  of  angular  momentum,  50 
of  energy,  20,37,72,123,135,156,138 
laws  of,  9 

of  mass,  17,72,135,155 
of  momentum,  1 9, 36, 72,118,1 20, 1 34, 1 36, 1 38 
Conservative  form  of  differential 
equations,  28,  29 
Constant- state  region,  63 
Contact  discontinuity,  71 
Contiguity,  assumption  of,  15 
Continuous  fluid,  approximation  of,  11 
definition  of,  1 4 
Coordinate,  cylindrical  32 
Eulerian,  1 6 
Lagrangian,  1 6 
Coordinate  systems,  31,58 
Corner,  flow  around,  >8l 
Curved  shock,  71 
Cylindrical  coordinates,  32 
Cylindrical  shell,  expansion  of,  35 


due  to  viscosity,  1 43 
Damping  of  sound  waves,  156 
Decay  of  shock  pulse,  91 
Decay  rate,  shock  pressure,  102,105 
Decomposition  of  velocity  field,  51 
Deformation  tensor,  52 
Degrees  of  freedom,  124 
active,  125 
Density,  1 6 

Derivative,  along  fluid  motion,  21 
of  volume  integral,  26 
Detonation,  underdriven,  1 78 
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Detonation  theory,  1 76 
Diaphragm,  fractured,  1 65 
Diatomic  molecule,  124 
Differential  equations,  conservative 
form  of,  28,29 

Diffusion,  of  molecules,  113 
of  velocity,  115 
of  viscous  effect,  140 
Dilatation  pressure,  119,  124 
Dilatative  dissipation,  123 
Discontinuity,  71 
Dissipation,  dilatative,  123 
Distribution  function,  molecular,  130 
Disturbance,  small,  53 
Drag,  viscous,  1 1 4 

Energy,  1 6 

conservation  of,  20,37,72,123,135, 

136.138 

Energy  conservation  equation,  21,72, 

123.135.136.138 

Energy  equation,  general  derivation 
of,  120 

internal,  1 6, 40 
kinetic,  16 
potential,  38 

Energy  equation,  general  derivation 
of,  120 

Entropy,  1 0, 1 6, 39,42,43 
change  across  shock,  81 
Equation,  energy  conservation,  21 , 

72,123,135,136,138 

mass  conservation,  17,72,133,135 
momentum  conservation, 

20. 72. 118.120.1 34. 1 36. 1 38 
of  state,  9,34 

of  state,  adiabatic,  42 
of  state,  molecular  theory,  136 
of  state,  polytropic,  40,126,136 
Equations,  Eulerian  form,  10,28,30, 
41,84,87 

hydrodynamic,  molecular  derivation 

of,  133 

integral  form,  30 
inverted,  106 

Lagrangian  form,  10,21,26,29,30,93 
molecular  derivation,  1 33 
nonlinearity  of,  76 


Equations,  two  dimensional,  138 
viscous,  conducting,  138 
viscous,  two-dimensional,  1 38 
Equilibrium  state  of  atmosphere,  56 
Escape  speed,  61,67 
Eulerian  coordinate,  1 6 
Eulerian  form  of  equations,  10,28, 
30,41,84,87 

Expansion,  into  a  vacuum,  111 
of  cylindrical  shell,  33 
Explosions,  7 
Explosive,  176 
External  forces,  35 

Field,  local  decomposition  of 
velocity,  51 
Field  variables,  15*73 
First  coefficient  of  viscosity, 1 1 5,1 18 
Flow,  accelerated  viscous,  1 46 
incompressible,  188 
Prandtl-Meyer,  182 
Fluid/  characterisation  of,  9*15 
definition  of,  1 3 
isentropic,  42,46,53 
simple,  43,58 
Fluids,  cohesive,  113 
isothermal,  157 
Flux,  29 

Flux  tensor,  1 1 6, 1 18 
Force  due  to  viscosity,  142 
Forces,  external,  35 
nonlocal,  35 

Fractured  diaphragm,  1 65 

Gas,  ideal,  34,39 
isentropic,  58 

polytropic,  40,42,65,108,126 
simple,  58 
Gas  constant,  35 
Gas  viscosity,  113 
General  solution  of  inverted 
equations,  1 09 

General  solution  of  simple  fluid,  65 

Heat,  coefficient  of,  122 
Heat  conduction,  10,12,44 
rate  of,  120 

Helmholtz  instability,  187*191 
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Hodograph  transf ormation,  1 0 6 
Hugoniot  relation,  82 
Hydrodynamic  equations,  molecular 
derivation,  1 33 

Ideal  gas,  34,39 
Impulsive  motion  of  plate,  140 
Impulsive  piston  motion,  92 
Incompressible  flow,  188 
Instability,  Helmholtz,  1 87 y  191 
Taylor,  187,188 

Integral,  time  derivative  of,  2 6 
Integral  derivation  of  shock  rela¬ 
tions,  73 

Integral  form  of  equations,  30 
Integration,  along  fluid  motion,  18 
Interaction  of  shock  and  rarefaction, 
91,  102 

Internal  energy,  l6,40 
behind  shock,  80 
molecular  definition  of,  130 
Internal  pressure,  34 
Inverted  equations,  general  solution 
of,  109 

transformation  to,  1 06 
Inverted  hydrodynamic  equations,  106 
Isentropic  fluid,  42,46,53 
Isentropic  gas,  58 
Isothermal  fluid,  1 57 
Isothermal  rarefaction,  1 59 
Isothermal  shock,  1 60 
Isothermal  sound  speed,  155 

Jacobian,  26,27 

Kinetic  energy,  1 6 

Lagrangian  coordinate,  l6 
Lagrangian  form  of  equations,  10,21,26 
29,30,93 

Lagrangian  shock  speed,  75 
Laws,  conservation,  9 
Liouville  Equation,  132 
Liouville*s  Theorem,  131 

Mach  line,  1 84 

Mach  number,  78,86,164 

maximum  behind  shock,  80 


Mass,  conservation  of,  17,72,133,135 
Material  velocity,  1 6 
Molecular  derivation  of  equations,  1 33 
Molecular  diffusion,  1 1 3 
Molecular  distribution  function,  130 
Molecular  dynamics,  11,  13,  14,  127 
Molecular  relaxation,  1 26 
Molecular  velocity,  113 
mean,  130 

Molecule,  diatomic,  124 
Molecules,  collision  with  wall,  127 
diffusion  of,  1 1 3 
velocity  of,  1 1 3 

Momentum,  1 6 
angular,  49 

conservation  of,  19,3^,72,118,120, 

154,136,158 

Momentum  flux  tensor,  11 6, 1 1 8 

Nonlinearity  of  equations,  76 
Nonlocal  forces,  35 
Numerical  analysis,  7 
51 

Oblique  shock,  1 62 
One-Dimensional  flow,  pseudo,  17 
true,  17 

Piston,  arbitrary  motion,  65 
uniformly  accelerating,  68 
withdrawal  from  gas  of,  60,66,90 
Polytropic  gas,  40,42,65,108,126 
scaling,  86 

Polytropic  equation  of  state, 40, 

1 26, 1 36 

one-dimensional,  135 
Potential  energy,  38 
Potential  function,  1 89 
Prandtl-Meyer  Flow,  182 
Pressure,  1 6 
dilatation,  1 1 9 
total,  119,136 

Progressive-wave  solution,  57 
Propagation  speed,  55 

Radiation,  150 

Rankine-Hugoniot  relations,  72 
Rarefaction,  1 1 
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Rarefaction,  due  to  piston  withdrawal,  90 
interaction  of  shock  and,  91,102 
isothermal,  159 
overtaking  shock,  1 02 
Rates  of  strain,  117 
Reactive  shock,  1 7 6 
Relaxation,  molecular ,  1 26 
Reynolds  number,  1 1 8 , 1 42 
Rotation  tensor,  52 

Scaling,  84,85*118 
polytropic,  86 

Second  coefficient  of  viscosity,  118 
Second  law  of  thermodynamics,  59 
Self-similar  solution,  87,93 
Shell,  cylindrical  expansion  of,  33 
Shock,  10,11,70 

always  compressive,  8l 
compression  by,  71 
maximum  from,  80 
curved,  71 

formation  of,  67,69,70,75 

bitting  a  density  discontinuity,  1 69 

infinite  strength  of,  79,80 

into  fluid  at  rest,  77 

isothermal,  160 

oblique,  1 62 

production  of,  IX, 67,75 

reactive,  176 

thickness  of,  in  air,  152 

tube  theory  of,  1 65 

weak:,  1 03 

Shock  pressure  decay  rate,  102,105 
Shock  pulse,  decay  of,  91 
Shock  relations,  fundamental  derivation 
of,  72 

general,  74,75,77 
integral  derivation  of,  73 
polytropic  gas,  77,78 
special  cases,  77 
Shock  signals,  weak,  10,53 
Shock  speed,  71 
Lagrangian,  75 
Shock  structure,  12,70 
conduction  only,  1 53 
theory  of,  1 48 
viscosity  only,  151 
Similarity  solution,  general,  87 


Similarity  solution, 

for  shock  pulse,  95 
Similarity  transformation,  86 
Simple  fluid,  43,58 
definition  of,  1 4 
Simple  gas,  58 
Simple  isentropic  fluid,  46 
Simple  wave,  63 

general  solution  for,  64 
Small  disturbance,  55 
Smoothing,  viscous,  1 43 
Solution,  characteristic,  58,89 
example  of,  for  viscous  flow,  1 44 
general,  for  7  =  3,  110,111 
of  inverted  equations,  109 
progressive-wave,  57 
self- similar,  87,93 
for  simple  wave,  64 
Sound  in  variable  density 
atmosphere,  55 
Sound  speed,  10,53,54,86 
general  formula  for,  55 
isothermal,  155 
viscous  conducting  gas,  154 
Sound  waves,  damping  of,  1 56 
Specific  field  variables, 
definition  of,  1 6 
Specific  heat,  4l 
Specific  volume,  1 6 
Speedgraph  transformation,  106 
Speed  of  escape,  61,67 
Speed  of  propagation,  55 
Speed  of  shock,  71 
Lagrangian,  75 
Speed  of  sound,  10,53,54,86 
isothermal,  155 
in  viscous,  conducting  gas, 154 
State,  equation  of,  9,34 
Stress,  viscous,  12,13 
Stress  tensor,  119,136 
Structure  of  shock,  12,70 
theory  of,  1 48 

Summation  convention,  51,1 16 
Symmetry  of  flux  tensor,  1l6 

Taylor  instability,  1 87,1 88 
Temperature  1 6, 39,40 

viscosity  dependence  on,  114 
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Tensor  notation,  51 
Thermodynamics,  second  law  of,  39 
Thickness  of  boundary  layer,  1 4l 
Thickness  of  shock,  70 
in  air,  1 52 
theory  of,  1 48 
Total  pressure,  119,136 
Transformation,  Eulerian  to 
Lagrangian,  22,24 
Lagrangian  to  Eulerian,  23 
similarity,  86 
to  inverted  equations,  106 
Transformation  determinant,  108 

Vacuum,  expansion  into  a.  111 
Variable  density  atmosphere,  55 
Variables,  field,  15,73 
Velocity,  l6 
angular,  44 

field,  local  decomposition  of,  51 
material,  1 6 
mean  molecular,  130 
of  molecules,  1 1 3 
Viscosity,  10,12, 44 
artificial,  150 


coefficient  of,  115,118 
damping  due  to,  1 43 
force  due  to,  142 
gaseous,  1 1 3 
theory  of,  113,146 
Viscous  drag,  114 

Viscous  equations,  two-dimensional,  1 38 
Viscous  flow,  with  acceleration,  146 
example  of  solution,  1 44 
Viscous  fluid,  definition  of,  1 3 
Viscous  force,  1 42 
Viscous  smoothing,  1 43 
Viscous  stress,  12,13 
Volume  integral,  time  derivative  of,  26 
Vorticity,  45 

Wave,  simple,  63 
Wave  equation,  55*56 

Weak  shock,  1 03 

Weak  signals,  10,53 

Wedge,  flow  past,  162 

Work  done  on  fluid  element,  120,122 

Underdriven  detonation,  1 78 
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